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Summary 


The  research  performed  under  the  grant,  during  the  period  1  May  1994  through  30 
April  1995,  can  be  divided  into  three  main  topics;  (1)  regional  source  parameters,  seismic 
energy,  and  discrimination,  (2)  determining  surface- wave  magnitudes  from  regional 
Nevada  Test  Site  data,  and  (3)  isotopic  and  deviatoric  moment  inversion  of  regional 
surface  waves  from  Nevada  Test  Site  explosions:  implications  for  yield  estimation  and 
seismic  discrimination. 

In  section  1,  we  have  examined  broadband  waveforms  from  a  large  number  of 
NTS  explosions  and  earthquakes  throughout  the  southwestern  United  States  in  order  to 
characterize  seismic  sources.  Explosions  were  found  to  be  richer  in  coda  energy  than 
earthquakes.  Most  earthquakes  show  relatively  little  long-period  (T>4  sec)  coda  energy 
and  tend  to  be  richer  in  long-period  and  shear-wave  energy  than  explosions.  We  have 
developed  several  seismic  discriminants  based  on  these  observations  and  our  modeling 
experience.  One  promising  discrimant  is  the  ratio  of  short-period  vertical  component,  P- 
wavetrain  energy,  to  long-period  surface  wave  energy,  averaged  over  three  components. 
Explosions  tend  to  have  a  higher  ratio  than  do  earthquakes,  essentially  an  extension  of 
Mb:Ms.  Magnitude  threshold  for  this  discriminant  is  about  3.5.  Another  useful 
discriminant  is  based  on  the  total  broadband  energy  to  moment  ratio  where  explosions  are 
distinguished  by  their  stronger  energy  levels  relative  to  their  long-period  amplitudes. 

This  approach  requires  Green’s  functions,  a  source  estimator  program,  and  processes  all 
events  as  earthquakes.  For  this  method  to  be  effective  requires  the  calibration  of  the 
region  using  relatively  large  earthquakes,  M>5,  but  does  not  require  calibrations  of 
explosions. 

In  section  2,  we  re-examine  the  use  of  surface-wave  magnitudes  to  determine  the 
yield  of  underground  nuclear  explosions  and  the  associated  magnitude-yield  scaling 
relationship.  We  have  calculated  surface-wave  magnitudes  for  190  Nevada  Test  Site 
(NTS)  shots  using  regional  long-period  seismograms  from  a  combined  super-network  of 
55  north  American  stations.  Great  effort  went  towards  making  the  data  set 
comprehensive  and  diverse  in  terms  of  yield,  source  location  and  shot  medium  in  order  to 
determine  the  portability  of  surface-wave  magnitude  scales.  In  particular,  we  examine 
Pahute  Mesa,  Rainier  Mesa  and  Yucca  Flat  explosions  detonated  above  and  below  the 
water  table,  and  which  range  over  three  orders  of  magnitude  in  yield.  By  observation  we 
find  a  low-yield  measure  threshold  of  approximately  one  kiloton  (kt)  for  (assumedly) 
moderately  well-coupled  explosions  recorded  at  near-regional  (<500  km)  stations,  which 
have  little  microseismic  noise.  In  order  to  utilize  regional  surface  waves  (A<15  °)  for 
quantifying  sources  and  for  discrimination  purposes,  we  have  developed  related  methods 
for  determining  time-domain  surface-wave  magnitudes  and  scalar  moments  from  regional 
Rayleigh  waves.  Employing  regional  surface-wave  data  lowers  the  effective  magnitude 
threshold.  One  technique  employs  synthetic  seismograms  to  establish  a  relationship 
between  the  amplitude  of  the  regional  Airy  phase,  or  Rayleigh  pulse  of  the  Rayleigh 
wavetrain  and  an  associated  surface- wave  magnitude,  based  on  conventional  Ms 
determinations,  calculated  from  synthetic  seismograms  propagated  to  40°.  The  other 
method  uses  synthetic  seismograms  in  a  similar  fashion,  but  the  relationship  used  is  a 
more  straightforward  one  between  scalar  moment  and  peak  Rayleigh  wave  amplitude. 
Path  corrections  are  readily  implemented  to  both  methods.  The  inclusion  of  path 
corrections  decreases  the  Ms  variance  by  a  factor  of  two  and  affects  the  absolute  scaling 
relationship  by  up  to  a  factor  of  0.1  magnitude  units.  This  latter  effect  is  attributed  to  the 
particular  station  network  used  and  the  Green's  functions  used  to  obtain  the  40°  Ms 
values.  Using  a  generic  structure  for  the  distance  traveled  past  the  actual  source-receiver 
path  minimizes  the  difference  between  magnitudes  determined  with  and  without  path 
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corrections.  The  method  gives  stable  Ms  values  that  correlate  well  with  other  magnitude 
scale  values  over  a  range  of  three  orders  of  magnitude  in  source  yield.  Our  Ms  value 
scale  very  similarly  to  more  standard  teleseismic  Ms  values  from  other  studies,  although 
the  absolute  Ms  values  vary  by  ±0.5  magnitude  units  about  ours.  Such  differences  are 
due  in  part  to  the  choice  of  Ms  formula  used.  For  purposes  of  future  user  comparisons, 
we  give  conversion  values  to  the  previous  studies.  Our  most  refined  Ms  values  give  the 
relationship  Ms  =  1.00  x  logio  (yield)  +  B,  where  B  is  dependent  upon  source  region  and 
shot  medium.  This  yield  exponent  of  unity  holds  for  events  of  all  sizes  and  is  in  line  with 
A/s-yield  scaling  relations  found  by  other  studies.  When  events  are  grouped  with  respect 
to  source  region,  significantly  better  fits  to  these  individual-site  linear-regression  curves 
are  obtained  compared  to  the  fits  obtained  using  a  single,  all-inclusive  model.  This 
observation  implies  that  shot-site  parameters  and  source  structure  can  significantly  affect 
surface-wave-magnitude  measurements.  We  present  these  Ms  values  primarily  to 
augment  the  extensive  historical  analysis  of  explosion  data  based  on  surface- wave 
magnitudes  by  using  regional  data  to  increase  the  number  of  events  with  surface-wave 
magnitudes.  These  magnitudes  are  consistent  with  the  teleseismically  determined 
magnitudes  of  larger  events. 

In  section  3,  seismic  moments  of  Nevada  Test  Site  (NTS)  explosions  are 
determined  from  regional  surface  wave  spectra.  Two  methods  are  used.  In  one  the 
moment  is  solved  for  assuming  only  an  explosive  source,  or  average  scalar  moment;  in 
the  other  a  joint  inversion  for  an  isotropic  (explosive)  source  plus  a  constrained  double 
couple  moment  component  representing  tectonic  strain  release.  Although  the  general 
moment  tensor  solution  to  this  joint  inversion  problem  is  non-unique,  if  some 
assumptions  are  made  concerning  the  non-isotropic  moment  components,  then  the 
remaining  source  parameters  can  be  solved  by  a  linear  least-squares  inversion  scheme. 

We  examine  the  errors  in  determining  the  isotropic  moment  component  (M\)  by  this  latter 
method  of  constrained  linear  inversion  solutions  in  a  canonical  study  using  a  theoretical 
network  of  long-period  (6-60  sec.)  surface  wave  data.  The  network  azimuthal  coverage 
was  chosen  to  represent  that  of  a  long-period  North  American  super-network  of  55 
stations  used  for  the  actual  NTS  events.  We  compare  these  errors  in  moment  estimate  to 
those  obtained  from  surface  wave  magnitude  (M$)  measurements  for  the  same  surface 
wave  observations.  For  a  ratio  of  M(expi/M(eq)  less  than  1.0  we  find  that  the  inverted  M\ 
solution  is  a  much  better  estimate  of  the  actual  isotropic  moment  than  either  Ms  or  Mo, 
and  the  standard  deviation  in  this  estimate  is  substantially  less  than  that  using  the  other 
two  methods  for  the  great  majority  of  isotropic  source  +  double  couple  sources.  Even 
when  the  inversion  constraints  are  off  in  dip  and  rake  each  by  30°,  the  mis-estimate  of  the 
isotropic  moment  is  less  than  35  percent  of  the  actual  value.  In  the  case  of  a  vertical 
strike-slip  fault,  the  inverted  isotropic  moment  solution  which  assumes  this  fault 
orientation  is  exact  to  three  figures,  whereas  Ms  and  Mo  under-estimate  the  moment  by 
45  percent  and  32  percent,  respectively  because  of  uneven  azimuthal  coverage.  This 
moment  tensor  inversion  method  is  applied  to  determine  the  isotropic  source  for  109  NTS 
underground  explosions  using  vertical  and  tangential  component  surface  wave  data  from 
this  regional  network.  We  also  calculate  M$  and  Mo  for  these  same  events  and  compare 
the  results.  Isotropic  source  errors  are  smallest  using  the  spectral  domain  inversion 
method.  However,  this  spectral  domain  method  cannot  attain  as  low  a  magnitude 
threshold  as  the  time  domain  moment  or  Ms  method.  The  extensive  moment  data  set 
analyzed  here  were  combined  with  larger  yield  explosions  from  prior  moment  studies  to 
create  a  comprehensive  data  set  with  which  to  obtain  conclusive,  well-constrained  long- 
period  explosion  source  scaling  relationships  at  the  separate  NTS  sub-sites. 


SECTION  1 


Regional  Source  Parameters,  Seismic  Energy, 

and  Discrimination 


REGIONAL  SOURCE  PARAMETERS,  SEISMIC  ENERGY,  AND 
DISCRIMINATION 


DON  V.  HELMBERGER 
Seismological  Laboratory  252-21 
California  Institute  of  Technology 
Pasadena,  CA91125 

and 


BRAD  WOODS 
Woodward  Clyde  Consultants 
566 El  Dorado  Street,  Suite  100 
Pasadena,  CA  91101 


1.  Abstract 

We  have  examined  broadband  waveforms  from  a  large  number  of  NTS  explosions  and 
earthquakes  throughout  the  southwestern  United  States  in  order  to  characterize  seismic 
sources.  Explosions  were  found  to  be  richer  in  coda  energy  than  earthquakes.  Most 
earthquakes  show  relatively  little  long-period  (T>4  sec)  coda  energy  and  tend  to  be  richer 
in  long-period  and  shear-wave  energy  than  explosions.  We  have  developed  several 
seismic  discriminants  based  on  these  observations  and  our  modeling  experience.  One 
promising  discriminant  is  the  ratio  of  short-period  vertical  component,  P-wavetrain 
energy,  to  long-period  surface  wave  energy,  averaged  over  three  components.  Explosions 
tend  to  have  a  higher  ratio  than  do  earthquakes,  essentially  an  extension  of  mb:Ms. 
Magnitude  threshold  for  this  discriminant  is  about  3.5.  Another  useful  discriminant  is 
based  on  the  total  broadband  energy  to  moment  ratio  where  explosions  are  distinguished 
by  their  stronger  energy  levels  relative  to  their  long-period  amplitudes.  This  approach 
requires  Green’s  functions,  a  source  estimator  program,  and  processes  all  events  as 
earthquakes.  For  this  method  to  be  effective  requires  the  calibration  of  the  region  using 
relatively  large  earthquakes,  M>5,  but  does  not  require  calibrations  of  explosions. 


2.  Introduction 

A  number  of  broadband  arrays  have  been  introduced  in  recent  years.  One  such  array, 
TERRAscope,  is  presently  being  installed  in  southern  California,  see  figure  1.  These 
stations  are  a  part  of  the  global  IRIS  (International  Research  Institution  for  Seismology) 
network.  Presumably  the  IRIS  network,  in  conjunction  with  short-period  arrays,  will 
provide  some  of  the  essential  data  necessary  for  worldwide  monitoring  of  seismic  activ¬ 
ity.  Unfortunately,  station  spacing  in  remote  areas  is  rather  sparse.  Thus,  we  may  have 
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to  rely  on  a  single  station  to  characterize  events  and  to  distinguish  an  explosion  from  an 
earthquake.  This  task  will  be  difficult,  but  may  be  possible  in  regions  that  have  an  abun¬ 
dance  of  earthquakes  that  can  be  used  to  calibrate  regional  paths.  We  envision  an  envi¬ 
ronment  not  unlike  that  of  western  United  States  and  thus  we  can  use  the  population  of 
NTS  events  and  the  natural  seismicity  to  construct  and  test  such  a  methodology.  Before 
discussing  regional  calibration  and  the  development  of  energy  discrimination  techniques, 
we  will  first  briefly  review  the  observational  differences  between  explosions  and  earth¬ 
quakes. 


Figure  1.  Map  of  southwutern  United  States  displaying  the  locations  of  a  number  of  recent  earthquakes  and 
NTS  explosions;  also  included  is  the  TEkRAscope  array  as  existed  in  early  1993. 


Most  of  the  useful  regional  discriminants  for  populations  of  explosions  and 
earthquakes  in  this  region  have  been  discussed  by  Taylor  et  al.,  (1989).  Their  results 
suggest  that  mb:M$  works  very  well  for  well-calibrated  paths  but  they  had  difficulty  in 
determining  Ms  for  explosions  smaller  than  about  mb=4.  In  contrast,  they  report  Ms  for 
earthquakes  with  mb’s  as  small  as  2.5.  The  characteristics  of  events  as  described  above 
are  quite  compatible  with  TERRA  scope  observations,  as  displayed  in  figure  2,  where 
seismograms  of  the  Kearsarge  explosion  (mb*5.6)  is  compared  with  the  Skull  Mountain 
earthquake  (mb=5.7).  Simulations  of  the  broadband  data  appropriate  for  various  instru¬ 
ments  are  included  since  this  type  of  data  normally  used  in  defining  mb  and  Ml 
(WASP),  and  MS(LP).  Note  that  the  ratio  of  peak  short-period  to  long-period  amplitudes 
(averaged  over  the  components)  is  an  order  of  magnitude  larger  for  the  explosion. 
Because  the  paths  are  nearly  identical,  this  difference  must  be  caused  either  by  the  source 
excitation  and/or  epicenter  depth.  Thus,  it  appears  that  mb:Ms  can  be  extended  to  small 
events  if  the  Rayleigh  waves  can  be  detected. 
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re  2.  Conparison  of  the  broadband  observation*  and  assorted  simulated  instrumental  responses  of  a  150  kl 
explosion  (JVE)  and  a  recent  earthquake  at  Skull  Mountain  as  recorded  at  PAS  (A=350  km).  Note  that  the 
paths  are  nearly  identical  but  the  motions  are  quite  distinct. 
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Hgure  i.  A  small  Yucca  Flat  explosion,  Floydada  (Ml=40),  recorded  at  TERRAscope  stations  and  simulated 
Press-Ewing  Oong-period)  and  Wood- Anderson  (short-period),  after  Woods  et  al.  (1993). 
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In  figure  3,  we  display  the  TERRAscope  data  from  a  small  NTS  explosion 
where  the  Rayleigh  waves  are  quite  apparent  even  at  long-pen ods.  Note  that  the  surface 
wave  amplitudes  indicate  that  these  signals  would  not  be  discernible  on  the  actual  analog 
instruments.  Ml  for  this  event  is  4.0  and  its  log  moment  is  14.20,  Woods  et  al.,  (1993). 
Assuming  it  is  a  shallow  explosion  above  the  water  table,  the  yield  can  be  inferred  to  be 
less  than  10  kt  from  the  moment-yield  scaling  relationships  determined  for  NTS  by 
Woods  and  Harkrider  (1994).  Were  it  detonated  in  hard  rock  below  the  water  table  it 
would  correspond  to  a  two  kiloton  explosion.  We  estimate  that  were  this  event  2.5  times 
lower  in  yield,  it  would  still  be  possible  to  obtain  its  moment,  yielding  a  magnitude 
threshold  of  about  3.5.  These  observations  are  quite  typical  of  small  NTS  events  where 
only  the  fundamental  Rayleigh  wave  Airy  phase  is  detectable  and  reasonably  predictable 
across  the  array. 

Thus,  it  appears  that  regional  seismograms  from  explosions  are  indeed  higher 
frequency  than  earthquakes  with  comparable  Rayleigh  wave  excitation.  This  feature  is 
explored  by  Patton  and  Walter  (1993),  in  a  study  of  well-calibrated  mt>:Ms  and  the 
mb:Mo  discriminants.  Their  results  produced  a  clear  separation  as  does  the  Ml:Mo 
discriminant  proposed  by  Woods  et  al.,  (1993).  The  latter  results  are  displayed  in  figure 
4.  Note  there  is  a  significant  separation  of  earthquakes  and  explosions  with  no  real 
overlap  over  all  scales.  However,  many  of  the  M0*s  for  these  explosions  were 
determined  by  assuming  shallow  isotropic  excitation  (Woods  et  al.,  1993),  and  thus  while 
this  approach  demonstrates  that  the  high-frequency  vs.  low-frequency  source 
characterizations  remain  valid  to  small  magnitudes,  it  may  have  limited  usefulness  as  a 
direct  discriminant 


Figure  4.  Plot  of  Ml  vs.  M0  for  a  population  of  explosions  and  earthquakes,  after  Woods  et  al.  (1993). 
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Another  difficulty  with  the  Ml:Mo  discriminant  is  in  the  development  of  a 
physical  basis.  While  the  Ml  measure  is  easily  simulated  for  earthquakes  it  proves  prob¬ 
lematical  for  explosions.  For  example,  in  examining  figures  2  and  3,  we  find  that  the 
peak  short-period  amplitude  usually  occurs  on  the  tangential  component.  This  is  difficult 
to  explain  with  the  conventional  symmetric  RDP  (Reduced  Displacement  Potential) 
formalism  and  require  some  type  of  mode  conversion.  Presumably,  the  large  amount  of 
local  Rayleigh  wave  energy  released  by  the  source  into  the  slow-velocity -source ;  region 
gets  scattered  into  the  crustal  wave  guide  (e.g..  Stead  and  Helmberger,  1988).  We  find 
the  broadband  (BB)  records  from  earthquakes  occurring  in  the  normal  seismogenic 
depths  of  4  to  15  kms  to  be  relatively  simple  as  displayed  in  figure  2.  Thus,  we  will 
explore  the  possibility  of  using  the  ratio  of  accumulated  energy  to  surface  wave 
magnitude  or  Mp  as  a  discriminant 

To  pursue  this  approach  we  will  assume  all  events  are  earthquakes  with  respect 
to  estimating  source  parameters.  Explosions  are  then  distinguished  by  their  excess  short- 
period  energy.  However,  to  obtain  meaningful  estimates  of  source  parametersfrom  re¬ 
gional  waveforms  does  require  a  crustal  model  and  useful  Green’s  functions.  These  can 
be  obtained  by  modeling  moderate  sized  earthquakes. 
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Figure  5.  Display  of  velocity  vs.  depth  models  for  a  southern  California  model  (SC)  and  a  basin-and-range 

model  (PB). 
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3.  Estimation  Of  Earthquake  Parameters 

Essentially,  three  new  methods  have  been  developed  to  take  advantage  of  the  new  in¬ 
strumentation:  the  CMT  (Centroid  Moment  Tensor)  solution  at  long-periods  (Ritsema 
and  Lay,  1993),  inversion  of  long-period  body  waveforms  (Dreger  and  Helmberger, 
1993),  and  a  broadband  cut  and  paste  method  by  Zhao  and  Helmberger  (1994). 

For  events  larger  than  Mw>5,  it  is  possible  to  invert  surface  wave  records  for 
periods  greater  than  SO  seconds  assuming  the  PREM  model  (Dziewonski  and  Anderson, 
1981)  by  employing  a  CMT  procedure.  At  shorter  periods  the  surface  waves  show  re¬ 
gional  variation,  and  corresponding  regionalized  models  are  required  (see  for  example, 
Patton  and  Zandt,  1991;  Tbio  and  Kanamori,  1992). 

The  second  method  uses  the  relative  strengths  of  the  observed  bodywaves  com¬ 
pared  with  synthetics  to  determine  mechanisms,  moment,  and  depth.  Often,  only  one  sta¬ 
tion  is  sufficient  to  fix  the  source  parameters,  since  S  (SV  and  SH)  and  sS  (SV  and  SH) 
are  strongly  dependent  upon  source  orientation.  Cycling  through  source  depths  the 
proper  timing  between  P  and  pP,  S  and  sS,  etc.  allows  accurate  depth  estimates.  This  ap¬ 
proach  works  best  at  periods  greater  than  a  few  seconds  and  therefore  we  usually  work 
with  long-period  bandpassed  records.  The  Southern  California  model  (SC),  Dreger  and 
Helmberger  (1991),  displayed  in  figure  5  works  well  in  terms  of  waveform  matching 
throughout  the  entire  region,  as  reported  by  Dreger  and  Helmberger  (1993). 

The  third  approach  uses  a  direct  grid  search  for  the  fault  parameters  (strike  (0), 
dip  (5),  rake  (X.)).  This  method  matches  complete  broadband  observed  seismograms 
against  synthetics  over  discrete  phases  so  that  timing  shifts  between  particular  wave 
groups  are  allowed.  That  is,  in  matching  a  synthetic  seismogram  to  the  observed  record, 
we  may  allow  the  Rayleigh  wave  to  be  shifted  relative  to  the  Pnj  wavetrain.  This  allows 
a  better  correlation,  thus  the  name  cut-and-paste  method.  This  feature  desensitizes  the 
effect  of  the  crustal  model  used  in  generating  the  synthetics  and  allows  stable  estimates 
of  the  source  parameters  with  imperfect  Green’s  functions.  We  demonstrate  this  conclu¬ 
sion  by  generating  fault  parameters  for  a  number  of  regional  events  using  two  strongly 
contrasting  crustal  models,  the  SC  model  and  a  basin  and  range  model  (PB)  by  Priestley 
and  Brane  (1978),  displayed  in  figure  5.  The  source  parameter  determinations  are  given 
in  Table  1  for  three  large  events  in  the  region:  the  Utah  event,  the  Little  Skull  Mountain 
event,  and  the  Arizona  event  Paths  connecting  these  events  to  the  TERRAscope  stations 
provide  a  good  sample  of  the  propagational  features  of  the  region,  see  figure  1 . 

A  comparison  of  the  waveform  fits  assuming  the  PB  model  is  displayed  in 
figures  6  and  7.  The  numbers  above  each  trace  indicates  the  peak  amplitude  and  the 
moment  estimate  comes  from  a  least  square  fit  where  the  individual  amplitude  compari¬ 
son  indicates  the  relative  contribution  of  that  trace  to  the  average  moment.  The  relative 
shifts  of  the  various  phases  is  discussed  in  Zhao  and  Helmberger  (1994). 


4.  Estimation  of  Seismic  Energy  and  Discriminants 

With  the  advent  of  broadband  instrumentation,  it  is  possible  to  make  some  useful 
estimates  of  the  energy  levels  of  sources.  This  is  especially  attractive  at  local  and 
regional  distances  before  the  strong  attenuation  of  the  upper  mantle  has  snipped  away  the 
higher  frequencies.  However,  we  must  be  able  to  correct  for  the  strong  propagational 
effects  produced  by  the  crust  before  obtaining  accurate  energy  estimates. 
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Figure  6.  Comparison  of  broadband  observations  (Utah  event)  and  corresponding  synthetics  (PB)  as 
determined  in  obtaining  the  best  fitting  mechanism.  The  numbers  indicate  the  peak  amplitudes  in  cm.  The 
numbers  below  the  stations  indicate  the  distance  from  that  station  to  the  event  in  kms. 


Figure  7.  Comparison  of  synthetics  vs.  observations  for  the  Utah  event  where  the  time  history  has  been 
adjusted  such  that  the  energy  ratio  between  the  synthetic  (WASP/LP)  in  the  Pnl  window  is  equal  to  that 

observed  as  averaged  over  the  5  stations. 
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We  proceed  by  making  some  useful  definitions  in  terms  of  a  source  time  history 
needed  to  make  definitive  broadband  moment  estimates.  Because  energy  depends  on 
velocity,  we  must  be  careful  in  defining  the  time  history,  especially  at  the  highest  fre¬ 
quency.  The  approach  followed  here  is  to  base  this  estimate  on  the  Pnl  window  which 
we  think  is  the  least  contaminated  by  the  complex  surface  layer.  Since  the  WASP 
records  are  difficult  to  match  in  waveform,  we  choose  to  base  the  time  history  on  that  tri¬ 
angle,  5ti,  which  best  predicts  the  energy  ratio  of  synthetic  (WASP/LP)  to  that  of  the 
observations.  The  broadband  moment  (Mg)  is  defined  as  the  best  fitting  synthetic  to  the 
observed  data  assuming  this  time  history  with  the  orientation  parameters  determined  by 
the  long-period  fit 

In  the  same  fashion  as  Mb  we  define  the  energy  moment.  Me,  to  be  the  ratio  of 
the  total  integrated  broadband  energy  (3  components)  divided  by  the  corresponding  inte¬ 
grated  Green’s  functions.  To  be  more  specific  wc  will  write  down  the  explicit  expres¬ 
sions  for  the  tangential  component,  we  define 


M  =  Vobs(t)l  Vg(t ) 


where 

M=  moment 

Vg  (t)  =  synthetic  for  a  particular  far-field  source  history,  D  (t) 

Vobs(0  =  observed  record 

We  define  M=  M0  when  these  amplitude  comparisons  are  performed  at  long-periods, 

and  Mb  when  performed  with  the  D(t)  fixed  by  the  (SPZ/LP)  ratio.  The  energy 
strength  is  defined  by 


Me  = 


ki  w<*  /  k2 

.0  o 


il/2 


(t)dt 


(2) 


where  T  is  the  length  of  the  records.  This  same  procedure  is  applied  to  all  components 
and  averaged  to  define  the  M’s  for  a  particular  event.  If  the  synthetics  fit  the  observed 
data  exactly,  we  would  obtain  Me  equal  to  Mb  or  a  ratio  (Me/Mb)  of  1 . 

Applying  this  formalism  to  the  Little  Skull  Mountain  main  event  we  obtain 
broadband  fits  nearly  identical  to  those  in  figure  6.  We  obtain  a  time  history  given  by  8ti 
—  o.3  secs.  The  corresponding  short-period  comparisons  of  synthetics  with  observations 
are  displayed  in  figure  8.  These  comparisons  are  quite  good  and  are  typical  of  results 
from  other  events,  see  Zhao  and  Helmberger  (1994). 
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Figure  8 .  Con*>arison  of  predicted  short-period  synthetics  vs.  observations  for  the  Utah  event  assuming  the 

broadband  time  history. 
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Figure  9.  Plot  of  accumulated  energy  for  all  three  components  (broadband)  vs.  distance  for  the  three 

calibration  events. 
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We  do  not  want  definitions  of  source  properties  such  as  moment  and  energy  to 
depend  on  the  range  or  receiver.  Thus,  we  can  check  the  usefulness  of  our  definitions  by 
plotting  (Me/Mb)  as  a  function  of  distance  and  model  as  in  figure  9.  Most  of  the  points 
fall  between  .5  and  1.5.  This  plot  shows  no  obvious  distance  dependence  suggesting  that 
our  path  corrections  are  adequate.  However,  there  is  a  slight  baseline  shift  with  the  PB 
model  yielding  slightly  higher  values.  But  in  general,  the  various  moments  and  source 
parameter  appear  to  be  quite  independent  of  model  if  we  treat  deep  sources,  i.e.,  d>5  km. 
Note  that  the  more  detailed  the  upper  portion  of  the  model  becomes,  the  more  likely  it 
will  become  path  dependent.  If  we  want  to  use  the  same  model  for  a  large  region,  we 
want  to  keep  the  model  simple  and  restrict  the  source  depths  accordingly. 

With  this  brief  review  of  source  estimation,  we  return  to  the  observations  dis¬ 
played  in  figure  2.  If  we  simply  compare  the  top-traces,  we  see  that  while  the  peak 
short-period  (SP)  amplitudes  are  similar,  the  long-period  outputs  are  an  order-of- 
magnitudc  different.  This  comparison  is  typical  (Woods  et  al.,  1993).  Also,  note  that  the 
explosion  data  contains  many  more  arrivals  or  energy  than  do  the  earthquake  traces. 
Therefore,  if  we  simply  compare  the  ratio  of  (Me/Mb),  we  should  distinguish  the  two 
types  of  events.  But  to  do  this,  we  must  first  obtain  estimates  of  Mo,  Mb,  and  Me  from 
explosions.  We  do  this  by  assuming  all  events  are  earthquakes  or  double-couple's.  An 
example  calculation  is  given  in  figure  10. 
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Figure  10.  Comparison  of  data  (small  explosion)  with  synthetics  (assumed  to  be  double-couple).  In  this  case 
the  code  found  essentially  a  strike-slip  solution.  The  depth  of  5  km  was  determined  by  the  best  fitting  solution. 
Note  that  there  are  many  more  scattered  arrivals  in  the  data  that  are  not  in  the  nature  earthquake. 
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In  this  match  of  synthetics  to  observations,  we  have  used  the  (WASP/WALP) 
ratio  to  fit  the  time  history  because  of  the  noise  in  the  LP  bandpass.  A  value  of  (.1,-1) 
triangle  was  obtained.  While  the  short-period  details  are  not  well  explained,  the  overall 
estimate  of  long-period  waveform  fits  is  reasonable.  We  obtain  a  moment  of  1.1  x  102 
dyne-cm  and  a  source  orientation  of  (220°,  30°,  115°)  for  strike,  dip,  and  rake.  The 
source  depth  search  preferred  the  depth  of  5  km  which  is  the  shallowest  depth  allowed. 
Woods  et  al.,  (1993)  obtained  a  M0  =  3  x  1021  dyne-cm  for  this  event  or  about  4  times 
smaller  than  the  above  estimate.  This  difference  is  expected  because  of  the  relative 
strengths  of  Rayleigh  wave  Green’s  functions  for  two  reasons.  First,  since  the  excitation 
of  Rayleigh  waves  per  unit  of  moment  is  stronger  at  the  shallower  depth,  we  can 
understand  why  a  larger  moment  is  needed  to  fit  the  data  assuming  a  depth  of  5  Ion. 
Second,  since  the  radiation  pattern  for  an  earthquake  is  always  less  than  for  an  explosion, 
we  again  require  an  increased  moment  to  compensate.  The  Mb  moment  is  obtained  by 
matching  the  amplitudes  in  the  upper  two  plots  yielding  .5  x  1022  dyne-cm,  which  is 
lower  than  M0.  This  is  caused  by  the  short-period  spikes  riding  on  top  of  the  observed 
Rayleigh  waves;  something  which  does  not  occur  in  the  synthetics  but  holds  true  for 
most  observations  of  explosions.  Thus,  the  Mb  measure  from  explosions  is  not  a  very 
good  measure  of  source  strength  since  it  appears  to  be  affected  by  these  short-penod 
spikes.  The  Me  measure  is  also  strongly  affected  by  these  spikes  and  decreased 
accordingly.  However,  the  extra  arrivals  occurring  in  the  observations  increase  Me, 
yielding  Me  =  1  2  x  1022  dyne-cm  or  a  ratio  of  Me/Mb  =  2.3.  Earthquakes  yield  values 
near  unity  as  discussed  above,  thus  it  is  an  effective  depth  discriminant 


Figure  11.  Plot  of  (Me/Mb)  vs.  Mq  for  a  small  population  of  earthquakes  and  NTS  explosions  displaying 

good  separation. 
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Figure  11  shows  the  results  for  NTS  explosions  and  southwestern  U.  S.  earth¬ 
quakes.  Each  point  represents  one  event.  The  chemical  NPE  (non-proliferation 
explosion)  also  has  been  added  to  the  data  set  This  event  and  a  normal  Rainier  shot 
(denoted  with  an  “R”)  lie  at  the  lower  limit  of  the  explosion  population,  but  they  are  still 
separated  from  the  earthquake  population. 

A  more  empirical  approach  is  to  use  the  short-period:long-period  (SP:LP)  ratio. 
The  source  properties  that  we  want  to  quantify  are  the  short-period  (1  Hz)  P-wave  and 
long-period  (0.14  to  0.05  Hz)  surface  wave  energy  levels,  the  ratio  of  which  is  used  as 
the  domination  criterion.  The  short-period  bandpass  is  the  same  used  to  measure  tele- 
seismic  P-wave  amplitudes  for  the  discriminant.  The  long-period  bandpass  rep¬ 
resents  the  predominant  frequency  range  of  the  fundamental-mode  Airy  Phase  at  re¬ 
gional  distances  (Alewine,  1972).  This  short-period  vs.  long-period  energy  ratio 

(EsP:LP)  ^  defined  as* 


with  the  summation  being  for  the  three  components  and  fj  representing  the  windowing 
fimpc  determined  from  travel  path  length  and  the  wave  train  of  interest;  tpn  corresponds 
to  the  time  before  the  onset  of  the  P-wave  and  the  time  prior  to  the  S-wave  onset 
time,  and  f/  and  t2  bracket  the  time  window  of  the  fundamental  Rayleigh  and  Love 
waves.  Vsp  and  v/p  are  the  short-period  and  long-period  ground  bandpass  velocities, 
respectively.  vsp  is  obtained  by  convolving  the  broadband  velocity  record  convolved 
with  a  Wood-Anderson  (WA)  short-period  instrument  and  v/p  is  the  broadband  velocity 
record  convolved  with  a  long-period  instrument  (PE).  The  velocities  are  squared  in  order 
to  obtain  units  of  energy;  the  factor  of  m/2  in  the  numerator  and  denominator,  where  m  is 

the  unit  mass  of  the  particle  of  motion,  cancel  out 

Figure  12  displays  the  log  SP:LP  integrated  energy  ratio  vs.  distance  for  all 
data;  each  data  point  represents  one  event-station  pair.  Crosses  represent  earthquakes, 
circles  signify  nuclear  explosions,  and  stars  are  data  points  for  the  chemical  kt  test  The 
explosions  tend  to  have  higher  SP:LP  integrated  energy  ratios  than  do  the  earthquakes  at 
all  distance  ranges.  Although  there  isn’t  complete  separation  of  the  two  populations,  the 
portion  of  the  earthquake  population  which  overlaps  with  the  explosion  population  is 
small  (approximately  10  percent). 

5.  Discussion 

In  order  to  better  appreciate  the  robustness  of  these  energy  measurements  and  their 
application  to  particularly  small  events,  we  will  show  a  suite  of  regional  waveforms  and 
their  associated  integrated  energy  curves  for  small  explosions  from  the  three  NTS 
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subsites  and  for  earthquakes  near  NTS,  namely  the  Little  Skull  Mountain  sequence.  We 
will  also  briefly  discuss  some  2-D  scattering  effects  produced  by  surface  geology. 


SP-LP(3-comp)  Ratio  vs.  Distance,  North  America 


1000  1100  1200  1300  1400 


Figure  12.  Plot  of  the  short  periodslong  period  discriminant  for  a  population  of  events,  after  Woods  and 

Helmberger(1994). 

Figure  13  displays  vertical  component  records  of  NTS  shots  recorded  at  GSC 
(top  four  rows)  and  ISA  (bottom  four  rows).  The  first  column  is  the  broadband 
displacement,  followed  by  a  convolution  with  the  WA  instrument  on  the  right  The  test 
site  or  name  of  the  event  is  to  the  left  of  each  row.  At  both  these  stations,  long-period 
fundamental-mode  Rayleigh  wave  energy  is  evident  for  all  four  events,  with  T  <  7  sec. 
dispersed  coda  waves  having  the  largest  amplitudes.  In  making  comparisons  between 
events,  it  should  be  noted  that  the  Pahute  shot  is  larger  than  the  others,  so  that  it  has  a 
much  better  signal  to  noise  ratio.  On  the  broadband  recordings  there  are  no 
distinguishing  waveform  characteristics  between  the  three  test  sites.  This  holds  true  for 
all  TERRAscope  records,  Woods  and  Helmberger  (1994).  On  inspection  of  the  WA 
records,  certain  patterns  emerge.  Recordings  of  the  Pahute  shot  have  a  strong,  prominent 
P-wavetrain  followed  by  relatively  small  short-period  coda  compared  to  the  other  test- 
site  shots.  The  Yucca  shot  shows  a  great  deal  of  what  appears  to  be  scattered  shear-wave 
energy  at  both  stations,  which  equals  or  exceeds  the  peak  P-wave  amplitude.  The  Ranier 
shots  seem  to  be  an  intermediate  case,  for  which  the  shear-wave  energy  may  nearly  rival, 
but  not  exceed  the  peak  P-wave  amplitudes.  In  all  cases  the  shear-wavetrain  does  not 
have  a  sharp  onset,  but  rather  is  a  dispersive  wavetrain. 
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Figure  13.  Comparison  of  NTS  events  at  GSC  (A  -  330  Ian).  Note  that  the  seismograms  of  the  chemical 
explosion  fit  nicely  into  the  other  underground  explosion  population.  Pahute  events  appear  distinct  in  their 

simplicity  at  this  azimuth. 


These  NTS  records  are  in  sharp  contrast  to  natural  earthquake  records  as 
displayed  in  figure  14  where  the  onset  of  shear  waves  is  quite  clear.  These  recordings 
are  at  fairly  near  distances  (200  <  D  <280  km),  so  that  propagation  effects  are 
minimized.  However,  these  characteristics  persist  to  larger  ranges.  Woods  and 
Helmberger  (1994).  These  waveform  characteristics  are  displayed  in  the  integrated 
energy  curves  as  well. 

The  panel  on  the  left  in  figure  15  displays  the  short-period,  vertical  component, 
integrated  energy  curves  vs.  time  for  explosions  recorded  at  the  four  TERRAscope 
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stations.  All  curves  are  normalized  to  unity,  with  the  actual  integrated  energy  value 
given  to  the  right  in  the  legend.  We  will  refer  to  these  plots  as  E(t),  or  accumulated 
energy  up  to  time  t.  For  records  with  a  prominent  P-wavetrain  arrival,  the  curves 
resemble  step  functions.  This  is  particularly  true  of  the  Pahute  shot  and  expected  on 
theoretical  grounds,  i.e.,  RDP  source.  For  the  Yucca  shot,  the  P-wavetrain  energy 
comprises  less  than  half  of  the  short-period  energy.  The  Rainier  nuclear  explosion  and 
chemical  blast  are  intermediate  in  shape  to  the  Pahute  and  Yucca  energy  curves.  At  the 
more  distant  stations  (PAS  and  PFO),  the  Rainier  energy  curves  more  closely  resemble 
those  of  Yucca  than  Pahute,  whereas  at  the  closer  stations,  particularly  ISA,  the  opposite 
is  true. 
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Figure  14.  Comparison  of  the  Skull  Mountain  events  at  GSC  and  ISA. 
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Figure  16.  Upper  panel  displays  a  2D  cross-section  from  Imperial  Valley  (1  eft)  to  Pasadena  (right),  distance  is 
about  260  km.  Lower  panel  displays  the  comparison  of  a  relatively  simple  event  (AN)  compared  to  the  more 
complex-looking  event  which  occurred  in  the  basin  (IV).  These  are  long-period  comparisons  with  synthetics 

predicted  from  the  above  model. 

The  panel  on  the  right  in  figure  16  displays  analogous  plots  of  integrated  energy 
curves  for  earthquakes  from  the  Little  Skull  Mountain  sequence.  Here  the  onset  of  die  S- 
wave  energy  is  pronounced  and  sharp.  For  explosions,  in  general,  the  S-wave  onset  is 
much  more  gradual,  with  the  exception  of  the  Yucca  event  recorded  at  PFO,  which  looks 
very  much  like  the  earthquake  energy  curves.  It  is  not  clear  whether  these  differences 
are  due  to  near-source  or  propagation  effects.  In  the  case  of  Yucca  Flat  for  which  the 
energy  curve  at  each  station  deviates  significantly  from  the  “cleaner”  Pahute  curve,  it 
seems  likely  that  shallow  structure  in  the  source  region  is  at  least  partly  responsible  for 
the  large  amount  of  scattered  energy  in  the  waveforms. 

Several  investigations  have  discussed  the  scattering  of  locally  trapped  Rayleigh 
waves  encountering  NTS  type  structures,  see  for  example,  McLaughlin  and  Jih  (1987), 
and  Stead  and  Helmberger  (1988).  However,  explaining  all  three  components  of 
regional  NTS  events  (figure  3)  proves  especially  difficult  because  of  the  amount  of 
tangential  energy  generally  observed.  Thus,  the  scattering  must  be  due  to  a  3-D  feature 
and/or  requires  secondary  sources  such  as  spall,  etc.  In  short,  it  is  difficult  to  predict 
regional  records  using  the  conventional  RDP  formalism. 
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An  easier  problem  that  has  been  studied  quite  successfully  involves  modeling 
events  occurring  along  a  corridor  from  Imperial  Valley  to  Pasadena,  essentially  events  on 
the  San  Jacinto  fault  zone,  see  Helmberger  et  al.,  (1992),  and  Ho  and  Helmberger  (1988). 
Events  occurring  in  the  Imperial  Valley  arrive  at  Pasadena  with  extensive  coda  compared 
to  events  occurring  outside  the  basin,  see  figure  16.  The  upper  panel  displays  a  2-D 
model  connecting  Imperial  Valley  to  Pasadena  (260  km).  The  (AN)  event  occurred  near 
the  edge  while  the  (IV)  event  is  located  well  into  the  basin.  Basin  events  not  only  have 
well  developed  dispersion  but  many  times  have  secondary  arrivals.  The  shallower  the 
event,  the  stronger  these  later  arrivals.  These  signals  can  be  modeled  as  shallow  surface 
waves  propagating  slowly  in  the  upper  layer  and  re-radiating  at  the  basin  edge,  and  can 
be  the  strongest  signals  on  the  record  if  the  source  extends  to  the  surface,  see  Ho  and 
Helmberger  (1988).  Thus,  it  is  relatively  easy  to  explain  the  excess  energy  associated 
with  shallow  events  occurring  in  soft  materials. 

Applying  this  experience  to  the  NTS  data,  we  suggest  that  the  Yucca  Flat  basin 
is  responsible  for  the  complicated  records  occurring  in  figure  13  and  that  natural 
earthquakes  (figure  14),  which  normally  occur  at  deeper  depths,  are  generally  easier  to 
model.  Thus,  the  similarity  in  waveshapes  and  energy  curves  (figure  15)  allows  us  to 
model  the  large  earthquake  (M>5)  where  the  signals  are  above  the  noise,  identify  the 
paths  of  short-period  arrivals  and  predict  their  behavior  for  small  events.  In  short,  we 
think  it  is  easier  to  understand  the  short-period  phases  produced  by  the  Little  Skull 
Mountain  events  (figure  14)  than  those  produced  by  the  NTS  shots  (figure  13).  Thus,  we 
propose  to  identify  earthquakes  by  their  waveshapes  and  energy  distributions  and 
identify  explosions  as  being  not-like-earthquakes. 

In  conclusion,  we  demonstrated  that  regional  seismograms  from  earthquakes 
can  be  used  to  estimate  their  fault  parameters,  moment,  and  depths  applying  a  procedure 
developed  in  Zhao  and  Helmberger  (1994).  Next  we  examined  the  energy  content  of  the 
various  phases,  defined  Mb  (broadband  moment)  and  Mg  (energy  strength),  and 
introduced  a  new  method  of  discrimination.  In  this  method  all  events  are  processed  as 
earthquakes,  and  explosions  are  distinguished  by  their  stronger  energy  levels  relative  to 
their  long-period  amplitudes.  This  was  followed  by  a  discussion  of  a  discriminant  based 
on  the  ratio  of  short-period  (1  Hz),  vertical  component,  Pni  wavetrain  energy  to 
intermediate-period  (0.05  to  0.16  Hz),  three  component,  surface  wave  energy,  for  which 
explosions  tend  to  have  a  higher  ratio  than  earthquakes.  This  discriminant  works  on  the 
same  premise  as  the  teleseismic  Ms;fnb  ratio,  for  which  earthquakes  are  richer  in  long- 
period  surface  wave  energy  relative  to  explosions. 
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SUMMARY 

We  re-examine  the  use  of  surface-wave  magnitudes  to  determine  the  yield  of 
underground  nuclear  explosions  and  the  associated  magnitude-yield  scaling  relation¬ 
ship.  We  have  calculated  surface-wave  magnitudes  for  190  Nevada  Test  Site  (NTS) 
shots  using  regional  long-period  seismograms  from  a  combined  super-network  of  55 
North  American  stations.  Great  effort  went  towards  making  the  data  set  com¬ 
prehensive  and  diverse  in  terms  of  yield,  source  location  and  shot  medium  in  order 
to  determine  the  portability  of  surface-wave  magnitude  scales.  In  particular,  we 
examine  Pahute  Mesa,  Rainier  Mesa  and  Yucca  Flat  explosions  detonated  above 
and  below  the  water  table,  and  which  range  over  three  orders  of  magnitude  in  yield. 
By  observation  we  find  a  low-yield  measure  threshold  of  approximately  one  kiloton 
(kt)  for  (assumedly)  moderately  well-coupled  explosions  recorded  at  near-regional 
(<500  km)  stations,  which  have  little  microseismic  noise.  In  order  to  utilize  regional 
surface  waves  (A  <  15°)  for  quantifying  sources  and  for  discrimination  purposes,  we 
have  developed  related  methods  for  determining  time-domain  surface-wave  mag¬ 
nitudes  and  scalar  moments  from  regional  Rayleigh  waves.  Employing  regional 
surface-wave  data  lowers  the  effective  magnitude  threshold.  One  technique  employs 
synthetic  seismograms  to  establish  a  relationship  between  the  amplitude  of  the 
regional  Airy  phase,  or  Rayleigh  pulse  of  the  Rayleigh  wavetrain  and  an  associated 
surface-wave  magnitude,  based  on  conventional  Ms  determinations,  calculated  from 
synthetic  seismograms  propagated  to  40°.  The  other  method  uses  synthetic  seismo¬ 
grams  in  a  similar  fashion,  but  the  relationship  used  is  a  more  straightforward  one 
between  scalar  moment  and  peak  Rayleigh  wave  amplitude.  Path  corrections  are 
readily  implemented  to  both  methods.  The  inclusion  of  path  corrections  decreases 
the  Ms  variance  by  a  factor  of  two  and  affects  the  absolute  scaling  relationship  by  up 
to  a  factor  of  0.1  magnitude  units.  This  latter  effect  is  attributed  to  the  particular 
station  network  used  and  the  Green’s  functions  used  to  obtain  the  40°  Ms  values. 
Using  a  generic  structure  for  the  distance  travelled  past  the  actual  source-receiver 
path  minimizes  the  difference  between  magnitudes  determined  with  and  without 
path  corrections.  The  method  gives  stable  Ms  values  that  correlate  well  with  other 
magnitude  scale  values  over  a  range  of  three  orders  of  magnitude  in  source  yield. 
Our  Ms  values  scale  very  similarly  to  more  standard  teleseismic  Ms  values  from 
other  studies,  although  the  absolute  Ms  values  vary  by  ±0.5  magnitude  units  about 
ours.  Such  differences  are  due  in  part  to  the  choice  of  Ms  formula  used.  For 
purposes  of  future  user  comparisons,  we  give  conversion  values  to  the  previous 
studies.  Our  most  refined  Ms  values  give  the  relationship  Ms  —  1.00  X  log10  (yield)  + 
B ,  where  B  is  dependent  upon  source  region  and  shot  medium.  This  yield  exponent 
of  unity  holds  for  events  of  all  sizes  and  is  in  line  with  M^-yield  scaling  relations 
found  by  other  studies.  When  events  are  grouped  with  respect  to  source  region, 
significantly  better  fits  to  these  individual-site  linear-regression  curves  are  obtained 
compared  to  the  fits  obtained  using  a  single,  all-inclusive  model.  This  observation 
implies  that  shot-site  parameters  and  source  structure  can  significantly  affect 
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surface-wave-magnitude  measurements.  We  present  these  Ms  values  primarily  to 
augment  the  extensive  historical  analysis  of  explosion  data  based  on  surface-wave 
magnitudes  by  using  regional  data  to  increase  the  number  of  events  with 
surface-wave  magnitudes.  These  magnitudes  are  consistant  with  the  teleseismically 
determined  magnitudes  of  larger  events.  We  present  our  preferred  surface-wave 
moment  values  in  a  sequel  paper. 

Key  words:  North  America,  nuclear  explosions,  surface  waves. 


INTRODUCTION 

We  re-examine  the  use  of  surface  waves  for  determining 
underground  nuclear-explosion  magnitudes,  particularly  for 
smaller  yield,  Y  (Y<20kt),  events.  The  surface-wave 
magnitude-yield  scaling  relationship  for  such  low-yield 
events  is  not  well  defined.  Even  for  larger  yield  explosions 
there  is  some  debate  as  to  the  scaling  relation  between  yield 
and  the  long-period  energy  radiation,  as  well  as  the 
relationship  between  Ms  and  mb.  Evernden  &  Filson  (1971) 
found  that  Ms 1.4  +  1.3  X  log  (Y)  for  hard  rock  sites  in 
North  America,  where  log  is  understood  to  be  log10. 
Marshall,  Douglas  &  Hudson  (1971)  found  that  Ms  scales 
with  yield  to  the  first  power,  with  consolidated  rock  (tuff, 
salt,  granite,  andesite  and  sandstone)  coupling  10  times 
more  efficiently  than  detonations  in  alluvium.  More  recently 
Marshall,  Springer  &  Rodean  (1979)  found  that  for  events 
detonated  in  hard  rock  (salt  or  granite)  or  water-saturated 
material  (below  the  water  table)  that  Ms  =  2.16  +  0.97  X 
log  (Y)  and  that  Ms  =  1.88  +  1.06  X  log  ( Y)  for  explosions 
above  the  water  table.  Taken  together,  these  two  populations 
yield  the  relationship  Ms  =  2-05  +  log  (Y)  (Bache  1982). 
Basham  &  Horner  (1973)  found  the  scaling  relationship  for 
explosions  in  consolidated  rock  at  sites  throughout  the  world 
(a  majority  of  the  events  being  from  NTS)  to  be 
Ms  =  1*56  +  1*24  X  log  (Y).  Sykes  &  Cifuentes  (1984)  found 
a  worldwide  empirical  relationship  of  Ms  =  2.16  -1-  0.95  x 
log  (Y)  for  explosions.  Murphy  (1977)  found  that  the  scaling 
law  varied  between  events  larger  than  100  kt  [Ms  = 
1.2  + 1.33  X  log  (Y)]  and  smaller  events  [Ms  —  2.14  +  0.84  X 
log  (Y)]. 

The  above  studies  utilized  data  from  a  suite  of  sites  to 
determine  magnitude-yield  relationships.  Doing  so  is  likely 
to  add  scatter  to  the  results,  for  the  shot  medium,  the  source 
region,  and  regional  propagation  effects  may  all  affect 
surface-wave  amplitudes.  We  subgrouped  our  data  set  into 
specific  source-region  data  subsets  in  order  to  ascertain 
whether  or  not  the  separated  explosion  populations  have 
different  magnitude-scaling  relationships. 

The  data  for  this  study  are  long-period  vertical- 
component  surface -wave  records  for  190  Nevada  Test  Site 
(NTS)  events.  The  stations  used  are  from  several  North 
American  networks.  Their  respective  instruments  all  have 
passbands  that  lie  within  the  6  to  60  s  range.  Surface  waves 
are  very  useful  for  yield  estimation  purposes,  for  Ms  is 
determined  from  relatively  long-period  seismic  waves  that 
are  insensitive  to  high-frequency  near-source  effects,  which 
may  be  caused  by  asymmetries  in  the  shot  cavity  (Zhao  & 
Harkrider  1991),  spall  (Taylor  &  Randall  1989;  Day  & 
McLaughlin  1991)  or  other  possible  mechanisms.  These 


high-frequency  source  effects  may  cause  appreciable  bias  in 
magnitudes  that  are  based  on  higher  frequency  waves,  such 
as  the  mb  and  Lg  scales.  There  are  advantages  to  using 
body-wave  measurements.  Teleseismic  body  waves  traverse 
mantle  paths  that  are  relatively  homogeneous,  whereas 
surface  waves  travel  in  crustal  and  upper-mantle  waveguides 
that  are  known  to  have  strong  lateral  inhomogeneities. 

Evernden  &  Filson  (1971)  suggest,  based  on  their 
observations  of  body-wave  and  surface-wave  magnitudes  of 
U.S.  underground  explosions  detonated  both  within  and 
outside  of  NTS,  that  the  change  in  Ms  —  mb  relationship 
from  site  to  site  is  due  to  abnormal  mb  values,  rather  than 
abnormal  Ms  values,  and  that  regional-crustal  and 
upper-mantle  attenuation  variations  near  the  source  (A t*) 
are  responsible  for  the  larger  scatter  in  mb  -yield 
correlations.  Ms  measurements  are  also  less  sensitive  to 
source-depth  effects  than  are  body-wave  measured  mag¬ 
nitudes  (Marshall  &  Basham  1972).  If  it  was  not  for 
contamination  due  to  tectonic  release,  which  has  a  more 
pronounced  effect  on  long-period  surface  waves  than  body 
waves,  and  lateral  inhomogeneity  along  the  surface-wave 
propagation  path  near  the  surface  of  the  earth,  the 
long-period  energy  measured  from  surface  waves  might  be  a 
more  stable  measure  of  seismic  yield  than  teleseismic 
body-wave  measurements.  It  is  the  purpose  of  this  paper  to 
develop  and  apply  a  technique  for  reducing  the  contaminat¬ 
ing  effect  of  lateral  propagation  on  Ms  measurements. 

Another  advantage  of  using  seismic  moment  or  Ms  is  that 
empirical  evidence  and  theoretical  studies  show  that  the 
scaling  relationship  between  Ms  (or  log  moment)  and  yield 
has  an  approximate  slope  of  unity,  i.e.  Ms  —  log  (yield)  +  B, 
whereas  the  mb -yield  and  mb(Lg)~ yield  relationships  have 
slopes  between  0.65  and  0.90.  As  Evernden  &  Filson  (1971) 
point  out,  a  0.3  error  in  mb  corresponds  to  a  three-fold  error 
in  yield  determination,  while  an  equivalent  error  in  Ms 
results  in  only  a  two-fold  error  in  the  yield  estimate.  Thus 
the  error  in  yield  estimation  is  inherently  larger  when 
obtained  from  higher-frequency  magnitude  measurements. 

For  lower-yield  events  it  becomes  necessary  to  include  the 
data  from  regional  stations  (A  <  25°),  for  teleseismic 
surface-wave  recordings  have  too  low  a  signal-to-noise  ratio 
(SNR),  which  makes  them  unusable.  At  regional  distances 
surface  waves  are  not  well  dispersed,  having  a  prominent 
Airy  phase  pulse  with  a  period  between  6  and  20  s  (Alewine 
1972),  so  that  it  is  not  possible  to  measure  Ms 
conventionally  (that  is  measuring  the  amplitude  of  a  stable, 
prominent  20  s  surface  wave).  For  North  America  in 
general,  there  is  a  minimum  in  the  group  velocity  curve  near 
12  s  for  the  fundamental  Rayleigh  wave  (Marshall  et  al 
1979).  To  make  accurate  surface-wave  magnitude  measure- 
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ments,  this  energy  bandwidth  ought  to  be  modelled  as  well 
as  possible,  for  it  is  the  predominant  signal  of  the 
wavetrain. 

To  measure  Ms  we  first  model  regional  Rayleigh  waves 
with  theoretical  seismograms.  These  synthetic  wavefields  are 
then  propagated  out  to  40°,  at  which  distance  stable  20  s 
surface-wave  magnitudes  can  be  measured.  In  using  this 
procedure  several  propagation-path  models  were  tested  to 
determine  the  effect  of  attenuation  and  seismic-velocity 
structure  upon  the  Ms  values.  These  calculated  Ms  values 
remain  stable,  have  reasonably  small  errors  and  correlate 
well  with  associated  mb  magnitudes  and  log  yield  for  the 
event  data  set.  The  Ms  —  mb  relationships  are  determined  by 
a  weighted  least-square  linear  regression;  both  free-  and 
fixed-slope  curves  were  fitted  to  the  data. 

We  also  determine  time-domain  moment  measurements 
from  the  same  data.  The  moment  is  determined  from  the 
ratio  of  the  maximum  peak-to-peak  amplitude  of  the  surface 
wave  train  to  that  of  a  synthetic,  with  a  given  input-step 
moment,  propagated  to  the  same  distance  as  the  data.  These 
two  time-domain  magnitude  measurements  ( Ms  and  log  M0) 
give  very  consistent  scaling  results. 

Besides  comparing  the  Ms  results  with  several  different 
independent  magnitude  scales,  the  data  have  also  been 
separated  with  respect  to  source  region  and  shot  material. 
No  corrections  were  made  in  magnitude  for  shot-medium 
coupling  effects,  although  such  effects  can  be  considerable, 
even  for  long-period  energy  (Werth  &  Herbst  1963), 


because  such  shot-site  information  would  not  necessarily  be 
available  for  events  detonated  in  other  countries.  This  study 
is  meant  to  test  the  effectiveness  and  portability  of  a 
surface-wave  magnitude  scale  in  the  most  general  case. 

We  do  not  account  for  tectonic-release  effects  upon  the 
magnitude  measurements.  Such  effects  are  best  accounted 
for  with  moment-tensor  inversions  of  sources  which  involves 
more  sophisticated  data  analysis.  Standard  Ms  measurement 
techniques  ignore  this  factor  as  well.  The  effects  of  tectonic 
release  are  considered  in  a  sequel  paper. 

DATA 

The  data  are  long-period  vertical  seismograms  recorded  at 
North  American  stations  for  190  explosions  at  NTS  and 
consist  of  digitized  World  Wide  Standard  Seismographic 
Network  (WWSSN)  and  Canadian  Seismographic  Network 
(CSN)  records,  Digital  World  Wide  Seismographic 
Network  (DWWSN),  Lawrence  Livermore  Regional 
Seismic  Network  (LLNL)  and  Regional  Seismic  Test 
Network  (RSTN)  digital  data.  The  analogue  WWSSN  and 
CSN  data  were  digitized  by  ENSCO  Inc.  Fig.  1  shows  a  map 
of  this  58  station  super-network.  Epicentral  distances  range 
from  220  km  for  NTS  to  GSC  (Goldstone,  California),  to 
4350  km  for  NTS  to  MBC  (Mould  Bay,  Northwest 
Territories).  For  the  smaller  events,  particularly  Rainer 
Mesa  explosions,  only  the  nearer  stations  (distance  < 
1000  km)  had  usable  data.  Station  coverage  varies  widely 


Figure  1.  Map  of  the  North  American  station  network  used  in  this  study.  The  ‘spoked  wheel’  is  the  Nevada  Test  Site. 
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between  events.  22  of  the  smaller  events  (or  very  early 
events)  only  had  one  usable  station  seismogram  each,  while 
some  events  had  over  30.  The  average  number  of  stations 
reporting  per  event  is  approximately  10.  For  current  and 
future  geographical  areas  of  monitoring  interest  it  is 
reasonable  to  assume  that  only  sparse  networks  will  be  able 
to  record  any  given  event,  particularly  smaller  explosions 
(below  10  kt)  or  intentionally  ‘muffled’  explosions,  so  it  is 
important  to  see  how  well  an  explosion  magnitude  can  be 
estimated  with  only  a  few  observations. 

Because  our  methods  for  determining  magnitudes  are 
done  by  means  of  time-domain  measurements,  analogue 
records  can  be  readily  used  as  well.  We  took  advantage  of 
this  fact  to  add  considerably  more  events  (72  of  the  190)  to 
the  sample  population.  These  events  were  chosen  with  a 
mind  to  filling  out  the  data  set  with  respect  to  yield,  depth  to 
water  table  and  geographic  location. 

We  chose  to  confine  our  study  to  surface  waves  travelling 
solely  along  continental  paths,  i.e.  within  North  America. 
Surface  waves  that  propagate  across  oceanic-continental 
margins  undergo  significant  modification  in  their  waveforms 
because  of  the  great  lateral  variation  in  crustal  and 
upper-mantle  structure  at  such  boundaries.  These  propaga¬ 
tion  effects  are  not  straightforward  to  model,  and  hence 
meaningful  Green’s  functions,  or  transfer  functions,  are 
difficult  to  obtain.  Without  robust  Green’s  functions  it  is 
hard  to  infer  accurate  source  information  from  the  data. 
Also,  smaller  events  are  not  likely  to  be  observed  at  the 


distant  stations,  which  often  include  oceanic  structure  along 
their  propagation  path,  and  make  these  longer  paths  even 
less  attractive  to  include  in  the  monitoring  network. 

Of  the  190  events,  48  are  from  Pahute  Mesa,  30  are  from 
Rainer  Mesa,  105  are  from  Yucca  Flat  and  seven  others  are 
from  other  sites  in  or  around  NTS,  but  outside  of  the  three 
major  test  sites.  Of  these  seven  events,  PILEDRIVER 
(detonated  at  Climax  stock)  was  the  only  one  for  which 
digital  data  were  available.  For  some  specific  stations, 
waveforms  varied  somewhat  between  events,  depending 
upon  source  location.  The  PILEDRIVER  data  from  a  given 
station  look  appreciably  different  from  those  of  any  other 
events  recorded  at  that  same  station.  This  was  true  for  every 
station  recording  PILEDRIVER  and  probably  is  caused  by 
differences  in  the  source  region  for  this  explosion. 
PILEDRIVER  was  detonated  in  a  granitic  source  region, 
north  of  the  other  sites.  The  source-to-receiver  geometries 
for  this  event  are  approximately  the  same  as  those  at  the 
other  NTS  events,  so  the  difference  in  waveforms  appears  to 
be  a  source  effect  rather  than  a  propagation  effect.  Because 
PILEDRIVER  was  the  only  Climax  Stock  event  with 
readily  available  data,  no  further  examination  of  this  site 
was  carried  out. 

Figure  2  compares  representative  NTS  vertical-component 
long-period  data  with  synthetic  Rayleigh  waves  for  each 
source-receiver  path.  More  than  one  event  was  used  since 
no  one  event  was  observed  at  every  station.  The  darker 
traces  are  the  observations  and  the  lighter  trace  below  each 
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Figure  2.  Comparison  of  vertical-component  Rayleigh-wave  waveforms.  The  data  time  series  are  the  upper,  thicker  traces;  the  lower  trace  in 
each  case  is  the  fundamental-mode  synthetic.  All  time  series  band-passed  between  60  and  6  s. 
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is  a  synthetic  seismogram  made  with  the  fundamental-mode 
Rayleigh  wave  only.  The  seismograms  were  bandpassed 
filtered  between  6  and  60  s  to  suppress  the  long-period  and 
short-period  noise  that  would  otherwise  affect  the 
peak-to-peak  measurement  of  the  Rayleigh  pulse. 


MAGNITUDE  CALCULATION  TECHNIQUE 


We  have  developed  a  method  to  measure  surface 
magnitudes  indirectly.  Because  a  large  portion  of  the  data 
for  low-yield  events  is  from  stations  recording  at  regional 
distances  (A  <25°),  it  is  not  possible  to  calculate  Ms 
conventionally,  for  the  Rayleigh  wave  is  pulse-like,  which 
precludes  measuring  a  well-dispersed  20  s  phase  (Alewine 
1972).  We  address  this  problem  with  the  use  of  synthetic 
seismograms  of  the  fundamental  Rayleigh  wave  using  an 
asymptotic  relation  for  mixed-path  surface  waves. 

For  our  mixed  path  expressions,  we  follow  Levshin  (1985) 
or  Yanovskyaya  (1989,  page  44)  and  write  the  spectral 
Rayleigh  wave  vertical  displacement  for  approximate 
propagation  in  a  slowly  varying  laterally  inhomogeneous 
media  (e.g.  Burridge  &  Veinberg  1977;  Babich,  Chikhachev 
&  Yanovskyaya  1976;  Woodhouse  1974;  Yomogida  1985)  as 
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p(z)  is  the  local  density  distribution  in  the  medium  and  we 
have  used  Saito’s  (1967)  Rayleigh  wave  eigenfunction 
notation,  y,(z).  The  eigenfunctions  are  normalized  in  such 
a  way  that  the  vertical-displacement  eigenfunction,  yi(z)  is 
equal  to  1  at  the  free  surface,  z—  0.  This  results  in  the 
horizontal  displacement  eigenfunction,  y3(z),  being  equal  to 
the  Rayleigh  mode  surface  ellipticity  at  this  boundary.  U 
and  c  are  respectively  the  local  group  and  phase  velocities. 
By  local  we  mean  the  eigenvalues  and  eigenfunctions  that 
one  would  obtain  for  a  laterally  homogeneous  half-space 
consisting  of  the  vertical  elastic  and  density  distribution  at 
that  location.  P  is  the  receiver  location  and  PQ  is  the  point- 
source  location  and  quantities  within  the  P  or  P0  subscripted 
square  brackets  are  evaluated  at  these  locations.  The 
integrals  are  taken  along  the  phase-velocity-determined  ray 
path  between  the  two  surface  locations.  J  describes  the 
geometrical  spreading  of  the  surface-wave  energy,  y  is  the 
frequency-dependent  attenuation  coefficient  due  to  the 
anelastic  structure  of  the  path,  i.e.  y  =  a>l(2 QU),  where  Q 
is  the  attenuation  quality  factor.  The  above  expression  is 
applicable  in  the  absence  of  foci  or  shadow  zones  in  the 
vicinity  of  the  receiver.  If  there  are  foci  along  the  path  an 
additional  phase  factor  of  exp  {in  12)  should  be  included  for 
each  foci.  For  an  explosion,  W  is 
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where  M(co)  is  the  isotropic  or  explosion  spectral  seismic 
moment.  We  also  assume  a  step  for  our  explosion  history, 
i.e.  M(o))  ~  M0/(iw). 

Since  we  will  assume  that  the  directions  of  the  horizontal 
gradients  of  the  material  properties  are  approximately 
aligned  in  the  direction  of  the  source  to  receiver,  the  ray 
path  is  a  straight  line  and  7  =  r,  which  is  the  distance 
between  the  two  locations.  We  further  assume  that  the 
lateral  inhomogeneity  can  be  considered  to  be  made  up  of  n 
homogeneous  segments  of  radius  rif  i.e.  2r;-  =  r.  For 
comparison  with  Stevens  (1986),  who  used  a  similar 
expression  to  estimate  seismic  moments  for  explosions,  and 
earlier  works  on  which  his  expressions  were  based  (e.g. 
Bache,  Rodi  &  Harkrider  1978;  Harkrider  1981),  we  write 
W  in  terms  of  K  where 
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and  y2  is  the  normalized  vertical  normal  stress  eigenfunction. 
The  relation  between  K  and  W  is  obtained  by  substituting 
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into  the  previous  W  expression. 

Now  we  can  write  the  multipath  displacement  as 
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where  the  summation  convention  of  repeated  subscripts  is 
used.  The  subscript  T  denotes  the  local  quantities  for  the 
source  medium  and  the  subscript  V  the  local  quantities  at 
the  receiver.  The  shear  velocity  is  denoted  by  /3  and  the 
compressional  velocity  by  a.  For  a  given  moment,  M0,  the 
ratio  of  the  square  of  these  two  quantities  plays  a  key  role  in 
determining  the  amplitude  effect  of  various  shot  media.  To 
this  order  of  approximation  the  spectral  amplitude 
neglecting  attenuation  is  only  dependent  on  the  local 
properties  at  the  source  and  receiver.  The  attenuation  and 
phase  are  dependent  on  the  local  properties  along  the  whole 
path. 

With  the  substitution 
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(Harkrider  &  Anderson  1966;  Harkrider  1981)  and 
multiplying  by  — Vc„/Ci,  we  obtain  the  same  expression  as 
used  by  Stevens  (1986)  to  obtain  his  path  corrections  from 
NTS  to  24  WWSSN  station  in  United  States  and  Canada  and 
to  12  SRO  stations.  For  his  models  n  =  2.  The  negative  sign 
results  from  the  differences  in  our  sign  criteria  for  vertical 
displacement.  In  Stevens  (1986)  vertical  displacement  is 
positive  up  while  in  this  article  it  is  positive  down.  The  phase- 
velocity  factor  is  due  to  the  use  of  wavenumber  spreading  by 
Bache  et  al  (1978),  Harkrider  (1981),  and  Stevens  (1986) 
compared  to  geometric  spreading  by  the  others.  Bache  et  al 
(1978)  based  their  expressions  on  the  conservation  of 
lateral-energy  flux  while  these  expressions  are  from  the  main 
term  in  an  asymptotic  expansion. 
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Glover  &  Harkrider  (1986)  performed  numerical  tests  in 
order  to  estimate  the  frequency  range  for  which  these 
approximations  were  valid  for  Rayleigh  waves  generated  at 
NTS  where  the  source  region  may  be  limited  by  sharp 
boundaries  such  as  in  the  low-velocity  basin  at  Yucca  Flat. 
Rayleigh  wave  seismograms  were  calculated  for  explosive 
sources  at  depth  in  a  finite  vertical  cylinder  with  contrasting 
elastic  properties  representative  of  the  various  test  areas  at 
NTS  embedded  in  a  vertically  stratified  propagation  media. 
The  technique  couples  laterally  inhomogeneous  finite- 
element  calculations  of  the  source  region  with  Green’s 
functions  for  teleseismic  Rayleigh  waves  using  the 
elastodynamic  representation  theorem.  The  details  of  the 
technique  can  be  found  in  Harkrider  (1981)  and  Bache,  Day 
&  Swanger  (1982).  The  spectra  for  these  Rayleigh  waves 
were  then  compared  with  those,  which  used  the  two 
approximations  to  cross  the  sharp  boundary.  It  is  surprising 
that  both  approximations  worked  as  well  as  they  did  since 
they  are  based  on  a  gradual  transition.  It  was  found  that 
both  approximations  worked  equally  well  for  periods  greater 
than  four  seconds,  and  that  for  shorter  periods  the 
asymptotic  approximation  used  in  this  paper  is  better.  The 
period  range  is  dependent  on  the  material  contrast  and  the 
vertical  extent  of  the  contrast  but  this  mixed  path 
approximation  is  certainly  adequate  for  the  determination  of 
long-period  moments  and  surface-wave  magnitudes  from 
NTS  Rayleigh  wave  observations  at  continental  stations. 

It  is  interesting  to  note  that  for  this  geometry,  i.e.  n  =  2, 
the  Rayleigh  wave  transmission  coefficient,  r(a>),  of  Bache 
et  al  (1978) 


is  identical  to  the  factor  R  of  Levshin  (1985)  and  was  used  in 
both  articles  to  illustrate  the  effect  of  mixed  paths  on  the 
amplitude  of  Rayleigh  waves. 

For  each  source-to-receiver  path,  a  theoretical  Rayleigh 
wave  is  generated.  The  earth  model  used  to  create  this 
synthetic  is  meant  to  reflect  the  average  earth  structure 
between  NTS  and  the  given  station.  The  earth  models  used 
in  this  study  were  determined  from  inversions  of  dispersion 
and  attenuation  data  as  well  as  forward  modelling  of  the 
waveform  to  fine  tune  the  models.  The  criteria  for 
determining  the  goodness  of  fit  of  the  synthetic  to  the  data 
are  dispersion,  absolute  traveltime  and  waveform  fit 
(relative  amplitude  of  different  dispersed  phases).  Hence  the 
synthetic  seismogram  displays  the  same  spectral  and 
time-domain  waveform  characteristic  as  the  data  which  it 
simulates  (see  Fig.  2).  This  was  done  for  all  paths.  The 
paths  to  WWSSN  and  Canadian  stations  were  taken  from 
the  explosion  moment  study  by  Stevens  (1986).  We 
determined  the  RSTN,  LLNL  and  DWWSN  path  structures. 

To  determine  Ms  for  a  particular  source-receiver 
geometry  two  synthetics  are  generated.  One,  which  is 
propagated  the  actual  path  distance,  that  is  meant  to 
simulate  the  data  and  one  which  is  propagated  to  40°.  At  40° 
the  surface  wave  train  is  well  dispersed  and  stable,  so  that  a 
conventional  Ms  value  can  be  calculated.  Fig.  3  illustrates 
tips  method.  The  upper  set  of  seismograms  are  a  comparison 
of  data  to  its  corresponding  synthetic  seismogram.  For  this 
particular  example  the  station  COR  (Corvallis,  Oregon)  and 
the  event  LOWBALL  are  used.  The  data  is  the  solid  line 


and  the  dashed  line  is  the  synthetic  time  series.  Note  that 
the  waveform  fit  (dispersion  and  relative  amplitude)  is 
exceptional.  This  feature  is  important  in  order  to  make 
maximum  peak-to-peak  amplitude  comparisons.  The  middle 
figure  schematically  shows  the  propagation  paths  for  the 
synthetic  seismograms.  The  path  of  length  R  is  the  actual 
source-to-receiver  distance.  The  longer  path  is  of  length  40°. 
The  bottom  figures  are  of  the  two  synthetic  seismograms. 
The  left  one  is  calculated  for  the  distance  R  (10.4°  in  this 
case)  and  the  right  seismogram  is  the  one  propagated  out  to 
40°.  They  are  plotted  to  the  same  time-scale.  Note  the  much 
better  dispersed  wave  train  in  the  40°  case.  The  arrows  in 
the  right-hand  figure  mark  the  cycle  or  phase  of  the  record 
which  is  used  to  obtain  a  Ms  value. 

To  calculate  Ms  we  use  a  modified  version  of  the  von 
Seggern  formula  (von  Seggern  1977): 

Ms  =  log  (AIT)  +  1.08  X  log  (A)  +  4.38,  (9) 

where  A  is  one-half  the  maximum  peak-to-peak  amplitude 
(in  microns)  for  periods  between  17  and  23  s  of  a 
well-dispersed  wave  train  measured  from  the  vertical  record, 
T  is  the  period  of  the  arrival  measured  in  seconds,  and  A  is 
the  propagation  distance  in  degrees.  The  original  formula 
was  modified  to  include  the  period  correction  and  is  the 
same  at  a  period  of  20  s.  This  formula  was  chosen  because 
the  distance  coefficient  (1.08)  more  closely  approximates  the 
effect  of  attenuation  along  continental  paths  (Basham  1971; 
Marshall  &  Basham  1972).  Evernden  (1971)  found  the 
distance  coefficient  to  be  0.92  for  Ms  measurements  at  less 
than  25°  and  1,66  for  measurements  at  greater  distances. 
This  latter  attenuation  coefficient  is  more  characteistic  of 
mixed  continental-oceanic  path  20  s  surface  waves.  Basham 
(1971)  also  found  that  the  surface-wave  magnitude  distance 
coefficient  for  regional  Rayleigh  waves  (10  <  T  <  14  s)  is 
between  0.7  and  0.8.  Marshall  &  Basham  (1972)  make 
similar  assertions,  but  employ  a  distance  correction  which  is 
a  function  of  distance  as  well. 

A  vertical-component  measurement  has  two  advantages 
over  horizontal-component  measurements.  The  horizontal 
components  usually  have  lower  signal-to-noise  ratios  than 
the  vertical  component  and  generally  are  more  likely  to  be 
contaminated  by  Love  wave  signals  which  may  be 
generated  by  tectonic  release,  source  effects,  or  scattering 
due  to  lateral  variations  in  the  earth’s  structure. 

Both  the  regional  and  teleseismic  (40°)  synthetics  are 
generated  with  the  same  site  and  source  function,  so  that  the 
peak-to-peak  amplitude  of  the  Rayleigh  pulse  of  the 
regional  synthetic  can  be  directly  related  to  the  Ms  value 
determined  for  a  theoretical  Rayleigh  wave  train  propagated 
out  to  40°.  In  order  to  reduce  the  effect  of  the  various 
network  instruments  on  the  Ms  measurement  and  any  path 
correction,  which  might  be  made,  the  magnitude  measure¬ 
ment  was  made  on  synthetics  generated  with  the  same 
long-period  WWSSN  instrument  while  the  amplitude  of  the 
synthetic  Rayleigh  pulse  was  generated  using  the  actual 
station  instrument  response.  The  relationship  between  the 
data  peak-to-peak  amplitude  and  its  indirectly  determined 
Ms  is: 

Af5(data) 

=  Ms (synth \w)  +  log  {(PPA  \daJI(PPA  jsynth)],  (10) 
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Figure  3.  Schematic  of  the  Ms  calculation  method:  (a)  the  top  figure  is  a  comparison  of  observed-to-synthetic  seismograms  for  the  event 
Lowball  recorded  at  the  WWSN  station  COR  (A  =  10.4°).  This  record  shows  a  prominent  Airy  phase  with  a  dominant  period  that  is 
considerably  less  than  20  s.  The  solid  line  is  the  observed  time  series,  and  the  dashed  line  is  that  of  the  synthetic  seismogram.  Both  time  series 
have  been  bandpass  filtered  between  60  and  6  s.  (b)  The  middle  figure  shows  the  paths  for  which  synthetic  Rayleigh  waves  are  calculated. 
There  are  two  receiver  distances.  One,  R,  is  the  distance  between  the  actual  receiver  and  the  source.  The  other  distance  is  40°.  A  synthetic 
generated  for  the  distance  R  is  made  with  a  structure  that  best  models  the  regional  seismogram,  (c)  The  bottom  two  figures  show  synthetic 
seismograms  calculated  for  the  two  receiver  distances  for  the  COR  path  model.  The  left-hand  one  is  for  the  actual  regional  path  distance  R;  it 
is  flipped  in  polarity  with  respect  to  the  same  synthetic  in  the  top  figure.  The  right-hand  one  is  a  Rayleigh  wave  propagated  to  40°;  the  arrows 
denote  the  pulse  that  is  used  to  calculate  Ms.  Notice  that  the  dominant  period  for  this  case  is  17.5  s.  This  pulse  is  considerably  closer  in  period 
to  20  s  than  that  of  the  regional  seismogram  which  has  a  period  near  12  s. 


where  PPA  is  the  peak-to-peak  amplitude  of  the  Rayleigh 
pulse.  A  path  correction  may  be  included  in  this  expression. 
This  path  correction  is  the  difference  between  the  individual 
path  synthetic-derived  Ms  and  the  average  theoretical  Ms 
for  the  entire  network.  It  also  differs  from  a  classical  station 
correction  that  comes  from  the  data  and  not  the  synthetics. 
Secondary  station  corrections  based  on  the  data  were  not 
used  in  this  paper.  For  each  source-receiver  pair,  a  Ms  is 
calculated  from  a  synthetic  seismogram  propagated  to  40°. 
Each  such  synthetic  has  the  same  site  and  source  size,  so 
ideally  one  would  want  all  Ms  values  so  measured  to  be 
equal  in  value.  Yet  this  is  not  so,  for  each  path’s  dispersion 
and  effective  attention  at  the  periods  measured  may  be 
different.  The  difference  between  the  mean  network 


synthetic  Ms  and  a  particular  receiver  Ms  is  the  path 
correction.  A  negative  path-correction  value  implies  that 
the  theoretical  40°  station  Ms  is  larger  than  the  network 
average.  Table  1  lists  the  network  path  corrections  used. 

The  question  arises  whether  or  not  it  is  valid  to  use  the 
average  earth  structure  for  a  particular  path  to  propagate  a 
surface  wave  to  40°  when  the  earth  model  is  only  meant  to 
reflect  the  seismic  properties  of  the  earth  for  a  path  that  may 
only  be  a  small  fraction  of  this  distance.  This  is  particularly 
true  of  the  shortest  paths  for  which  the  seismic  waves 
traverse  only  western  North  America,  an  area  of  relatively 
high  attenuation  compared  to  the  continental  craton  and 
shield  areas.  A  surface  wave  propagated  40°  along  a 
characteristic  tectonic  North  American  crust  and  mantle 
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Table  1.  Network  path  corrections. 


Path 
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Correction 

Station 

Period 

Dispersion 

Single  Path 

Mixed  Path  1 

Mixed  Path  2 
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model  (NTS  to  DUG,  for  example)  for  40°  will  be  much 
more  attenuated  than  a  wave  propagated  the  same  distance 
through  an  average  structure  from  NTS  to  eastern  North 
America  (NTS  to  SCP,  for  example).  Hence  the  calculated 
Ms  for  the  NTS  to  DUG  structure  would  be  smaller  than 
the  NTS  to  SCP  Ms. 

There  are  several  methods  to  correct  for  this  path- 
dependent  effect.  As  explained  above  one  may  implement 
path  corrections  that  account  for  the  theoretical  difference 
in  attenuation  between  paths.  Another  means  is  to  make  a 
mixed-path  structure  that  has  the  appropriate  path  structure 
from  the  source  to  the  actual  station  distance,  with  the  rest 
of  the  path  out  to  40°  being  a  generic  seismic  velocity  and 
attenuation  model.  For  the  cases  in  this  study  where  the 
structures  that  comprise  the  mixed  path  are  both 
continental  structures  (i.e.  not  too  dissimilar)  the  ap¬ 
proximation  is  robust  enough  for  the  synthetic-seismogram 
calculations. 

We  have  implemented  both  procedures  individually  and 
in  conjunction  to  see  what  their  effects  are.  Another  method 
would  be  to  include  empirical  station  corrections  (Yacoub 
1983;  Given  &  Mellman  1986).  Path  correction  effects  are 
discussed  in  the  results  section. 

Besides  the  Ms  determination,  we  also  calculated  a 
time-domain  moment  for  the  same  data.  This  time-domain, 
scalar  moment  is  determined  as  follows: 

M0(data)  =  M0(synth)  X  [(PPA  \dMJ/(PPA  |synth)],  (11) 

where  PPA  is  the  peak-to-peak  amplitude  of  the  Rayleigh 
pulse  or  Airy  phase.  This  method  is  simpler  than  the  Ms 
method  and  has  the  added  advantage  that  the  synthetic 
involves  only  two  structures:  the  source  region  and  the 
propagation  path  to  the  station.  Path  corrections  were  not 
incorporated  into  the  time-domain  M0  determinations  since 
the  propagation-path  synthetic  takes  the  place  of  a  path 
correction  and  we  are  not  correcting  to  a  generic  (RSSD) 
structure.  Making  a  correction  based  on  the  difference 
between  the  average  station  value  and  some  mean  for  a 
collection  of  events  is  a  form  of  the  classical  empirical 
station  correction  and  is  most  useful  when  there  are  only  a 
few  stations  reporting  since  a  zero  sum  of  the  corrections  is 
the  usual  constraint  (Given  &  Mellman  1986).  The  mean 
moment  can  then  be  converted  to  an  Ms  using  the 
moment -Ms  relation  for  the  generic  structure  (RSSD) 
propagated  to  40°,  i.e.  the  theoretical  RSSD  station 
magnitude 

MS{PPA)  =  log  M0(PPA)  -  11.38.  (12) 

Figure  4  plots  Ms  versus  M0(PPA)  for  the  entire  data  set. 
The  correlation  between  the  two  types  of  magnitude 
measurements  is  extremely  good.  The  regression  constant 
11.43  is  very  close  to  the  theoretical  value  11.38  given  above. 
Thus  the  difference  between  our  best  mixed-path  Ms 
regression  with  moment,  and  the  Ms  moment  relation  for  a 
pure  path  of  the  generic  RSSD  model  is  only  0.05 
magnitude  units.  On  first  glance  it  might  appear  that  both 
techniques  are  identical.  However,  even  if  all  the  stations 
had  the  same  moment,  the  individual  station  magnitudes 
would  be  different  due  to  different  extended  path  lengths 
and  possibly  different  structures  from  the  station  out  to  40°. 
We  try  to  reduce  this  difference  by  making  the  additional 
correction  to  a  mean  of  the  theoretical  values  for  all  stations 


Ms(PPA)  vs.  log  Mo(PPA),  NTS 


(s.d.=0.00) 

Figure  4.  Time-domain  Ms  values  regressed  versus  time-domain 
moments  for  all  the  events  in  this  study.  Note  the  extremely  good 
correlation  between  the  two  scales.  Vertical  error  bars  are  the 
variance  for  the  individual  events. 

in  the  mixed-path  evaluation.  These  observations  imply  that 
the  RSSD  model  is  a  good  average  model  for  the  network 
and  that  our  Ms  calculations  are  sound  and  result  in  robust 
measurements  of  surface-wave  magnitude,  which  are  not  too 
dependent  on  which  of  the  two  techniques  we  use.  In  our 
analysis  of  the  data  we  will  use  only  Ms  measurements.  A 
table  of  the  logA/0  and  their  standard  deviations  will  be 
given  in  Woods  &  Harkrider  (1995). 

Source  structure  significantly  affects  absolute  surface-wave 
amplitudes,  and  hence  surface-wave  magnitudes  or  mo¬ 
ments.  For  all  of  the  synthetic  seismograms  generated,  we 
used  the  Stevens  (1986)  and  Given  &  Mellman  (1986) 
NTS-source  elastic  structure,  which  is  basically  a  Pahute 
Mesa  velocity  structure.  By  numerical  simulations  using  a 
variety  of  different  NTS  structures,  we  found  that  for  the 
frequencies  of  interest  and  sources  in  the  upper  6  km  the 
primary  effect  was  due  to  the  difference  in  shot-point 

velocity  ratios.  The  size  of  the  effect  can  be  predicted 

extremely  well  from  their  explicit  presence  in  the 

mixed-path  expression,  eq.  (6).  As  an  example,  our  Green’s 
functions  are  computed  for  an  explosive  source  at  a  depth  of 
600  m.  In  the  Stevens  (1986)  source  structure,  the  second 
layer  starts  at  a  depth  of  500  m.  There  is  a  significant 
difference  between  the  Poisson’s  ratio  of  the  surface  and 
second  layer  in  the  source  earth  structure.  The  log  difference 
between  the  square  of  their  compressional  to  shear 
velocities  would  predict  from  eq.  (6)  an  Ms  difference  of 
0.17.  The  actual  difference  between  the  Ms  of  a  surface 
explosion  and  our  Green’s  function  is  0.16  with  the 

near-surface  explosion  smaller  as  predicted.  In  order  to 
reduce  the  effect  of  differing  shot-point  velocity  ratios, 
Stevens  (1986)  suggested  a  new  explosion  moment,  Mq, 
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defined  by 

B2 

M^3~M0.  (13) 

a 

For  a  shot-point  medium  with  Poisson’s  ratio  of  0.25 
(a2 /p2  —  3)  the  value  of  the  moment  is  unchanged. 

In  Fig.  2,  we  see  that  for  the  WWSSN  stations,  denoted  by 
three  letters,  the  period  of  the  dominant  phase  is 
significantly  lower  than  the  recommended  lower  cut-off  of 
17  s  for  the  standard  Ms  formula.  This  period  was 
determined  by  taking  twice  the  time  difference  between  the 
arrival  of  the  largest  peak  and  trough.  We  also  calculated 
the  ‘instantaneous  period’  of  this  arrival  and  found  it  to  be 
essentially  the  same  value.  The  dominant  period  at  each 
station  is  given  in  Table  1.  For  the  WWSSN  stations,  the 
periods  are  between  10  and  15  s.  Most  are  near  11s.  For  the 
digital  stations,  denoted  by  four  letters,  the  dominant  period 
is  between  15  and  19.5  s,  with  the  average  being  16.5  s.  An 
alternative  approach  for  using  the  maximum  amplitude  of 
Rayleigh  wave  observations  where  the  dominant  period  is 
significantly  different  from  20  s  was  developed  by  Marshall  & 
Basham  (1972).  Using  the  stationary-phase  approximation 
they  determined  a  path  correction,  which  corrected  for  the 
dispersive  characteristics  of  the  path.  Using  observed 
dispersion  curves  for  North  America,  Eurasia,  mixed 
ocean-continent,  and  pure  ocean  paths,  they  were  able  to 
determine  an  Ms  correction  based  on  the  period  of  an 
observed  Airy  phase  to  the  20  s  period  arrival  in  North 
America  or  Eurasia.  The  North  American  dispersion 
correction  appropriate  for  the  dominant  period  measured  at 
each  station  in  our  network  is  also  given  in  Table  1  (column 
2).  An  advantage  of  our  technique  is  that  our  path 
corrections  are  independent  of  recording  instrument 
whereas  Marshall  &  Basham’s  correction  depends  on  the 
dominant  period,  which  depends  not  only  on  dispersion  but 
also  instrument  response.  As  an  example,  the  station  ALQ 
has  a  dominant  period  of  12  s  and  ALQD  has  a  dominant 
period  of  19.5  s. 

An  advantage  that  time-domain  estimates  of  Ms  or  M0 
have  over  spectral  estimates  can  be  seen  in  Fig.  2.  Except 
for  the  work  of  Patton  (e.g.  Patton  1991),  the  Green’s 
functions  used  for  spectral  estimates  of  explosions  have  been 
fundamental  Rayleigh  and  Love  waves.  As  can  be  seen  from 
the  figure,  it  is  very  important  to  isolate  the  fundamental 
surface  wave  in  the  data  for  taking  its  spectra  for  moment 
estimates.  The  Rayleigh  waves  at  almost  every  station  show 
the  additional  presence  of  higher  modes.  The  higher  modes 
are  primarily  due  to  constructive  interference  of  multiple 
reflected  shear  waves  and  are  therefore  very  sensitive  to 
lateral  variations  in  crust  and  upper-mantle  structure.  This  is 
especially  true  for  non-parallel  layers  with  sharp  contrasts. 
Therefore,  in  the  presence  of  nearby  signals  or  noise,  it 
makes  more  sense  to  use  the  larger  time-domain  amplitudes 
of  the  fundamental-mode  Airy  phases  at  regional  distances. 
Because  of  the  possibility  of  tectonic  release,  it  is  also 
necessary  to  determine  the  polarity  of  the  surface  wave. 
Again  this  is  best  done  in  the  time  domain,  especially  for 
Love  waves. 

If  a  spectral  estimate  is  desired,  comparing  the  Green’s 
function  with  the  data  in  the  time  domain  should  allow  one 
to  determine  time  windows  and  tapers  so  as  reduce  the 
contamination  of  spectral-amplitude  estimates  with  higher 


modes  and  spurious  scattered  arrivals  at  intermediate 
ranges.  And  at  close  ranges  where  this  may  not  be  possible, 
it  should  help  in  deciding  which  time-domain  amplitude 
measurements  best  represent  the  spectral  amplitudes  of  the 
fundamental  modes. 

The  question  remains  how  well  do  either  of  these  two 
measurements  compare  to  spectral-moment  estimates.  For 
the  events  for  which  digital  data  were  available, 
spectral-domain  moments  were  also  determined.  Spectral 
moments  were  calculated  using  the  method  of  Stevens 
(1986),  with  the  exception  that  station  corrections  were  not 
included  in  our  moment  calculations.  Spectral  moments 
were  calculated  in  the  bandwidth  between  6  and  60  s.  These 
spectral  moments  will  be  referred  to  as  M0.  Moments  were 
also  obtained  by  inverting  for  an  isotropic  (explosion 
source)  component  (M{)  and  a  deviatoric  component 
(caused  by  tectonic  release  or  an  asymmetric  source  cavity) 
of  moment  (M#).  Details  of  these  moments  are  the  subject 
of  a  later  paper  by  the  authors  (Woods  &  Harkrider,  in 
preparation). 

We  compare  the  time-domain  moments  with  these  two 
types  of  spectral-domain  moments.  Fig.  5(a)  compares 
M0(PPA)  to  Mq(co)  and  Fig.  5(b)  compares  M0(PPA)  to 
M0((o)  refers  to  the  average  spectral  scalar  moment 
and  M7(o))  refers  to  the  isotropic  source  component 
determined  from  a  moment-tensor  inversion  scheme 
(Woods  &  Harkrider,  in  preparation).  M0(PPA)  correlates 
well  with  the  two  types  of  spectral  moments.  The  advantage 
of  time-domain  moments  is  that  analogue  data  can  be  used 
directly  and  the  effective  SRN  is  lower  than  for 
spectral-moments  methods,  thus  smaller  events  can  be 
measured. 

In  the  top  figure  there  is  some  scatter  in  the  moment 
correlation  for  several  of  the  smaller  events,  with  the 
time-domain  moments  being  significantly  larger  than  the 
spectral-domain  moments.  Most  of  these  outlying  events  are 
Rainier  shots,  none  is  from  Yucca  and  only  one,  REX,  is 
from  Pahute.  REX  (M0(PPA)  =  15.35,  M0(a>)  =  14.87)  was 
detonated  below  the  water  table  and  had  an  anomalously 
large  deviatoric  moment  component  (Woods  &  Harkrider, 
in  preparation).  Since  many  more  stations  were  used  in  the 
determination  of  the  time-domain  moments  with  a  more 
complete  azimuthal  coverage,  one  would  expect  them  to  be 
less  affected  by  the  deviatoric  component,  which  at  NTS  is 
speculated  to  be  such  that  it  will  average  out  with  adequate 
azimuthal  coverage.  The  scatter  is  somewhat  less  in  the 
M0(PPA):Mi((d)  curve  (Fig.  5b).  In  particular,  REX  no 
longer  stands  out,  although  several  of  the  Rainer  events  lie 
well  off  the  scaling  curve.  These  outlying  Rainier  events  can 
be  explained  in  several  ways.  First,  these  events  are 
relatively  small  and  are  only  measured  at  very  few  stations 
(sometimes  only  one  to  three  stations),  thus  the  scatter,  or 
error,  in  the  moment  measurement  is  larger.  Also  the 
spectral  moment  is  more  susceptible  to  noise  contamination 
since  it  requires  the  isolation  of  a  segment  of  the  surface 
wave  train  from  noise  and  the  time-domain  measurement 
only  requires  that  the  maximum  amplitude  not  be  too 
distorted  by  noise.  One  problem  with  this  explanation  is  that 
there  are  other  small  events  recorded  at  Pahute  Mesa  and 
Yucca  Flat  that  lie  right  on  the  moment  scaling  curve  (Fig. 
5a)  and  these  events  are  no  better  recorded  than  the  Rainier 
events.  Another  possibility  is  that  these  outlying  events 
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Mo(PPA)  vs  Mo(spectral) 


Log  Mo(spectral)  N-M 

Log  Mo(PPA)  N-M  =1.00  x  Log  Mo(spectral)  N-M  +  0.03 
(s.d.=0.1 1) 


Mo(PPA)  vs  Mi(spectral) 


Log  Mo(PPA)  N-M  =1.00  x  Log  Mi(spectral)  N-M  +  0.02 
(s.d  =0.10) 

Figure  5.  Time-domain  log  moments  regressed  against  spectral- 
domain  log  moments.  In  the  top  figure  spectral  moments  were 
determined  assuming  an  isotropic  source  only,  while  in  the  bottom 
figure  the  spectral  moments  were  determined  by  inverting  for  a 
isotropic  source  +  a  double-couple  source.  The  regressions  were 
constrained  to  a  slope  of  unity. 

reflect  differences  in  source  spectra.  As  discussed  previously, 
the  time-domain  moments  measure  energy  predominantly  in 
the  10-14  s  range,  the  period  range  of  the  continental  Airy 
phase,  whereas  the  spectral  moment  is  an  average  of  the 


spectral  ratio  between  6  and  60  s.  So,  it  is  possible  that  the 
Rainier  test  sites  excite  more  high-frequency  energy  than 
do  either  the  Phaute  or  Yucca  site.  This  effect  was  seen  in 
data  at  several  of  the  closer  stations,  in  particular. 

DATA  ANALYSIS  AND  RESULTS 

The  seismograms  were  bandpassed  filtered  between  6  and 
60  s  to  minimize  contaminating  noise  as  described 
previously.  The  vertical  records  were  visually  inspected  to 
ensure  that  the  correct  time  window  was  used  and  that  their 
signal-to-noise  ratio  was  about  2.0  (approximately).  M$ 
values  were  then  calculated  for  the  data  as  per  the  method 
described  above  (eq.  9)  with  several  variations.  The 
synthetic  seismograms  were  also  bandpassed  filtered 
between  6  and  60s  for  consistency.  The  Ms  values  are 
plotted  against  seismic  magnitudes  of  several  scales  for  the 
same  set  of  events.  It  should  be  noted  that  no  complete 
magnitude  list  was  available  for  all  190  events. 

We  chose  to  compare  or  plot  our  data  primarily  with 
body-wave  magnitudes  determined  by  Lilwall  &  McNeary 
(1985).  The  Lilwall-McNeary  (LM)  data  set  contains  143  of 
the  190  events  examined  by  us  and  is  believed  to  be  a 
well-determined  and  self-consistent  list  of  mb  values  that 
have  small  errors  due  to,  among  other  things,  the  inclusion 
of  network  station  corrections.  Fig.  6  shows  the  mb -yield 
relationship  for  events  in  this  study  for  which  mb  and  yield 
information  were  available.  It  is  important  to  notice  that 
events  above  and  below  the  water  table  separate  into  two 
distinct  populations.  For  this  data  set  this  separation  is  only 
apparent  near  the  cluster  of  events  with  mb  s  around  5.4. 


mb(LM)  vs.  Log  Y 


Figure  6.  Lilwall  mb  versus  log  Yield  for  events  from  this  study. 
The  solid  line  is  the  best-fitting  regression  line.  The  dashed  lines 
show  the  95  per  cent  confidence  interval  of  this  line.  Lines  through 
the  data  points  represent  one  standard  deviation  in  a  datum 
measurement.  Blackened  circles  represent  sources  beneath  the 
water  table,  open  circles  are  events  above  the  water  table  and  open 
squares  are  events  for  which  this  information  is  not  known.  Events 
detonated  below  the  water  table  have  a  larger  mb  for  a  given  yield. 
Besides  this  separation  of  data,  there  is  little  scatter  to  the  data.  The 
consistency  of  the  mb -yield  relationship  makes  it  reasonable  to  use 
these  mb  values  to  plot  our  Ms  measurements  against. 
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Also  notice  the  very  small  error  bars  for  this  data;  for  many 
events  the  error  bars  are  smaller  than  the  symbols 
demarking  a  data  point.  The  solid  line  is  the  best-fitting, 
least-squares  curve,  with  the  dashed  curves  being  the 
two-sigma  confidence  interval  of  the  regression  relationship. 
The  correlation  between  mb  and  yield  is  good,  with  the 
scatter  mostly  being  due  to  the  above-water  table  shots.  The 
slope  of  the  regression  curve  is  0.67. 

This  scaling  curve  slope  is  slightly  lower  than  that  found 
in  other  studies  of  teleseismic  mb  -yield  scaling  relation¬ 
ships.  Marshall  et  al  (1979)  found  that  mb  was  proportional 
to  Y074  for  well-coupled  Yucca  flat  explosions,  and 
proportional  to  Y078  for  explosions  throughout  NTS  and 
Amchitka.  Longer  period  telesiesmic  body-wave  magnitudes 
mLPP  introduced  by  Basham  &  Horner  (1973)  show  that  for 
events  in  tuff  and  rhyolite  the  amplitude  of  the  arrivals  is 
proportional  to  Y0'72.  Murphy  (1977)  compared  theoretical 
mh  -yield  scaling  relations  for  cube-root  scaling  models  and 
the  modified  Mueller  &  Murphy  (1971)  source  model.  He 
found  that  the  yield  exponent  varies  between  0.6  and  1.0  for 
the  cube-root  model  in  the  yield  range  of  interest,  whereas 
the  exponent  is  a  constant  0.85  for  their  modified  source 
model.  Schlittenhardt  (1988)  found  mb  to  be  proportional  to 
Y0'82  for  NTS  explosions.  The  empirically  derived  curves 
have  errors  in  their  slopes  of  the  order  of  0.05  to  0.1  units 
and  are  based  on  small  sampling  populations.  The  LM 
mb  -yield  scaling  relationship  is  determined  from  a 
significantly  larger  data  set,  making  it  at  least  as  reliable  as 
any  other  empirical  scaling  curve. 

The  same  scaling  law  slope  (-0.67)  holds  for  the  LM  data 
when  they  were  separated  with  respect  to  test-site  and 
shot-medium  coupling  (whether  detonated  above  or  below 
the  water  table).  There  is  consensus  in  the  literature  that 
that  seismic  coupling  is  a  function  of  the  percentage  dry  (or 
gaseous)  porosity  of  a  material.  In  a  study  of  small-scale, 
high-explosive  experiments  with  15  rock  types,  Larson 
(1981)  found  for  a  given  size  explosion  that  the  elastic  radius 
of  a  porous  material  (such  as  tuff)  increased  with  increasing 
water  content.  The  dominant  non-linear  mechanism  (within 
the  plastic  radius)  working  at  low  yields  appears  to  be  pore 
crushing  of  the  surrounding  material  (Stevens  1991). 
Non-linear  finite-difference  calculations  (Bache  1982)  also 
indicate  that  porosity  is  the  most  important  characteristic  of 
NTS  tuff  for  seismic-coupling  purposes.  In  the  same  study, 
source  functions  for  Yucca  Rat  wet  and  dry  tuff  are 
significantly  different,  with  the  long-period  amplitude  of  wet 
tuff  being  larger  by  50  per  cent  and  its  corner  frequency 
being  lower.  Springer  (1966)  has  observed  this  effect  for 
teleseismic  P-wave  amplitudes.  Patton  (1988),  Gupta  et  al 
(1989),  and  Vergino  &  Mensing  (1989)  have  observed  this 
coupling  effect  in  regional  phases  such  as  Lg,  Pn  and  P  . 

Several  sets  of  synthetic  Rayleigh  waves  were  generated 
at  the  40°  distance  for  calculating  Ms.  One  set  was 
propagated  along  the  single-structure  model  (hereafter 
referred  to  as  the  single-path  case)  which  reflects  the 
average  earth  structure  between  NTS  and  a  given  station. 
We  also  generated  mixed-path  synthetics  for  which  that 
part  of  the  path  beyond  the  actual  source-receiver  distance, 
out  to  40°,  the  surface  wave  is  propagated  along  a  generic 
earth  structure.  The  NTS-RSSD  Earth  structure  was  chosen 
for  this  generic  path  section,  as  it  is  a  relatively  simple 
structure  that  generates  stable  surface  waves  and  it  is 


roughly  an  intermediate  range  station  (distance  <  1900  km), 
so  that  its  structure  can  be  considered  to  be  an  ‘average’ 
North  American  structure  for  the  network. 

Surface  magnitudes  were  first  calculated  from  the  40° 
synthetics  generated  with  a  single-structure  propagation 
path.  Figs  7(a)  and  (b)  display  single-path  ms  values, 
calculated  as  described  above,  versus  body-wave  magnitude 
(j mb ).  These  mb  s  are  those  of  Lilwall  &  McNeary  (1985).  In 
the  upper  figure  (7a)  Ms  is  calculated  without  path 
corrections,  whereas  Ms  is  calculated  with  path  corrections 
in  the  lower  figure.  The  vertical  error  bars  represent  the  one 
standard-deviation  confidence  interval  for  each  Ms  value. 
The  solid  line  is  the  best-fitting  weighted  least-squares 
regression  of  the  data,  with  the  weighting  factor  being 
inversely  proportional  to  individual  event  variances.  The 
dashed  lines  represent  the  two  standard-deviation  error 
(assuming  a  Student  t  distribution)  of  the  fit  of  the  line  to 
the  data.  Since  Student  t  statistics  (Lapin  1983)  are  functions 
of  the  sample  size,  the  confidence  level  will  vary  with  data 
set  size.  Solid  black  circles  are  shots  below  the  water  table, 
open  circles  are  shots  above  the  water  table,  and  open 
squares  are  shots  for  which  this  information  is  not  known. 
Note  the  error  bars  are  approximtely  50  per  cent  larger  for 
the  uncorrected  Ms' s  (Fig.  7a)  than  for  the  case  of 
path-corrected  Mf  s  (Fig.  7b).  The  scatter  in  the  data  is  also 
significantly  less  for  the  path-corrected  Msy s  as  evidenced  by 
the  reduction  of  the  standard  deviation  of  the  free-slope 
regression.  Although,  it  should  be  noted  here  and  for  the 
following  discussion  that  the  success  of  these  techniques  is 
not  in  how  well  they  reduce  the  standard  deviation  from 
some  assumed  linear  relation  between  surface  and 
body -wave  magnitudes,  but  how  well  they  reduce  the 
variance  or  standard  deviation  of  the  magnitude  determina¬ 
tion  of  an  individual  event,  and  even  more  important,  their 
ability  to  include  small  events  in  the  data  base  where 
teleseimic  data  was  too  sparse  for  determining  a  classical 
Ms.  On  the  other  hand,  it  is  interesting  that  the  standard 
deviation  for  the  LM  body-wave  magnitudes  versus  yield  is 
0.15  and  our  best  regression  between  Ms  and  their  mb  for 
the  same  explosion  data  set  is  0.26.  Both  regressions  are 
improved  when  the  events  are  separated  into  different 
populations.  When  we  get  to  Figs  15  and  16,  which  are  our 
regressions  on  log  yield,  you  will  see  that  our  regression 
standard  deviations  are  quite  competitive  with  those  of  mb 
for  events  restricted  to  NTS  below  the  water  table  and 
Rainier. 

There  are  two  significant  effects  of  including  path 
corrections.  One  is  the  reduction  in  variance  of  individual 
magnitudes.  Without  path  corrections  the  individual  station 
magnitudes  have  a  bi-modal  distribution  reflecting  the  two 
generic  earth  models  of  North  America:  the  tectonic  western 
and  stable  cratonic  eastern  crust  and  upper-mantle 
structures.  The  path  corrections  bring  in  the  outlying 
station  magnitude  values  towards  the  mean  value.  Includ¬ 
ing  path  corrections  for  the  single -path  derived  A/5’s 
increases  the  average  value  by  0.14  units  (or  32  per  cent). 
This  effect  can  be  attributed  to  the  smaller  events  that  are 
brought  more  in  line  with  the  curve  containing  larger  events. 
This  in  turn  is  due  to  the  fact  that  the  smaller  events  are 
only  observed  at  nearer  stations  in  tectonic  North  America 
(TNA),  for  which  path  structures  exhibit  higher  attenuation 
than  do  more  cratonic  or  shield-like  models,  so  that  surface 
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Ms(npc)  vs.  mb,  single  path 


(s.d.=0.35) 


Ms(pc)  vs,  mb,  single  path 


(s.d.=0.25) 


Ms(npc)  vs.  mb,  mixed-path  1 


(s.d.=0.25) 


Ms(npc)  vs.  mb,  mixed-path  2 


(s.d.=0.28) 


Ms(pc)  vs.  mb,  mixed-path  1 


(s.d.=0.26) 


Ms(pc)  vs.  mb,  mixed-path  2 


(f)  Ms  *1.61  x  mb -4.75 
(s.d.=0.25) 


Figure  7.  Here  Ms  is  plotted  versus  Lilwall  mb' s.  For  the  left  figures  Ms  is  calculated  with  single-path  Green’s  functions,  without  path 
corrections  (npc)  and  with  path  corrections  (pc).  The  best-fitting  regression  model  is  the  solid  line  running  through  the  data  points.  The  dashed 
lines  are  the  two-sigma  confidence  intervals  of  the  line.  The  Ms-mb  relationship  and  the  rms  error  of  the  data  are  at  the  bottom  of  each  figure. 
The  middle  figure  Ms  values  are  determined  using  the  mixed-path-1  synthetics  and  in  the  right  two  figures  mixed-path-2  Green’s  functions  were 
used. 


waves  propagated  along  a  TNA  path  for  40°  will  be 
significantly  more  attenuated  than  waves  propagated  along  a 
craton  or  shield  path  for  that  same  distance.  Path 
corrections  reduce  this  effect  significantly  for  the  single-path 
derived  magnitudes.  Table  1  lists  these  network-path 
corrections.  The  third  column  lists  the  corrections  for 
single-path  synthetics.  A  positive  value  denotes  that  the  Ms 
for  a  station  is  smaller  than  the  network  theoretical 
average. 

We  next  explored  the  effect  of  mixed-path  transfer 
functions  upon  the  Ms  calculations.  As  described  above,  we 
chose  the  path  to  RSSD  as  a  generic  structure  for  the  second 
portion  of  the  mixed-path  synthetic  seismogram  calculations. 
We  generate  two  sets  of  these  synthetics.  The  difference 
between  these  two  mixed-path  earth  structures  is  their 
spectral  attenuation  coefficients,  with  y(a>)  being  twice  as 
large,  at  a  given  frequency,  for  the  mixed-path-2  case  as  for 
the  mixed-path-1  case.  Fig.  8  is  a  plot  of  the  attenuation 
factor  (y(w))  as  a  function  of  period.  The  line  labelled 
RSSD  X  2  is  that  of  the  increased  attenuation  model.  It  is 
referred  to  as  ‘mixed-path  T  throughout  this  study.  The 
lower,  dashed  curve  is  the  attenuation  curve  for  the  RSSD 
structure.  Synthetics  made  with  this  RSSD  generic  structure 
for  the  latter  portion  of  the  40°  travel  path  will  be  referred 


to  as  ‘mixed-path  1’.  Table  1  gives  the  path  corrections  for 
each  station  for  these  two  cases,  also. 

The  network  averages  and  their  standard  deviations  for 
the  three  different  path  assumptions  are 

Single  path 

MS(PPA)  =  log  M0~  11.50, 

with  a  standard  deviation  of  0.38. 

Mixed-path  1 

MS(PPA)  =  log  M0-  11.44, 

with  a  standard  deviation  of  0.17. 

Mixed  path  2 

MS(PPA)  =  log  M0  —  11.64, 
with  a  standard  deviation  of  0.22. 

These  standard  deviations  are  only  a  measure  of  the 
spread  of  the  theoretical  network  path  corrections  for  a 
given  moment.  The  above  standard  deviations  could  also  be 
found  from  squaring  and  summing  the  path-correction 
values  in  their  respective  columns  in  Table  1  after  removing 
the  three  digital  stations  that  have  the  same  correction  as 


(14) 

(15) 

(16) 
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Explanation 


RSSDx2 


- RSSD 


Figure  8.  The  two  attenuation  models  for  the  generic  portion  of  the  hybrid  propagation  model  are  plotted  versus  period.  Gamma  is  the 
attenuation  coefficient  (eq.  1)  at  a  given  period.  Model  RSSD  X  2’s  attenuation  is  twice  that  of  the  RSSD  model.  Ms  values  calculated  with 
synthetics  using  model  RSSD  are  referred  to  as  mixed-path  1,  while  values  determined  from  synthetics  created  using  attenuation  model 
RSSD  X  2  are  referred  to  as  mixed-path  2. 


their  analogue  entries.  It  should  be  remembered  that  the  40° 
Ms  were  all  calculated  using  the  same  WWSSN  long-period 
instrument  response  and  thus  the  path  corrections  are 
instrument-independent.  In  relating  theoretical  station 
amplitudes  to  Ms,  the  actual  station  instrument  response  is 
removed.  The  previously  discussed  regression  constant  of 
11.43  for  the  entire  data  set  differs  slightly  from  the 
theoretical  network  value  of  11.44  given  in  eq.  (15)  primarily 
because  the  individual  observed  values  are  not  always 
determined  using  the  full  network  of  stations. 

In  Figs  7(c)  and  (d)  the  Ms  magnitudes  were  calculated 
using  synthetic  seismograms  using  the  mixed-path-1  model. 
In  Fig.  7(c)  the  Ms* s  are  calculated  without  path-correction 
terms,  while  in  Fig.  7(d)  path  corrections  are  included.  The 
addition  of  the  path-correction  terms  cuts  the  data  variance, 
but  by  no  more  than  25  per  cent,  and  then  not  in  all  cases. 
Assuming  a  fixed-slope  regression  (m  =  1.50),  there  is  no 
offset  in  the  intercept  between  the  uncorrected  and 
path-corrected  Mf  s.  So  using  a  generic  structure  for  the 
remainder  of  the  40°  path  acts  as  a  path  correction  as  well. 

Figures  7(e)  and  (f)  are  Ms  versus  mb  plots  for  the 
mixed-path-2  case  without  and  with  path  corrections, 
respectively.  It  should  be  noted  that  the  slope  of  the 
regression  line  is  nearly  the  same  (1.61)  for  all  three 
path-corrected  cases  (Figs  7b,  d  and  f)  and  that  the 
difference  in  the  regression  intercept  between  the  three 
cases  is  essentially  what  one  would  predict  from  the  three 
theoretical  network  intercepts  given  in  eqs  (14),  (15)  and 
(16).  This  should  not  be  surprising  since  the  determined 
moment  for  each  event  is  the  same  from  figure  to  figure  and 
on  average  their  slopes  should  be  equal  and  the  differences 
in  Ms  between  the  three  figures  should  be  close  to  the 


difference  in  the  theoretical  network  values.  If  explosions  at 
only  one  site  are  plotted,  the  slope  of  the  curve  is  closer  to 
1.5,  so  that  we  will  take  the  Ms-mb  scaling  relationship  to 
be: 

Ms  =  1.50  Xmb+  B.  (17) 

For  this  fixed-slope  scaling  relationship,  the  uncorrected  and 
path-corrected  mixed-path-1  Ms  curves  have  the  same 
intercept,  whereas  for  the  mixed-path-2  case  the  intercept  is 
0.10  units  larger  for  the  path-corrected  curve  than  for  the 
uncorrected  curve.  For  path-corrected  Ms' s,  the  mixed-path- 
1  intercept  is  0.21  log  units  greater  than  that  of  the 
mixed-path-2  curve.  Because  mixed-path  synthetics  are 
propagated  to  40°  along  a  more  attenuative  path,  the  Ms 
measured  also  will  be  smaller. 

For  the  single-path  case,  path-corrected  Ms  values  give 
the  same  relationship  (slope  =  1.61),  but  the  slope  is  larger 
(1.66)  for  the  uncorrected  magnitudes,  although  the 
difference  lies  within  the  errors  bounds.  It  would  seem  that 
both  path  corrections  and  mixed-path  Green’s  functions 
improve  Ms  determinations  for  the  method  used  here.  The 
most  consistent,  reliable  magnitudes  are  obtained  using 
mixed-path-generated  synthetics  in  conjunction  with  path 
corrections  for  the  40°  Ms  measurements.  The  variance 
among  the  mixed-path-based  Ms  values  for  the  network  is 
smaller  than  that  when  Ms  is  derived  from  single-path 
synthetics,  so  that  magnitude  measurements  will  be  more 
consistent  when  they  are  determined  from  mixed-path 
synthetics.  This  is  particularly  important  for  events  with  few 
reporting  stations.  All  further  plots  of  Ms  in  this  study  use 
values  obtained  from  the  mixed-path-1  case  with  path 
corrections,  unless  stated  otherwise. 
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Figure  9.  Ms  of  this  study  regressed  against  those  determined  by  other  studies. 


How  well  the  final  Ms  values  reflect  the  actual  seismic 
magnitude  of  these  events  necessitates  having  another 
measure  of  their  size.  In  the  event  of  anomalously  high  or 
low  seismic-source  coupling,  for  example,  both  body  waves 
and  surface  waves  should  be  affected  similarly  by  coupling 
effects.  A  magnitude  parameter  independent  of  seismic 
observations  would  be  useful  to  plot  the  Ms  against,  so  we 
have  also  fitted  our  results  to  estimate  log  yields.  These 
relationships  are  shown  and  discussed  later  in  the  paper. 
Yield  values  are  estimated  to  be  within  10  per  cent  of  the 
actual  yield  (Springer  &  Kinnaman  1971).  Yield  information 
was  available  for  174  of  the  events,  thus  yields  make  up  the 
most  comprehensive  data  set  with  which  to  compare  our 
results.  The  yields  for  this  data  set  range  over  three  orders 
of  magnitude  in  size.  The  greatest  scatter,  as  in  the  case  of 
mb  versus  log  yield,  is  due  to  shots  above  the  water  table.  It 
should  also  be  kept  in  mind  that  the  scatter  would  be 
further  reduced  if  the  data  were  separated  into  populations 
based  on  their  location  at  NTS  (i.e.  Pahute  Mesa,  Rainier 
Mesa  and  Yucca  Flat). 

Since  our  magnitude  values  are  based  on  theoretical 
continental  structures,  as  well  as  the  particular  network 
used,  we  wanted  to  compare  our  Ms  values  to  those 
obtained  from  more  standard  Ms  methods.  In  addition,  the 
reference  distance  of  40°  is  arbitrary  and  along  with 
differing  magnitude  formulae  will  result  in  an  offset  from 
previous  studies.  In  order  to  make  it  possible  for  readers  to 


convert  our  values  to  the  results  of  others,  we  regressed  on 
magnitudes  of  mutual  events.  Fig.  9  shows  our  Ms  values  (x 
axis)  versus  those  from  six  other  studies  (y  axis)  (Basham 
1969;  Marshall  &  Basham  1972;  Marshall  et  al.  1979; 
Basham  &  Horner  1973;  Yacoub  1983;  von  Seggern  1973). 
The  overlap  in  data  sets  varies  between  8  and  16  events. 
We  performed  a  fixed-slope  (slope  =  1.0),  linear  regression 
of  our  Ms  values  to  those  of  the  six  outside  studies;  in 
general  the  correlation  is  very  good.  It  is  important  to  note 
that  with  our  method  we  are  able  to  measure  Ms  for  events 
one  tenth  the  size  of  the  smallest  events  measured  in  the 
other  studies  (i.e.  Ms  =  1.75  to  2.0).  This  is  after  having 
corrected  for  differences  in  absolute  Ms  scales.  We  are  able 
to  measure  Ms  for  these  smaller  events  because  we  are  able 
to  make  use  of  near-regional  (<500  km)  records  with  the 
method  described  in  this  paper. 

The  offsets  in  Ms  values  vary  considerably.  This  offset  is 
due  in  part  to  the  difference  in  definition  of  Ms  for  each 
study,  in  particular  the  distance-correction  term.  As 
discussed  earlier,  we  chose  the  distance-correction  term 
(1.08  X  log  (A)),  whereas  the  other  studies  use  a  variety  of 
terms.  Yacoub  (1983)  and  Basham  (1969)  use  variations  of 
the  Prague  formula:  (1.66  X  log  (A))  B&th  (1967).  von 
Seggern  (1973)  used  a  slightly  smaller  distance  factor 
(0.9  x  log  (A))  than  that  of  his  later  study  which  we  use. 
The  other  three  studies  use  distance  corrections  developed 
by  Marshall  &  Basham  (1972)  and  all  are  approximately 
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0.45  magnitude  units  smaller  than  ours.  If  we  had  used  the 
Prague  formula  at  a  distance  of  40°,  our  magnitudes  would 
have  been  only  0.15  units  smaller.  The  difference  in 
distance-correction  factor  is  believed  to  be  the  primary 
cause  of  the  offset  in  magnitude  between  their  results  and 
ours.  Our  distance  correction  is  also  dependent  upon  the 
generic  path  structure  chosen  to  generate  the  40°  synthetics. 

These  three  studies,  as  well  as  that  of  Basham  (1969),  use 
mostly,  if  not  all,  data  recorded  at  Canadian  stations;  thus 
their  networks  have  strong  azimuthal  and  distance  biases  as 
well,  which  may  also  affect  magnitude  measurements.  It 
should  be  noted  that  the  method  described  in  this  study  to 
calculate  Ms  also  is  based  upon  a  theoretical  network 
average  Ms,  so  it  will  have  a  bias  attached  to  it  which  is 
dependent  upon  the  network  used.  This  network  bias  may 
be  responsible  for  part  of  the  offset,  as  well.  Our  network 
does  have  considerably  better  azimuthal  coverage  than  these 
other  studies,  so  that  tectonic-release  effects  upon  the 
long-period  radiation,  assuming  strike-slip  faulting,  should 
be  mitigated,  thus  giving  more  accurate  Ms  measurements. 

A  significant  difference  between  our  Ms  calculations  and 
those  of  the  other  studies  is  that  we  include  data  from 
close-in  stations.  Since  the  40°  synthetics  used  to  calculate 
Ms  travel  further  along  an  arbitrary  path  model  for  these 
nearer  stations,  it  is  important  to  consider  whether  or  not 
our  Ms  values  have  some  functional  dependence  upon 
distance.  Fig.  10(a)  plots  relative  station  Ms  versus  distance 
for  the  entire  data  set.  No  apparent  distance  dependence  is 
observed.  We  also  examined  this  relation  for  individual 
events  and  found  the  evidence  more  compelling  that  there  is 
no  distance  dependence  for  the  Ms  values,  which  makes  this 
Ms  method  very  attractive,  particularly  for  small  events,  for 
which  Rayleigh  wave  amplitudes  are  measurable  only  at 
near  distances,  since  there  will  be  no  bias  in  magnitude 
values  between  large  and  small  events.  Fig.  10(b)  shows  the 
relative  station  Ms  versus  azimuth.  There  is  some  variation 
with  azimuth.  This  is  to  be  expected  for  we  do  not  take  into 
account  tectonic  release  in  our  Ms  calculations.  Azimuthal 
variations  in  propagation  paths,  caused,  perhaps,  by 
different  tectonic  regions  may  also  contribute  to  this  effect. 
The  potential  bias  due  to  azimuthal  averaging  of  our 
network  and  those  of  previous  studies  is  discussed  in  the 
sequel  paper  as  a  function  of  the  orientation  and  strength  of 
tectonic  release  at  NTS. 

As  pointed  out  by  an  anonymous  reviewer,  there  is  no 
physical  significance  to  the  magnitude-offset  constant,  and  it 
only  matters  for  comparison  with  other  studies.  It  was 
suggested  that  we  could  renormalize  our  Ms  values,  i.e. 
change  the  constant  in  the  our  modified  Ms  formula,  eq.  (9), 
so  that  they  then  could  be  used  together  with  other  Ms 
measurements,  e.g.  on  an  mb  versus  Ms  plot  with 
earthquakes  where  the  earthquake  is  determined  by 
standard  means.  Unfortunately,  very  few  explosion  studies 
have  been  made  using  standard  means.  We  recommend  the 
use  of  explosion  moments  for  that  purpose  since  they  have  a 
physical  interpretation  and  have  very  few  arbitrary 
assumptions  in  their  determination.  In  this  paper,  we  prefer 
to  give  the  regression-conversion  constants  between  the 
various  explosion  studies.  They  can  be  found  in  the  labels  on 
Fig.  9  but  for  convenience,  we  give  them  here.  In  order  to 
convert  our  values  to  those  of  Basham  (1969),  add  0.45; 
Marshall  &  Basham  (1972),  subtract  0.49;  Basham  &  Horner 


Delta 


Relative  Ms  vs.  Azimuth 


Figure  10.  Relative  Ms  (individual  station-network  average)  versus 
distance  (top  figure)  and  azimuth  (bottom  figure).  Ms  values  do  not 
appear  to  be  a  function  of  distance.  There  also  is  no  apparent 
functional  relationship  between  azimuth  and  Ms ;  some  azimuths  are 
not  covered,  however. 

(1973),  subtract  0.45;  von  Seggern  (1973),  add  0.28;  Marshall 
et  aL  (1979),  subtract  0.43;  and  Yacoub  (1983),  subtract  0.01. 
As  an  example,  if  we  modify  our  regression  relation 
between  Ms  and  log  M0  for  NTS  (Fig.  4)  to  Marshall  et  aL 
(1979)  Ms  values  using  the  above,  we  obtain  Ms  = 
log  M0~  11.85.  This  is  in  fair  agreement  with  the 
observation  by  Stevens  &  McLaughlin  (1986)  that 
Ms  =  log  M0  - 12.00  fit  well  for  NTS,  Novaya  Zemla  and 
Amchitka  explosions  using  Marshall  et  aL  (1979)  Ms  values 
for  NTS  and  Amchitka.  Only  six  NTS  events  were  on  the 
figure  showing  this.  Later  in  the  report,  15  NTS  explosions, 
with  11  Marshall,  Lilwall  &  Farthing  (1986)  and  4  NEIC 
determinations  of  Ms,  were  plotted  versus  log  M0  and  these 
were  in  much  better  agreement  with  our  converted  constant 
of  11.85  than  their  three-test-site  value  of  12.00.  In  Stevens 
(1986),  it  was  found  that  11.86  fitted  his  East  Kazakh 
explosion  moments  using  Sykes  &  Cifuentes  (1984)  Ms 
values  for  the  10  events  they  had  in  common. 

Table  2  lists  the  final  mixed-path-1,  path-corrected  Ms 
values  for  the  190  events  of  this  study.  The  first  column  lists 
the  number  of  stations  recording  the  event.  Next  the 
surface-wave  magnitude  and  associated  error  for  the  event 
are  given.  This  is  followed  by  a  three-letter  shot-information 
code.  The  first  letter  denotes  its  geographic  location:  Yucca 
(Y),  Pahute  (P),  Rainer  (R),  or  Climax  Stock  (C).  The 
second  is  whether  its  shot  depth  was  above  (A)  or  below 
(B)  the  water  table.  The  last  letter  describes  the  shot-site 
rock  as  tuff  (T),  rhyolite  (R),  granite  (G),  or  alluvium  (A). 
An  underscore  means  that  the  information  is  not  known. 
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Table  2.  Event  information. 


No. 

Sta. 

Ms 

Ms 

s.e. 

Shot 

Info. 

Event 

Julian 

Date 

No. 

Sta. 

Ms 

Ms 

s.e. 

Shot 

Info. 

Event 

Julian 

Date 

1 

2.49 

R 

Rainier 

57262 

5 

3.86 

0.04 

YAT 

Yard 

67250 

1 

3.25 

- 

R 

Logan 

58289 

1 

2.08 

- 

YAA 

Marvel 

67264 

1 

3.73 

- 

R 

Blanca 

58303 

9 

3.64 

0.05 

YBT 

Cobbler 

67312 

1 

2.96 

- 

NBG 

Hardhat 

62046 

10 

5.40 

0.06 

NBT 

Faultless 

68019 

1 

2.89 

- 

YAA 

Dormouseprime 

62095 

6 

3.53 

0.05 

RAT 

Dorsalfin 

68060 

1 

3.57 

- 

YA 

Aardvark 

62132 

4 

2.85 

0.06 

N _ 

Buggyl 

68072 

1 

3.12 

- 

YAA 

Haymaker 

62178 

3 

5.80 

0.01 

PBR 

Boxcar 

68117 

1 

3.69 

- 

YAA 

Sedan 

62187 

10 

4.75 

0.09 

PAT 

Rickey 

68167 

1 

3.66 

- 

YA 

Mississippi 

62278 

11 

4.48 

0.06 

PAR 

Chateaugay 

68180 

2 

4.85 

0.22 

YBT 

Bilby 

63256 

9 

3.43 

0.08 

RAT 

Hudsonseal 

68268 

1 

3.77 

- 

RAT 

Clearwater 

63289 

1 

3.14 

- 

YAA 

Crew 

68309 

2 

3.07 

0.28 

YAL 

Handcar 

64310 

4 

3.66 

0.10 

PAT 

Schooner 

68343 

2 

2.54 

0.08 

YAA 

Merlin 

65047 

4 

5.92 

0.08 

PBT 

Benham 

68354 

1 

2.82 

NAA 

Wishbone 

65049 

11 

4.03 

0.05 

RAT 

Wineskin 

69015 

12 

3.97 

0.06 

YBT 

Wagtail 

65062 

1 

3.20 

- 

RAT 

Cypress 

69043 

6 

3.96 

0.07 

YAT 

Cup 

65085 

15 

4.06 

0.04 

YBT 

Blenton 

69120 

1 

2.59 

- 

PAR 

Palanquin 

65104 

9 

5.58 

0.05 

PBT 

Jorum 

69259 

4 

2.40 

0.06 

PBT 

Buteo 

65132 

9 

4.68 

0.05 

PAR 

Pipkin 

69281 

4 

3.08 

0.09 

NAA 

Dilutedwaters 

65167 

3 

2.11 

0.06 

YAT 

Cruet 

69302 

7 

3.66 

0.05 

YAT 

Charcoal 

65253 

3 

2.53 

0.12 

YAT 

Pod 

69302 

5 

3.81 

0.09 

YBT 

Lampblack 

66018 

4 

4.36 

0.14 

YBT 

Calabash 

69302 

8 

3.91 

0.12 

PBT 

Rex 

66055 

2 

3.17 

0.02 

R _ 

Dieseltrain 

69339 

13 

3.96 

0.07 

PAR 

Duryea 

66104 

2 

3.34 

0.05 

RAT 

Dianamist 

70042 

5 

3.47 

0.11 

NAT 

Pinstripe 

66115 

2 

3.46 

0.05 

YAT 

Cumarin 

70056 

3 

2.81 

0.05 

YAA 

Cyclamen 

66125 

2 

3.22 

0.18 

YAA 

Yannigan 

70057 

11 

4.01 

0.10 

PAR 

Chartreuse 

66126 

3 

2.14 

0.24 

YAT 

Cyathus 

70065 

15 

4.34 

0.04 

YBT 

Piranha 

66133 

3 

1.66 

0.20 

YAT 

Arabis 

70065 

6 

3.50 

0.06 

YAT 

Discus  thrower 

66147 

2 

2.22 

0.27 

YAA 

Jal 

70078 

26 

4.27 

0.04 

CBG 

Piledriver 

66153 

17 

4.29 

0.05 

YBT 

Shaper 

70082 

18 

4.29 

0.03 

YBT 

Tan 

66154 

8 

5.56 

0.08 

PBT 

Handley 

70085 

3 

2.70 

0.13 

YAA 

Vulcan 

66176 

5 

3.35 

0.13 

RAT 

Mintleaf 

70125 

3 

5.06 

0.15 

PBR 

Halfbeak 

66181 

7 

3.76 

0.08 

YAT 

Cornice 

70135 

3 

5.37 

0.11 

PBT 

Greeley 

66354 

13 

3.54 

0.06 

YAT 

Morrones 

70141 

2 

2.39 

0.05 

YAA 

Ward 

67039 

2 

1.79 

0.33 

YAT 

Manzanas 

70141 

1 

2.38 

- 

YAA 

Persimmon 

67053 

5 

3.12 

0.06 

RAT 

Hudsonmoon 

70146 

5 

4.55 

0.09 

YAA 

Agile 

67054 

16 

4.12 

0.07 

YAT 

Flask 

70146 

7 

4.96 

0.06 

YBT 

Commodore 

67140 

2 

2.40 

0.10 

YAA 

Embudo 

71167 

10 

4.64 

0.07 

PAR 

Knickerbocker 

67146 

4 

3.01 

0.09 

YAT 

Laguna 

71174 

3 

3.46 

0.05 

RAT 

Midimist 

67177 

4 

3.21 

0.02 

YAT 

Harebell 

71175 

2 

3.03 

0.14 

RAT 

Doormist 

67243 

10 

3.16 

0.06 

RAT 

Camphor 

71180 

15 

4.13 

0.03 

YBT 

Miniata 

71189 

11 

4.40 

0.05 

YAT 

Baseball 

81015 

21 

3.63 

0.06 

YBT 

Algodones 

71230 

10 

4.18 

0.07 

YBT 

Rousanne 

81316 

2 

2.59 

0.06 

YAT 

Pedemal 

71272 

9 

4.67 

0.05 

YBT 

Jornada 

82028 

4 

2.59 

0.04 

YAT 

Cathay 

71281 

10 

4.52 

0.06 

PBR 

Molbo 

82043 

4 

2.29 

0.17 

YAA 

Longchamps 

72110 

8 

4.41 

0.08 

PAR 

Hosta 

82043 

7 

3.44 

0.10 

RAT 

Mistynorth 

72123 

4 

2.99 

0.10 

YAT 

Tenaja 

82107 

5 

3.27 

0.06 

YBT 

Monero 

72140 

6 

4.44 

0.05 

PAT 

Gibne 

82115 

7 

3.35 

0.06 

RBT 

Diamondsculls 

72202 

15 

4.29 

0.04 

YBT 

Bouschet 

82127 

1 

2.50 

- 

YA_ 

Delphinium 

72270 

9 

4.52 

0.05 

PAR 

Nebbiolo 

82175 

12 

4.09 

0.04 

YBT 

Miera 

73067 

9 

3.04 

0.10 

YAT 

Monterey 

82210 

22 

4.09 

0.03 

YBT 

Starwort 

73116 

9 

4.73 

0.05 

YBT 

Atrisco 

82217 

8 

3.35 

0.08 

RAT 

Didoqueen 

73156 

8 

3.38 

0.05 

RAT 

Huronlanding 

82266 

5 

5.10 

0.12 

PBR 

Almendro 

73157 

5 

3.54 

0.11 

RAT 

Frisco 

82266 

15 

4.35 

0.05 

YBT 

Latir 

74058 

2 

2.49 

0.29 

YAA 

Seyval 

82316 

8 

3.43 

0.08 

RAT 

Mingblade 

74170 

6 

3.13 

0.09 

YAA 

Manteca 

82344 

20 

4.59 

0.05 

YBT 

Escabosa 

74191 

2 

1.80 

0.01 

YAA 

Cerro 

82245 

13 

3.96 

0.05 

YBT 

Stanyan 

74269 

3 

2.52 

0.08 

YBT 

Borrego 

82272 
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No. 

Sta. 

Ms 

Shot 

Julian 

Info. 

Event 

Date 

Shot 

Julian 

Info. 

Event 

Date 

The  final  two  columns  give  the  name  and  Julian  data  of  each 
event,  respectively.  The  events  are  listed  in  chronological 
order.  The  error  listed  is  the  standard  deviation  of  the  mean 
or  standard  error  and  the  standard  deviation  of  the 
observations  can  be  obtained  by  multiplying  by  the  square 
root  of  the  number  of  observations  or  stations  used.  The 
average  standard  deviation  for  the  fist  is  0.15.  It  is 
questionable  as  to  what  that  statistic  means,  however,  since 
there  are  events  in  the  table  with  Ms  determined  with  only 
one  station  or  infinite  standard  deviation,  and  these  were 
obviously  not  used  in  determining  the  average.  There  are 
not  a  lot  of  NTS  studies  that  give  the  error  for  Ms  but  the 
standard  deviations  or  the  more  frequently  used  standard 
error  of  path-corrected  log  M0,  which  is  very  similar  to  our 
Ms  method,  for  NTS  can  be  found  in  Stevens  (1986)  and 
Stevens  &  McLaughlin  (1988).  Stevens  (1986)  notes  that  the 
network  standard  deviations  in  (log)  moment  for  his  study 


of  40  NTS  events  are  quite  small,  about  0.1,  and  that  even 
for  recent  NTS  explosions,  which  included  data  from  several 
distant  SRO  stations,  the  standard  deviations  are  only  0.15. 

To  determine  the  portability  of  this  Ms  calculation 
method  the  events  need  to  be  separated  into  groups  based 
on  their  source  regions  and  then  compared,  one  group  to 
another,  in  order  to  see  if  there  are  systematic  differences  in 
Ms  values  relative  to  any  other  magnitude  scale.  Three  main 
geographic  source  regions  comprise  the  event  data  set: 
Pahute  Mesa,  Rainier  Mesa  and  Yucca  Flat. 

Whether  or  not  a  shot  occurs  within  saturated  material  is 
another  criterion  by  which  to  separate  events  in  order  to 
look  for  systematic  differences  in  Ms  values.  Other  studies 
have  found  significant  seismic-coupling  differences  between 
explosions  detonated  above  and  below  the  water  table 
(Marshall  et  al  1979;  Gupta  et  al  1989;  Vergino  &  Mensing 
1989).  It  is  important  to  quantify  this  seismic-coupling  effect. 
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Ms  vs.  mb,  NTS 


mb 

(a)  Ms  =1.61  x  mb -4.54 
(s.d.=0.26) 


Ms  vs.  mb,  NTS  bwt 


Ms  vs.  mb,  NTS  awt 


mb 

(e)  Ms  =1.42  x  mb -3.59 
(s.d.=0.27) 


Ms  vs.  mb,  NTS 


(s.d.=0.26) 


Ms  vs.  mb,  NTS  bwt 


Ms  vs.  mb,  NTS  awt 


(s.d.=0.27) 


Figure  11.  Ms  regressed  versus  Lilwall  mb  for  all  NTS  events.  The  data  are  also  separated  with  respect  to  shot  point  being  above  or  below  the 
water  table.  The  bottom  figures  are  constrained  least-square  regressions  assuming  a  slope  of  1.5. 


Reviewing  Fig.  6,  it  is  also  apparent  that  shots  fired-off 
below  the  water  table  have  a  larger  body-wave  magnitude 
than  those  detonated  above  the  water  table. 

Figure  11(a)  shows  the  Ms:mb  (Lilwall)  relationship  for 
all  NTS  events.  The  surface-wave  magnitudes  were  all 
calculated  using  mixed-path  Green’s  functions  (with  the 
RSSD-1  structure  for  the  generic  path  part)  and  path 
corrections.  Figs  11(b)  and  (c)  divide  the  data  populations 
into  above  and  below  the  water  table,  respectively;  shots  for 
which  water-table  information  was  not  available  were  left 
out.  Although  all  but  one  Rainier  Mesa  events  were 
detonated  above  the  water  level,  we  found  that  their 
coupling  (Ms  versus  log  yield)  was  diagnostic  of  explosions 
detonated  below  the  water  table.  Taylor  (1983)  notes  that 
Rainier  Mesa  sports  a  perched  aquifer.  We  believe  that  the 
Rainier  Mesa  events  are  detonated  within  this  zone,  hence 
they  are  assumed  to  be  well-coupled  events,  i.e.  the  pore 
space  of  the  shot  medium  is  filled  with  water  and  thus 
pore-space  crushing  will  not  be  a  strong  effect. 

The  bottom  three  figures  (lld,e  and  f)  plot  the  same 
data,  but  a  constrained  least-squares  fit  was  performed  with 
the  slope  =  1.50.  The  offset  in  curves  between  events 
detonated  above  and  below  the  water  table  is  0.08.  This 
amount  is  within  the  scatter  of  the  data  (i.e.  statistically 
insignificant),  but  it  would  appear  that  shot-medium 


coupling  effects  associated  with  pore-filling  phenomena  are 
similar  for  surface  waves  and  P  waves. 

Figures  15(a)-(c)  are  Ms  versus  log-yield  plots  analogous 
to  Figs  ll(a)-(c).  It  is  important  to  note  that  the  individual 
explosion  variances  are  about  the  same  size  for  the  entire 
range  of  yields,  so  that  our  predicted  yield  values  for  small 
events  should  be  as  accurate  as  for  the  larger  events.  The 
slope  of  the  Ms  and  log-yield  scaling-relation  curves  was 
found  to  be  near  unity  for  all  populations.  Assuming  the 
scaling  relationship  has  a  slope  of  1,  BWT  shots  couple  more 
strongly  than  AWT  shots  by  0.52  units— a  substantial 
amount;  for  mb  -yield  scaling  the  coupling  effect  found  in 
this  study  is  0.28  units.  This  coupling  factor  depends  on  the 
slope  of  the  scaling  curve  and  has  been  found  to  be  as  large 
as  0.7  to  0.9  throughout  the  literature.  For  individual  source 
regions,  the  offset  in  the  M5-yield  and  mb-yit\d  scaling 
curves  for  shots  fired  above  and  below  the  water  table  vary 
slightly  from  these  values  determined  from  the  entire  data 
set.  There  is  some  scatter  in  the  data  which  is  not  surprising 
considering  the  diversity  of  the  sampled  populations. 
However  the  best-fitting  Ms-mb  curves  are  well  constrained, 
for  the  population  covers  a  wide  range  of  magnitudes. 

The  various  Ms:mb  relationships  for  Yucca  events  are 
shown  in  Figs  12(a)— (c).  The  scatter  in  the  data  is  reduced 
by  25  per  cent  over  that  of  the  general  population  (Figs 
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Ms  vs.  mb,  Yucca 


(a)  Ms  =1.29*  mb -2.93 
(s.d.=0. 1 8) 


Ms  vs.  mb,  Yucca  bwt 
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Figure  12.  Ms  regressed  versus  Lilwall  mh  for  Yucca  events.  The  data  are  also  separated  with  respect  to  shot  point  being  above  or  below  the 
water  table.  The  bottom  figures  are  constrained  least-square  regressions  assuming  a  slope  of  1.5. 


lla-c).  The  free-slope  regression  curves  are  not  as  well 
constrained  as  those  of  Fig.  11  because  the  yield  range  for 
Yucca  explosions  is  smaller  than  that  of  the  entire  data  set. 
Fixing  the  slope  to  1.5  (Figs  12d-f)  leads  to  AWT  shots 
coupling  0.13  Ms  units  more  strongly  than  BWT  shots  for  a 
given  mb.  The  error  in  the  fit  to  the  curve  is  larger  than  this 
variation,  so  it  is  not  a  statistically  significant  result.  It  would 
appear  the  pore-filling  coupling  affects  surface  wave  and 
body  waves  similarly.  When  the  Yucca  Ms  data  are 
regressed  with  respect  to  log  yield,  as  shown  in  Figs 
15(d)— (f),  it  is  found  that  BWT  events  couple  four  times 
more  efficiently  than  AWT  shots  (AA75=0.61).  This  is  a 
significant  amount  and  the  data  set  on  which  it  is  based  is 
more  extensive  than  that  of  the  Ms-mb  regression.  Springer 
(1966)  found  that  high  dry  porosity  (60  per  cent)  shot 
mediums  coupled  four  to  five  times  less  effectively  than  in 
saturated  alluvium.  Most  Yucca  Flat  shots  are  detonated  in 
alluvium. 

Figure  13(a)  plots  all  Pahute  event  Mfs  versus  their 
respective  mb  s.  The  scaling  relationship  is  significantly 
different  than  that  of  the  Yucca  data  above.  Comparing  the 
unconstrained  below-  and  above-water-table  curves  (Figs 
13b  and  c)  to  their  Yucca  counterparts  (Figs  12b  and  c),  it  is 
apparent  that  explosions  at  the  two  sites  do  not  display  the 
same  scaling  relationships.  One  possible  explanation  for  this 
difference  is  that  there  is  not  enough  data  to  constrain  the 


scaling  curves,  particularly  for  Yucca  BWT  and  Pahute 
AWT  events.  Another  possible  explanation  is  that  this 
scaling  relationship  difference  is  real  and  may  be  caused  by 
differences  in  the  source  medium,  source  structure  or 
tectonic-strain  release  associated  with  the  sites.  Figs 
12(d)-(f)  and  13(d)-(f)  show  constrained  (slope  =  1.5) 
regression  curves  for  the  Yucca  and  Pahute  data, 
respectively.  For  a  given  mb ,  surface-wave  magnitudes  for 
events  at  Pahute  Mesa  are  larger  than  those  at  Yucca  Flat 
by  0.39  and  0.18  log  units  for  BWT  and  AWT  shots, 
respectively.  There  is  also  an  appreciable  difference  in  the 
Ms :  log-yield  relationship  between  Yucca  and  Pahute  events 
detonated  in  water-saturated  material  (0.23  units).  The 
Pahute  data  are  plotted  in  Figs  16(d)— (f).  For  the  case  of 
events  exploded  in  dry  material  there  is  a  significant 
difference  with  Yucca  events  having  a  Ms  0.44  units  smaller 
than  Pahute  events. 

Figures  14(a)-(f)  display  the  Ms-mh  regression  curves 
and  data  for  Rainier  Mesa  events  in  combination  with  and 
without  Pahute  Mesa  data.  Figs  16(a)-(c)  are  analogous 
figures  for  the  Ms  versus  log-yield  data  to  Figs  14(a)-(c). 
Although  the  clustering  of  Rainier  data  near  ^=5.0 
causes  the  curve  to  be  poorly  constrained,  a  slope  is 
obtained  that  is  close  to  that  for  Pahute  and  Yucca  BWT 
shots.  Comparison  of  the  equations  at  the  bottom  of  Figs 
13(e)  and  14(e)  give  an  offset  of  0.47  between  Ms  estimates 
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Figure  13.  Ms  regressed  versus  Lilwall  mh  for  Pahute  events.  The  data  are  also  separated  with  respect  to  shot  point  being  above  or  below  the 
water  table.  The  bottom  figures  are  constrained  least-square  regressions  assuming  a  slope  of  1.5. 


at  Pahute  and  Rainier  (for  a  given  mh).  Either  the  Pahute 
site  is  more  efficient  at  producing  surface  waves  or  the 
Rainier  site  is  more  efficient  at  coupling  body-wave  energy. 
Rainier  events  are  tunnel  shots.  The  immediate  source 
region  ( R  <  200  m)  may  behave  like  an  asymmetric  cavity, 
resulting  in  a  source  that  is  non-isotropic  (Zhao  & 
Harkrider  1992)  and/or  seismic  coupling  that  has  strong 
frequency  dependence.  Either  of  these  effects  may  account 
for  this  difference.  The  difference  in  the  Af^-log-yield 
scaling  relationship  is  somewhat  less  (0.31  units),  implying 
that  Rainier  more  efficiently  couples  short-period  energy 
than  Pahute. 

Comparing  Fig.  11(b)  with  15(b),  we  see  that  regressing 
against  log  yield  has  reduced  the  regression  standard 
deviation  from  0.25  to  0.18  for  NTS  BWT  events  plus 
Rainier  events.  This  is  probably  caused  by  the  reduced 
standard  deviation  of  the  Pahute  BWT  events,  0.18  to  0.13 
(Figs  13b  and  16e)  and  the  combined  Pahute  BWT  plus 
Rainier  events,  0.22  to  0.16  (Fig  14b  and  16b).  The  opposite 
is  true  for  the  ‘all’  NTS  events.  This  is  because  of  the 
increase  in  standard  deviation  for  Yucca  AWT  events  when 
regressed  against  yield  due  primarily  to  the  inclusion  of  a 
lot  of  small-yield  Yucca  events  for  which  we  did  not  have 
mb’ s. 


CONCLUSION 

In  this  study  we  have  determined  surface-wave  magnitudes 
for  small  as  well  as  large  underground  nuclear  explosions. 
Our  technique  allows  us  to  include  smaller  events  in  a 
consistant  manner  with  the  historic  set  of  large  events  for 
which  surface-wave  magnitudes  have  been  determined  by 
classical  means.  Thus  it  was  not  primarily  an  attempt  to 
improve  Ms  for  large  events  but  to  extend  it  to  lower-sized 
events  by  including  regional  stations  not  usually  used  in  NTS 
Ms  determinations.  The  Ms  formula  used  was  one  that  had 
previously  been  found  appropriate  for  NTS  explosions  and 
not  the  Prague  formula  used  for  earthquakes.  In  the  process 
of  making  these  determinations,  we  also  calculated  station 
and  network  moments.  Since  the  assumptions  and 
corrections  used  in  the  moment  determinations  were  more 
straightforward,  we  feel  that  future  estimates  of  surface- 
wave  energy  should  be  moment  and  until  that  time  we  feel 
that  our  technique  is  best  for  including  small  events  in  the 
historical  surface-wave  magnitude  data  base.  The  moment 
values  will  be  given  in  the  sequel  paper. 

The  method  we  have  described  to  calculate  surface-wave 
magnitudes  allows  the  measurement  of  Ms  for  nuclear 
explosions  over  a  wider  magnitude  distribution  than  was 
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Figure  14.  Ms  regressed  versus  Lilwall  mb  for  Pahute  and  Rainier  events.  Regression  of  Rainier  data  alone  (middle  figures),  all  Rainier  and 
Pahute  events  (left  figures),  and  all  Rainier  events  with  Pahute  shots  below  the  water  table  (right  figures)  are  shown.  The  bottom  figures  are 
constrained  least-square  regressions  assuming  a  slope  of  1.5. 


previously  possible.  These  Ms  values  scale  consistently 
(within  a  constant  factor)  with  other  seismic-magnitude 
scales.  Using  our  technique,  it  is  now  possible  to  use 
near-regional  (A  <  8°)  long-period  records,  as  well  as  more 
conventional  far-regional  (A  <  15°)  and  teleseismic  observa¬ 
tions,  in  order  to  measure  surface-wave  magnitudes.  As  it  is 
a  time-domain  measurement,  it  is  easy  to  calculate  Ms  from 
historical  analogue  waveforms,  since  it  is  only  necessary  to 
measure  the  peak-to-peak  Rayleigh  wave  amplitude. 

This  Ms  method  is  very  useful  for  quantifying  small 
explosions,  because  time-domain  magnitude  measurements 
of  regional  waveforms  lowers  the  effective  magnitude 
threshold.  Small  events,  for  which  teleseismic  surface  waves 
are  not  observed,  can  now  be  analysed  with  regional 
surface-wave  data,  thus  lowering  the  effective  measuring  Ms 
threshold.  Fig.  17  illustrates  this  point,  showing  unrotated 
three-component  data  for  FLO  YD  ADA  (8/15/91,  mb  =  4.2) 
detonated  at  Yucca  Flat  and  recorded  by  three 
TERRAscope  stations  convolved  with  a  Press-Ewing  30-90 
response.  The  source-to-receiver  distances  are  between  210 
and  390  km.  The  maximum  peak-to-peak  amplitudes  are 
quite  small  (<0.5  mm).  On  the  actual  analogue  instrument 
it  would  not  be  possible  to  measure  the  surface-wave 
amplitude.  Because  of  the  low  signal-to-noise  ratio  a 


spectral  moment  would  be  of  dubious  value.  However,  the 
Ms  and  A/0(PPA)  methods  described  in  this  paper  would 
furnish  an  accurate  surface-wave  magnitude  with  which  to 
estimate  its  yield. 

These  small  surface-wave  magnitudes,  based  on  near- 
regional  data  would  also  be  of  considerable  value  for 
discrimination  methods  that  make  use  of  the  difference 
between  the  long-period  and  short-period  spectral  content 
of  earthquakes  and  explosions,  for  it  is  possible  to  lower  the 
discrimination  threshold  using  such  data. 

The  increase  in  nearer  observations  has  several  other 
advantages.  Station-network  coverage  is  enhanced  in  terms 
of  overall  numbers  as  well  as  in  azimuthal  coverage.  In  this 
study  stations  a  few  hundred  kilometers  away  from  NTS  in 
the  south-western  U.S.  were  included  in  the  network  that 
otherwise  would  have  no  coverage  to  the  west  or 
south-west.  These  improvements  make  the  network  Mf  s 
more  stable  and  statistically  robust.  Potential  monitoring 
areas  may  well  have  similar  geographical  constraints 
requiring  the  use  of  near-regional  (A<8°)  seismic  data. 
Also,  the  effect  of  inaccuracies  in  estimating  Q  are  negligible 
for  very  near-regional  recordings. 

From  the  results  obtained  with  the  data  set  used  here, 
there  do  appear  to  be  significant  differences  in  seismic 
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Figure  15.  Ms  regressed  versus  log  yield  for  all  NTS  events  (top  figures)  and  for  Yucca  events  (bottom  figures).  Event  populations  have  also 
been  grouped  with  respect  to  shot-point  water-table  location. 
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Figure  16.  Ms  regressed  versus  log  yield  for  Pahute  and  Rainier  events  (top  figures)  and  for  Pahute  events  alone  (bottom  figures).  Event 
populations  have  also  been  grouped  with  respect  to  shot-point  water-table  location. 


Determining  surface -wave  magnitudes  497 

Floydada  (Yucca  Flat,  Y  <  10  kt) 


Figure  17.  TERRAscope  streckeisen  recordings  of  an  NTS  explosion  Floydada  at  Yucca  Flat  on  8/15/91  with  an  estimated  yield  of  <10  kt. 
(mb  =  4.2,  Ml  ~  4.0,  and  logM0  =  14.16  N-M).  The  broad-band  records  have  been  convolved  with  a  Press-Ewing  30-90  instrument  response. 
All  four  stations  record  the  surface  wave  train  well  enough  to  measure  the  Airy-phase  peak  amplitudes.  Records  from  an  actual  30-90 
long-period  instrument  would  be  unusable. 


coupling  between  NTS  subsites,  with  events  at  Pahute  Mesa 
producing  larger  surface-wave  magnitudes  for  a  given  mb 
than  at  Rainier  Mesa  or  Yucca  Flat.  For  well-coupled  events 
this  discrepancy  is  largest  for  Rainier  Mesa  events.  Ms 
values  at  Yucca  Flat  tend  to  be  larger  than  those  at  Rainier 
Mesa  by  0.08  magnitude  units  for  a  given  mb.  There  also 
appears  to  be  some  difference  in  waveforms  between  events 
of  these  two  source  regions.  Pahute  Mesa  events  are  0.39 
magnitude  units  larger  than  those  at  Yucca  flat  for 
explosions  set  off  below  the  water  table  and  with  the  same 
mb. 

Although  Lg  measurements  with  a  calibration  shot  give 
more  accurate  estimates  of  explosion  yields,  there  may  be 
cases  where  Lg  ‘blockage’  caused  by  strong  lateral  variations 
in  the  propagation  path  may  occur,  and  one  must  use  other 
methods,  such  as  surface-wave  magnitudes,  to  estimate 
yields  or  for  discriminating  the  event. 
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Summary 

Seismic  moments  of  Nevada  Test  Site  (NTS)  explosions  are  determined  from  regional  sur¬ 
face  wave  spectra.  Two  methods  axe  used.  In  one  the  moment  is  solved  for  assuming  only 
an  explosive  source,  or  average  scalar  moment;  in  the  other  a  joint  inversion  for  an  isotropic 
(explosive)  source  plus  a  constrained  double  couple  moment  component  representing  tectonic 
strain  release  (TSR).  Although  the  general  moment  tensor  solution  to  this  joint  inversion 
problem  is  non-unique,  if  some  assumptions  are  made  concerning  the  non-isotropic  moment 
components,  then  the  remaining  source  parameters  can  be  solved  by  a  linear  least-squares 
inversion  scheme.  We  examine  the  errors  in  determining  the  isotropic  moment  component 
(Mi)  by  this  latter  method  of  constrained  linear  inversion  solutions  in  a  canonical  study  using 
a  theoretical  network  of  long-period  (6-60  sec.)  surface  wave  data.  The  network  azimuthal 
coverage  was  chosen  to  represent  that  of  a  long-period  North  American  super-network  of  55 
stations  used  for  the  actual  NTS  events.  We  compare  these  errors  in  moment  estimate  to 
those  obtained  from  surface  wave  magnitude  (Ms)  and  spectral  scalar  moment  (M0)  mea¬ 
surements  for  the  same  surface  wave  observations.  For  a  ratio  of  M(exp;)/M(e9)  less  than 
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1.0  we  find  that  the  inverted  Mj  solution  is  a  much  better  estimate  of  the  actual  isotropic 
moment  than  either  Ms  or  Mo,  and  the  standard  deviation  in  this  estimate  is  substantially 
less  than  that  using  the  other  two  methods  for  the  great  majority  of  isotropic  source  -f 
double  couple  sources.  Even  when  the  inversion  constraints  are  off  in  dip  and  rake  each  by 
30°,  the  mis-estimate  of  the  isotropic  moment  is  less  than  35  percent  of  the  actual  value.  In 
the  case  of  a  vertical  strike-slip  fault,  the  inverted  isotropic  moment  solution  which  assumes 
this  fault  orientation  is  exact  to  three  figures,  whereas  Ms  and  Mo  under-estimate  the  mo¬ 
ment  by  45  percent  and  32  percent,  respectively  because  of  uneven  azimuthal  coverage.  This 
moment  tensor  inversion  method  is  applied  to  determine  the  isotropic  source  for  109  NTS 
underground  explosions  using  vertical  and  tangential  component  surface  wave  data  from 
this  regional  network.  We  also  calculate  Ms  and  M0  for  these  same  events  and  compare 
the  results.  Isotropic  source  errors  are  smallest  using  the  spectral  domain  inversion  method. 
However,  this  spectral  domain  method  cannot  attain  as  low  a  magnitude  threshold  as  the 
time  domain  moment  or  Ms  method.  The  extensive  moment  data  set  analyzed  here  were 
combined  with  larger  yield  explosions  from  prior  moment  studies  to  create  a  comprehen¬ 
sive  data  set  with  which  to  obtain  conclusive,  well-constrained  long-period  explosion  source 
scaling  relationships  at  the  separate  NTS  sub-sites. 

Introduction 

The  purpose  of  this  work  is  to  determine  the  precision  of  seismic  moment  (Mo)  estimates 
of  underground  nuclear  explosions  determined  from  spectral  measurements  of  intermediate- 
period  (6  to  60  sec.),  regional  surface  waves.  Seismic  moment  is  a  measure  of  low  frequency 
seismic  source  spectrum.  Besides  being  useful  as  a  quantifier  of  earthquake  sources,  it  also  can 
be  used  to  estimate  seismic  yields  of  explosions  and  employed  in  long-period:short-period 
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discriminants  (  Woods  et  al. ,  1993).  Woods  &  Harkrider  (1994)  dealt  with  time-domain 
moment  and  magnitude  (Ms)  measurements  of  regional  surface  waves.  A  criticism  of  time 
domain  amplitude  measurements  of  regional  surface  waves  is  that  lateral  variations  in  the 
crustal  structure  cause  multi-pathing  and  scattering  which  affect  the  peak  amplitude  of  the 
Airy  Phase,  particularly  for  paths  shorter  than  1000  km  for  which  the  longer-period  wave- 
train  is  not  well  dispersed.  Spectral  moment  M(u)  measurements  reduce  this  problem  by 
averaging  over  the  frequency  band  of  interest.  They  are  also  insensitive  to  spall  effects  which 
are  significant  primarily  in  this  0.2  to  2  Hz  frequency  bandwidth  (Taylor  &  Randall,  1989), 
since  the  spectral  moment  is  calculated  between  0.0167  and  0.167  Hz. 

Unlike  short-period  body  wave  amplitudes,  surface  waves  are  relatively  insensitive  to  source 
region  attenuation  effects  which  give  rise  to  mj,  biases  in  yield-magnitude  relationships  be¬ 
tween  inter-regional  testing  sites. 

Previous  moment  studies  of  explosions  have  relied  on  distant  regional  (D  >  1000  km)  and 
teleseismic  surface  wave  measurements  and  as  such  only  the  larger  explosions  (mj,  >  5.5) 
could  be  examined.  Stevens  (1986)  determined  scalar  moments,  assuming  a  pure  isotropic 
(explosive)  source,  for  Nevada  Test  Site  (NTS)  and  East  Kazakhstan  Test  Site  events.  Given 
&  Mellman  (1986)  inverted  fundamental-mode  Rayleigh  and  Love  wave  spectral  data  to  ob¬ 
tain  the  isotropic  component  and  double-couple  component  (due  to  tectonic  strain  release) 
of  the  source  moment  tensor  for  virtually  the  same  data  set.  Both  studies  relied  on  seis¬ 
mic  networks  with  poor  azimuthal  coverage  and  included  mixed  continental-oceanic  paths 
which,  as  discussed  in  Woods  &  Harkrider  (1995),  can  significantly  alter  surface  wave-train 
waveforms  and  spectral  content.  Although  both  studies  made  use  of  uniquely  derived  path 
structures  to  account  for  propagation  effects,  they  also  made  use  of  empirically  derived  ”sta- 
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tion  corrections”  directly  or  indirectly,  but  which  are,  in  fact,  path  corrections  to  reduce  the 
scatter  in  the  network  averaged  moment  caused  by  apparent  errors  in  modeling  the  paths. 
These  two  studies  also  relied  heavily  on  hand-digitized  data  for  events  in  the  70  s;  the  digital 
data  available  was,  for  the  most  part,  from  stations  with  oceanic  paths.  One  would  rather 
use  digitized  data  from  pure  continental  paths,  since  waveforms  should  be  more  reliable 
and  the  signal  to  noise  ratio  (SNR)  should  be  better,  particularly  for  smaller  events.  The 
mixed-path  approximation  theory  for  synthetic  surface  wave  generation  discussed  in  Woods 
&  Harkrider  (1995)  works  best  for  paths  that  aren’t  radically  different,  so  it  is  not  clear  that 
this  approximation  would  work  well  for  ocean-continent  paths. 

Evernden  et  al.  (1971)  suggest,  that  with  high  dynamic  range  digital  seismometers,  it 
should  be  possible  to  observe  surface  waveforms  from  small  earthquakes  (mj,  «  4.0)  up  to 
6000  km  distance  away;  for  explosions  which  have  relatively  lower  level  long-period  source 
spectra,  a  higher  detection  threshold  would  be  expected. 

This  study  re- investigates  the  long-period  source  spectrum  (seismic  moment)  of  explosions 
using  surface  wave  spectral  amplitude  measurements  from  digitally  recorded,  reasonably  high 
dynamic  range,  regional  (375  km  <  5000  km)  waveforms.  Such  data  have  been  used  in  earth¬ 
quake  source  studies  (Patton  &  Zandt,  1991;  Yelin  &  Patton,  1991;  and  Thio  &  Kanamori, 
1992),  but  until  now  has  not  been  applied  to  underground  explosion  source  studies.  We 
wish  to  determine  whether  such  data  better  constrain  moment  estimates  of  explosions  and 
whether  the  observational  magnitude  threshold  for  measuring  surface  moments  can  be  low¬ 
ered.  This  study  is  confined  to  NTS  events;  however,  it  is  a  more  comprehensive  data  set 
in  that  it  includes  109  events-nearly  three  times  as  many  events  as  the  previous  studies  and 
events  from  Rainier  Mesa  are  examined  as  well.  As  in  Woods  &  Harkrider  (1995)  localized 


5 


site  effects  on  moment  estimates  can  be  studied.  The  digital  data  for  this  study  is  not  as 
comprehensive  as  that  in  Woods  &  Harkrider  (1995)  which  included  many  analog  data,  so 
that  it  will  not  be  possible  to  constrain  the  magnitude-moment  scaling  relationships  as  well 
as  the  in  time-domain  moment- Af$  moment  study. 

With  the  use  of  nearer  regional  stations,  better  azimuthal  coverage  of  NTS  is  attain¬ 
able.  This  is  important  for  obtaining  unbiased  average  scalar  moments  and/or  to  invert  for 
non-isotropic  moment  components.  Also  with  shorter  paths,  errors  in  propagation  effects, 
principally  because  of  attenuation,  Q,  will  be  smaller  than  for  the  longer  paths. 

Particular  attention  is  paid  to  error  analysis  in  this  study,  since  high  confidence  levels  in 
nuclear  monitoring  are  a  primary  concern.  Standard  deviations  in  scalar  moments  estimates 
by  Stevens  (1986)  ranged  between  12  and  200  percent  ,  with  most  being  under  25  percent. 
While  no  error  analysis  appears  to  have  been  conducted  in  the  moment  tensor  inversion 
study,  it  was  found  that  some  explosions  exhibited  radiation  patterns  which  have  double¬ 
couple  moment  components  as  large  as  50  percent  of  the  isotropic  moment  component.  The 
explosions  with  the  largest  apparent  tectonic  strain  release  moments  in  the  Given  &  Mellman 
study  also  tended  to  have  the  largest  error  estimates  in  the  Stevens  (1986)  scalar  moment 
study.  Hence  it  would  be  instructive  to  examine  how  well  correcting  for  tectonic  release 
effects  would  reduce  the  variance  in  moment  estimates. 

The  assumption  that  the  radiation  pattern  from  an  explosive  source  has  no  azimuthal 
dependence,  while  theoretically  sound,  is  not  born  out  in  practice.  Many  researchers  have 
reported  evidence  for  a  non-isotropic  component  for  explosion  sources  (Press  &  Archambeau, 
1962;  Brune  &  Pomeroy,  1963;  Toksoz  et  al .,  1971;  Toksoz  &  Kehrer,  1972;  Tsai  &  Aki, 
1972;  Wallace  et  al.,  1985)  and  has  been  attributed  predominantly  to  tectonic  strain  release 
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(Archambeau  &  Sammis,  1970)  and,  more  recently,  to  asymmetric  explosion  cavities  (Zhao 
&  Harkrider,  1992).  This  additional  source  can  have  appreciable  long-period  content,  and 
as  such  can  effect  surface  wave  radiation  patterns  (Helle  &  Rygg,  1984;  Given  &  Mellman, 
1986),  and  thus  bias  moment  estimates  if  it  is  not  accounted  for. 

The  most  general  moment  tensor  inversion  solution  for  a  joint  isotropic  +  double-couple 
source  is  non-unique  (Mendiguren,  1977;  Ekstrom  &  Richards,  1994).  And  for  very  shallow 
sources  some  off-diagonal  elements  of  the  moment  tensor  become  unresolvable  (Kanamori 
&  Given,  1981;  and  Patton,  1988).  If  one  is  to  invert  for  source  depth  and  source  time 
functions,  the  problem  becomes  non-linear  and  even  more  complicated.  In  the  inversion 
procedure  developed  herein,  several  assumptions  are  made  concerning  the  sources.  The 
tectonic  release  source  and  explosion  sources  are  both  modeled  as  step  moment  source  time 
functions  coincident  in  time  and  depth,  which  are  reasonable  assumptions  for  the  periods 
measured  (6-60  sec.).  Secondly,  as  described  later  in  the  Inversion  Method  Section,  certain 
moment  tensor  components  and  source  excitation  functions  are  assumed  to  be  negligible 
when  source  depth  is  taken  as  the  limit  h  — ►  0  (see  Given  &  Mellman,  1986,  for  details). 
The  tectonic  release  radiation  pattern  is  assumed  to  be  modeled  as  a  double-couple  source. 
Further,  the  dip  and  rake  of  the  double-couple  fault  plane  are  constrained  by  previous  studies 
of  tectonic  release  at  NTS  to  be  predominantly  vertical,  strike-slip  in  nature  with  a  sin2<t> 
radiation  pattern  (Toksoz  et  of.,  1965;  Aki  &  Tsai,  1972;  Wallace  et  ul .,  1983  and  1985, 
Lay  el  al,  1984;  Given  &  Mellman,  1986),  with  the  strike  angle  varying  between  N10°W  and 
N60°W. 

The  tectonic  release  strike  angle  ( <j> ),  moment  (Af# ),  and  the  explosion  moment  (Afj)  are 
then  inverted  from  fundamental-mode  Rayleigh  wave  and  Love  wave  spectra  in  what  is  then 


a  one-step  linear  least-squares  inversion.  The  errors  in  M/  are  compared  to  those  of  scalar 
moments  (Mo).  The  source  parameters  and  yield-magnitude  relationships  axe  compared 
between  the  three  major  NTS  sub-sites:  Pahute  Mesa,  Rainier  Mesa,  and  Yucca  Flats. 

Data  Processing  and  Analysis 

The  digital  and  digitized  data  of  Woods  &  Harkrider  (1995)  are  used  in  this  study.  Table  1 
lists,  in  chronological  order,  the  events  used  along  with  available  source  information.  In  the 
date  column  is  the  Julian  date.  The  three-letter  site  code  (3rd  column)  gives  the  NTS  sub-site 
location  (P=Pahute  Mesa,  Y=Yucca  Flat,  R=Rainier  Mesa  (Tunnel  shots),  or  C=Climax 
Stock),  whether  the  shot  was  detonated  above  (A)  or  below  (B)  the  water  table,  and  rock 
type  (T= Volcanic  Tuff,  G=Granite,  R=Rhyolite,  A=Alluvium).  When  available,  the  shot 
point  elevation  and  depth  of  burial  are  given.  Also  listed,  when  available,  are  body  wave 
magnitudes  determined  by  the  International  Seismic  Centre  (ISC)  and  Lilwall  &  Neary  (1985) 
(L),  and  local  magnitudes.  Data  for  events  before  1982  are  digitized  records,  while  all  later 
events  rely  on  digital  recordings.  In  most  cases  the  vertical  and  two  horizontal  components 
are  all  available  for  any  particular  event-station  pair,  although  for  a  small  minority  of  the 
data  pairs  only  the  vertical  or  horizontal  component  records  were  usable/existent.  Horizontal 
components  were  rotated  to  their  great  circle  path  to  obtain  the  transverse  component;  the 
radial  components  were  not  used  as  they  yield  only  redundant  P-SV  wavefield  (Rayleigh 
wave)  information  and  can  be  contaminated  with  scattered  transverse  component  energy. 
Horizontal  components  also,  in  general,  have  a  lower  SNR  than  vertical  components. 

The  SNR  for  some  observations  is  low,  particularly  for  the  smaller  events  (mj  «  5.0  and 
below),  with  long-period  (T  >  60  sec)  noise  overwhelming  the  transient  signal.  High-pass 
filtering  improved  the  time-domain  SNR  significantly,  making  positive  identification  of  Love 
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NAME 

Date 

Site 

Lat. 

MINGBLADE 

74170 

RAT 

37.21 

ESCABOSA 

74191 

YBT 

37.08 

STANYAN 

74269 

YBT 

37.13 

CABRILLO 

75066 

YBA 

37.13 

DININGCAR 

75095 

RAT 

37.19 

OBAR 

75120 

YBT 

37.11 

MEZEN 

75154 

YBT 

37.09 

HUSKYFUP 

75297 

RAT 

3722 

KEELSON 

76035 

YBT 

37.07 

MIGHTYEPIC 

76133 

RAT 

3721 

RUDDER 

76363 

YBT 

37.10 

BULKHEAD 

77117 

YBT 

37.09 

CREWLINE 

77145 

YBT 

37.09 

LOWBALL 

78193 

YBT 

37.08 

QUARGEL 

78322 

YBT 

37.13 

QUINELLA 

79039 

YBT 

37.10 

PYRAMID 

80107 

YBT 

37.10 

MINERSIRON 

80305 

RAT 

37.24 

BASEBALL 

81015 

YAT 

37.09 

JORNADA 

82028 

YBT 

37.09 

MOLBO 

82043 

PBR 

3122 

HOSTA 

82044 

PAR 

37.35 

TENAJA 

82107 

YAT 

37.02 

GEBNE 

82115 

PAT 

37.26 

KRYDDOST 

82126 

■a 

37.12 

BOUSCHET 

82127 

YBT 

37.07 

NEBBIOLO 

82175 

PAR 

3724 

MONTEREY 

82210 

YAT 

37.10 

Lon.  Elev. 


-116.21 


-116.07 


-116.08  1315 


-116.21  2265 


-116.03 


-116.04  1274 


-116.03  1285 


-116.21  2251 


-116.03 


-116.04  1252 


-116.04 


-116.08 


-116.05 


-116.03 


-11621  2239 


-116.04 


-116.05 


-116.46 


-116.01  1210 


-116.42  1964 


-116.13  1390 


-116.05  1244 


-116.07  I  1280 


4.97  0.04 


5.74  0.02  5.6 


5.52  0.03  5.0 


5.57  0.03  5.2 


4.94  0.03 


5.12  0.04  5.0 


5.71  0.02  5.6 


4.87  0.05 


5.72  0.02  5.8 


4.85 


5.69  0.02  5.5 


5.43 


5.40 


5.72  0.02 


5.71  0.02  5.2 


5.45 


4.97  0.06  4.9 


5.72 


5.92 


5.43  0.03 


5.55 


.59 


5.38  0.03  5.4 


Table  1  (c) 


NAME 


Lon.  Elev. 


ATRISCO 


82217 


HURONLANDING  82266 


BORREGO  82272 


FRISCO  82267 


MANTECA  82344 


CABRA  83085 


TORQUOISE  83104 


CROWDIE  83125 


FAHADA  83146 


DANABLU  83160 


CHANCELLOR  83244 


MTONTTEZEPHYR  83264 


TECHADO 


ROMANO 


MILAGRO 


TORTUGAS 


MUNDO 


CAFROCK 


DUORO 


KAPPELI 


CORREO 


DOLCETTO 


BRETON 


VILLITA 


EGMONT 


TERRA 


VAUGHN 


COTTAGE 


83265 


83350 


84046 


84061 


84122 


84152 


84172 


84207 


84215 


84243 


84257 


84315 


84344 


84350 


85074 


85082 


5.1 


NAME 

Date 

Site 

HERMOSA 

MISTYRAIN 

85096 

R_T 

TOWANDA 

SALUT 

85163 

PBR 

SERENA 

85206 

PAR 

PONIL 

EiSSESill 

ROQUEFORT 

85289 

YAT 

KINIBITO 

85339 

YBT 

GOLDSTONE 

85362 

P_ 

GLENCOE 

86081 

YBT 

MIGHTYOAK 

JEFFERSON 

PANAMINT 

86141 

YAT 

TAJO 

86156 

YAT 

DARWIN 

86176 

P_ 

CYBAR 

86198 

PAR 

CORNUCOPIA 

— «wi 

LABQUARK 

86273 

PAR 

BELMONT 

86289 

P_ 

GASCON 

86318 

YAT 

BODIE 

86347 

PAT 

DELAMAR 

87108 

P_ 

HARDIN 

87120 

PAT 

PANCHUELA 

87181 

YAA 

TAHOKA 

mis 

KZ1I 

LOCKNEY 

87267 

P_ 

BORATE 

87296 

ESI 

Table  1  (d) 


Lon.  Elev. 


-116.03 


-116.21 


-116.33  2112 


-116.49 


-116.44 


6.00 


-116.12 


-116.05 


-116.47 


-116.07  1260 


-116.18 


-116.44 


6.06 


-116.02  1316 


-116.50  1876 


-116.36  2044 


-116.07 


-116.31 


-116.05  1263 


-116.41  2018 


-11631  1902 


-116.42  1970 


-116.04  1206 


-116.05  1239 


-116.08  I  1321 


Depth 


-640 


-389 


-661 


-608 


-597 


-366 


-415 


-600 


-500 


-600 


-400 


-600 


-480 


-5 


-549 


-628 


-400 


-600 


-600 


mb(I )  mb{l) 


5.77 


4.98 

0.04 

5.68 

0.02 

5.37 

0.03 

5.24 

0.03 

433 

0.07 

430 

0.07 

5.67 

0.02 

.9  5.09 


53  535 


5.65 


5.1 


0.03  5.0 


3.9 


0.03  5.3 


5.3 


0.02  53 


53 

531 

0.03  5.2 

5.6 

5.56 

0.03  53 

5.8 

5.80 

0.03  5.6 

53 

530 

53 

5.46 

0.03  5.3 

53 

5.44 

0.03  5.3 

4.6 

4.64 

4.0 

5.9 

5.81 

0.03  53 

5.61 

52 

521 

0.03  5.0 
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waves  much  easier.  In  examining  the  ensemble  of  records,  it  became  apparent  that  Love 
waves  were  observed  for  almost  all  events.  The  largest  Love  waves  are  associated  with 
Pahute  Mesa  and  Rainier  events;  however,  most  Yucca  shots  also  showed  evidence  of  Love 
wave  generation.  Most  of  the  low-magnitude  events  were  from  Yucca  Flat,  so  that  the  lack  of 
Love  wave  observations  may  be  due  to  the  extremely  low  SNR.  Also,  the  smaller  events  can 
only  be  observed  at  the  nearer  stations,  so  that  the  Rayleigh  and  Love  wavefields  have  not 
separated  enough  in  time  by  dispersion  to  identify  the  Love  wave,  since  the  Love  wave  is  a 
pulse  and  arriving  near  and  within  the  Rayleigh  wave  time  window.  Figures  1  and  2  plot  the 
three  component  records  for  LON  (Longmire,  WA)  and  ALQ-ANMO  (Albequerque,  NM), 
respectively.  These  two  stations  were  the  best-reporting  for  the  network  used  in  this  study. 
Both  stations  were  upgraded  with  longer-period  instruments  for  events  after  1982,  which 
explains  the  oberved  difference  in  frequency  content  between  the  two  types  of  records.  Love 
waves  are  present  on  all  but  the  lowest  SNR  records,  implying  that  tectonic  strain  release  or 
some  other  effect  is  generating  a  non-isotropic  source.  Given  &  Mellman  (1986)  found  that 
all  events  they  studied  displayed  evidence  of  tectonic  strain  release.  However,  their  study 
was  confined  to  large  (mj,  >  5.5)  events. 

Poor  quality  data  such  as  ANMO:  85092,  82217(tang),  82127,  84152,  86318  and  87197(vert) 
were  removed.  Low  SNR  tangential  components  were  removed  in  cases  where  the  vertical 
component  was  kept,  such  as  ANMO:  87181  and  86141,  or  ALQ:  65062  and  72140.  All 
seismograms  were  bandpass-filtered  between  6  and  60  seconds  to  best  observe  the  wavetrain 
and  for  data  analysis.  Synthetic  seismograms  were  filtered  identically.  Next,  cosine-tapered 
time  windows  were  applied  to  vertical  and  tangential  records  in  an  attempt  to  remove  all  but 
the  fundamental-mode  Rayleigh  wave  and  Love  wave  respectively.  A  more  robust  alternative 


Figure  1:  Explosion-generated  surface  waves  observed  at  Albequerque,  N.M. 
(ALQ),  some  900  km  from  NTS. 
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Figure  2:  Explosion-generated  surface  waves  observed  at  Longmire,  WA  (LON), 
some  1150  km  from  NTS. 
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LON  (d) 


to  this  time-domain  windowing  would  be  to  phase-match-filter  (Herrin  &  Goforth,  1977)  the 
records,  particularly  the  tangential  component  and  all  low  SNR  records,  in  order  to  remove 
noise  and  Rayleigh  contamination,  particularly  for  nearer  stations  where  the  fundamental¬ 
mode  Rayleigh  wave  and  Love  wave  group  arrival  times  overlap. 

Next,  the  spectral  ratio  of  each  windowed  record  and  its  appropriate  synthetic  is  computed 
by  averaging  over  the  0.167  to  0.0167  Hz  bandwidth.  The  standard  deviation  in  the  spectral 
ratio  estimate  is  taken  as  the  error  in  the  observation.  Spectral  ratio  curves  for  Rayleigh 
waves  were  generally  smooth  and  flat,  implying  a  good  fit  between  observed  and  synthetic 
spectra.  Love  wave  data-to-synthetic  spectral  ratios  were  not  as  flat.  It  is  suspected  that 
the  Rayleigh  wavetrain  energy,  contaminating  the  Love  wavetrain  primarily  is  responsible 
for  this  effect. 

The  polarity  of  the  waveform,  used  as  phase  information,  was  determined  by  visual  in¬ 
spection.  This  was  primarily  of  importance  for  Love  wave  observations.  Only  in  the  case 
of  PILED  RIVER  were  there  any  reversed-polarity  Rayleigh  waves.  The  resultant  spectral 
ratios  and  phase  information  were  used  in  the  source  inversion.  The  inversion  technique 
is  described  in  the  ”  Synthetic  Seismogram  Generation  and  Moment  Inversion  Technique” 
section.  The  results  of  the  source  inversion  are  discussed  later. 

Path  Modeling 

To  determine  accurate  long-period  source  parameters  it  is  necessary  to  correct  for  propaga¬ 
tion  effects  along  the  surface  wave  travel-paths.  This  involves  determining  an  average  earth 
structure  for  each  path  from  a  general  source  region  to  a  particular  receiver.  The  paths  are 
modeled  as  1-dimensional  laterally  homogeneous  structures.  Paths  from  the  Nevada  Test 
Site  (NTS)  to  23  World  Wide  Seismic  Station  Network  (WWSSN)  and  Canadian  Seismo- 
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graph  Network  (CSN)  stations  were  determined  for  a  surface  wave  study  by  Stevens  (1986 
and  personal  communication).  These  stations  are  shown  as  circles  in  Figure  3.  The  path 
Green’s  functions  were  determined  by  a  two  step  inversion  of  shear  wave  velocity  and  Q  as 
discussed  in  Stevens  (1986).  For  other  stations  in  the  network  (shown  as  triangles  in  Figure 
3)  path  models  were  chosen  or  developed  through  a  variety  of  means.  Whenever  possible, 
pre-existing  models  were  used  for  path  structures  for  any  given  NTS-station  path.  For  some 
paths  several  models  were  tried.  To  determine  whether  or  not  a  particular  earth  model 
was  appropriate  for  a  particular  path,  a  synthetic  seismogram  was  generated  and  compared 
to  the  highest  SNR  observed  record.  The  quality  of  the  fit  was  determined  by  visual  in¬ 
spection  of  time  domain  records.  Absolute  travel  time  and  amplitude  comparisons  between 
observed  and  synthetic  intermediate-period  (6<T<60  sec.)  Rayleigh  waveforms  as  well  as 
fundamental-mode  dispersion  data  were  the  criteria  used  to  judge  the  goodness  of  fit  of  the 
paths. 

Most  of  the  intermediate-range  to  distant  station  paths  (A  >  1500  km)  were  successfully 
modeled  using  a  structure  from  Stevens’  ensemble  of  structures.  However,  for  regional  sta¬ 
tions  at  lesser  distances  (i.e.,  the  southwestern  United  States),  alternative  earth  models  were 
necessary.  Station  coverage  in  this  area  is  especially  important  for  several  reasons.  At  re¬ 
gional  distances  the  threshold  magnitude  for  observing  surface  waves  is  lower  than  that  of  a 
more  distant  receiver,  since  the  effects  of  geometrical  spreading  and  attenuation  (due  both 
to  anelastic  absorption  and  scattering)  are  less  for  shorter  paths.  Amplitude  errors  due  to 
using  incorrect  Q  values  for  path  corrections  will  also  be  less  for  shorter  paths  for  the  same 
reason.  Because  NTS  is  situated  in  the  southwest  corner  of  the  U.  S.,  the  azimuthal  coverage 
for  more  distant  paths  is  limited,  so  that  there  are  no  seismographic  stations  600  to  5000  km 


Figure  3:  Map  of  North  American  station  network  used  in  this  study.  The  "spoked 
wheel"  is  the  Nevada  Test  Site. 
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to  the  west  or  southwest  of  NTS  (i.e.,  in  the  Pacific  Ocean).  Station  coverage  to  the  south 
(for  pure  continental  paths)  is  also  very  limited,  as  there  were  few  long-period  seismographic 
stations  recording  in  Mexico  and  Central  America  during  the  time  period  in  which  events 
for  this  study  were  taken. 

Broad  azimuthal  coverage  is  also  important  for  constraining  the  seismic  moment.  For 
example,  a  vertical  strike  slip  double-couple  source,  which  has  a  (sin  2^)  radiation  pattern 
superimposed  with  an  isotropic  source  will  generate  a  Rayleigh  wave  radiation  pattern  with 
quadrants  of  increased  and  decreased  amplitude,  which  is  caused  by  the  constructive  and 
destructive  interference  of  the  two  sources.  For  the  case  of  perfect  azimuthal  coverage,  the 
network-averaged  scalar  moment  will  be  that  of  the  isotropic  source.  Lack  of  coverage  in  one 
of  these  quadrants  will  cause  the  observed  network-averaged  moment  to  be  larger  or  smaller 
than  the  actual  value,  depending  upon  whether  constructive  or  destructive  interference  oc¬ 
curs  in  the  unobserved  quadrant.  If  phase  and  amplitude  information  are  to  be  inverted  to 
obtain  the  moment  tensor,  broad  azimuthal  coverage  improves  the  constraints  of  the  moment 
tensor  inversion. 

A  large  proportion  of  the  regional  stations  in  this  study  have  paths  which  predominantly 
traverse  the  Basin  and  Range  province.  Fortunately  this  region  has  been  the  focus  of  many 
crustal  and  upper  mantle  structure  studies,  so  that  there  axe  a  good  many  earth  structures 
derived  from  both  surface  wave  and  body  wave  studies  available  from  which  to  choose  to 
model  the  paths  of  interest  in  this  study. 

Mooney  &  Braile  (1989)  review  P-wave  velocity  structures,  determined  from  P-wave  refrac¬ 
tion  studies,  for  the  different  physiographic  regions  of  North  America,  including  the  Basin 
and  Range  province.  Saikia  &  Burdick  (1991)  forward-modeled  P„/  waveforms  from  NTS 
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explosions  to  many  of  the  same  regional  stations  used  in  this  study  in  order  to  determine  a 
general  crustal  model  for  the  area. 

Press  (1960)  used  a  combined  refraction,  surface  wave  phase  velocity,  and  gravity  data  to 
infer  a  seismic  velocity  structure  for,  primarily,  the  Mojave  block.  Hadley  &  Kanamori,  in 
more  recent  studies,  have  interpreted  this  same  area  both  in  terms  of  P-wave  refraction  data 
(1977)  and  interstation  phase  velocities  from  teleseismic  Rayleigh  waves  (1979). 

Keller  et  al  (1976)  used  regional,  explosion-generated,  short-period  (4  to  22  sec.)  Rayleigh 
wave  group  velocities  to  determine  crustal  structure  for  several  areas  in  the  Western  U.  S., 
including  the  eastern  Basin  and  range  Province.  Priestly  &  Brune  (1978)  made  use  of 
teleseismic,  wide-band  (4  to  120  sec.),  interstation  fundamental-mode  Rayleigh  wave  and 
Love  wave  phase  velocities  to  determine  the  shear-wave  velocity  structure  for  the  Nevada 
and  Western  Utah  Great  Basin  area.  The  dispersion  results  of  this  study  are  corroborated 
in  a  study  by  Patton  (1982)  in  which  single-station  Rayleigh  wave  phase  velocities  were 
determined  using  regional,  digital,  broadband  data. 

Patton  &  Taylor  (1984)  determined  shear- wave  Q  structure  from  regional  Love  wave  and 
Rayleigh  wave  amplitude  data.  The  shear  velocity  model  used  in  their  method  is  a  hybrid 
of  the  Priestly-Brune  (PB,  1978)  upper  mantle  model  and  a  NTS-TUC  (Tuscon,  Arizona) 
crustal  model  developed  by  Bache  et  al.  (1978)  discussed  below.  Their  results  imply  a  much 
lower  Q  in  the  crust  and  upper  mantle  than  a  more  general  study  of  attenuation  in  western 
North  America  by  Mitchell  (1975),  which  they  attribute  to  partial  melting  in  the  upper 
mantle.  This  attenuation  model  was  used  in  conjunction  with  all  earth  models  examined  in 
this  study  for  paths  in  this  region,  which  don’t  have  an  associated  Q  structure. 

The  model  that  generally  best  fit  the  Rayleigh  wave  and  Love  wave  observed  waveforms 
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for  the  Basin  and  Range  stations  (DUG,  GSC,  ELK,  KNB,  LAC  and  MNV)  from  the  above 
cited  models  was  found  to  be  the  PB  model.  The  Patton  and  Taylor  (PT)  Basin  and  Range 
attenuation  model  was  used  in  conjunction  with  the  PB  velocity  model. 

Bache  et  al.  (1978)  determined  flat-layered  velocity  structures  for  paths  from  NTS  to 
ALQ  and  TUC  using  surface  wave  dispersion  data.  Attenuation  values  were  taken  from  the 
previously  cited  study  of  regional  Rayleigh  attenuation  study  by  Mitchell  (1975).  Langston 
&  Helmberger  (1974)  also  modeled  the  NTS-TUC  path  by  forward-modeling  body  wave 
phases  and  by  trial  and  error  inversion  of  group  velocity  data.  Both  studies  have  similar 
crustal  thicknesses  for  the  NTS-TUC  path,  as  body  wave  phases  were  used  to  constrain  the 
depth  to  the  Moho;  however,  the  shear  velocity  is  higher  at  the  top  of  the  mantle  in  the 
Langston-Helmberger  model  than  in  the  other  (vt  =  4.6,  vs.  v,  =  4.42).  The  Bache  model 
was  chosen  for  this  study.  Synthetic  fundamental-mode  surface  waves  generated  with  the 
Bache  models  match  observed  surface  wave  records  significantly  better  than  the  Langston- 
Helmberger  model.  Comparisons  of  observed  and  syntethic  surface  wave  waveforms  for  the 
NTS  to  ALQ  and  TUC  paths,  along  with  all  others,  are  shown  in  Figure  4  (for  Rayleigh 
waves)  and  Figure  5  (for  Love  waves). 

Path  structures  determined  in  this  study  were  arrived  at  through  a  combination  of  forward 
modeling  and  inversion  of  dispersion  data  for  shear-wave  structure.  Compressional  velocity 
and  density  were  calculated  from  functional  relationships  with  shear-wave  structure  as  per 
Stevens  (1986).  The  philosophy  behind  these  path  structures  is  to  obtain  earth  models  that 
will  give  accurate  path  corrections  or  Green’s  functions  with  which  to  generate  synthetic 
seismograms  which  fit  the  observed  waveforms  with  respect  to  amplitude  and  arrival  times, 
rather  than  to  recover  the  lithospheric  structure  or  the  earth  along  the  path,  although  the 


Data  vs.  Synthetic  (vertical  component) 


Figure  Comparison  of  vertical  component  fundamental  Rayleigh  wave  wave¬ 
forms.  The  data  time  series  are  the  upper,  thicker  traces;  the  lower  trace  in  each 
case  is  the  fundamental  mode  synthetic.  All  time  series  band-passed  between  60  and 
6  seconds. 


Data  vs.  Double-Couple  Synthetic  (tangential  component) 


Comparison  of  observed  and  synthetic  fundamental-mode  Love 
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models  also  should  not  diverge  too  drastically  from  realistic  gross  earth  structure.  NTS  sur¬ 
face  waves  for  many  of  the  stations  in  this  study’s  network  traverse  several  tectonic/geologic 
provinces.  A  flat-layered  earth  model  obviously  cannot  accurately  describe  a  crust  and  up¬ 
per  mantle  cross-section  traversing  various  portions  of  the  Western  Cordillera  (the  Basin  and 
Range,  Rocky  Mountains,  and  the  Colorado  Plateau)  and  the  North  American  craton.  How¬ 
ever,  one  can  obtain  a  model  which  describes  an  average  earth  model  for  a  given  path  which 
fits  the  observed  dispersion  data  and  spectral  amplitude  curves,  as  well  as  the  time-domain 
waveforms. 

The  starting  models  used  for  this  path  determination  scheme  come  from  studies  in  which 
the  earth  structures  for  similar/nearby  paths  were  inverted  for  using  rigorous  inversion  pro¬ 
cedures  and  resolution  analysis  and  geophysical  constraints.  The  first  step  was  to  see  how 
well  the  initial  structure  modeled  Rayleigh  waveforms.  As  it  turned  out,  none  of  these  struc¬ 
tures  fit  both  arrival  times  and  waveforms  very  well  for  the  desired  paths.  Next,  for  each 
path  the  initial  structure  was  perturbed  in  an  effort  to  better  fit  synthetic  seismograms  gen¬ 
erated  with  the  structure  to  the  observed  waveforms.  This  included  adding  and  removing 
low- velocity  surface  layers  (representing  sedimentary  basins),  increasing  and  decreasing  the 
gross  (or  average)  velocity  in  the  crust  and/or  upper  mantle,  and  increasing  and  decreas¬ 
ing  the  crustal  thickness.  This  exercise  also  gives  one  a  feel  for  the  effects  of  these  model 
changes.  Perturbations  to  the  initial  model  were  retained  if  they  improved  the  wave-train 
arrival  and/or  waveform. 

The  depth  to  the  moho  and  the  shear-velocity  contrast  there  are  the  most  important 
influences  upon  the  surface  wave-field.  Independent  geophysical  constraints  on  these  factors 
was  desired.  Mooney  &  Braile  (1989)  summarize  crustal  models  for  North  America  (N A) 
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inferred  from  seismic  refraction  and  reflection  studies.  Their  map  of  NA  crustal  thickness 
and  representative  P-wave  velocity  structures  for  the  various  geologic  regions  were  useful  as 
guidelines  for  estimating  average  depths  to  the  moho  for  the  various  paths  investigated,  as 
well  as  for  inferring  the  crust-mantle  shear- wave  velocity  contrast.  In  inverting  for  shear 
velocity  structure,  as  described  below,  it  was  found  that  using  a  starting  model  with  a 
correct  crustal  thickness  (constrained  by  reflection  and  refraction  data)  was  important  in 
order  to  obtain  earth  structures  which  produced  well-fitting  dispersion  values  and  time-series. 
Inverted  shear-wave  earth  models  obtained  from  starting  models  with  crustal  depths  which 
differed  from  what  was  believed  to  be  a  reasonable  estimate  of  the  average  crustal  thickness 
(constrained  by  independent  geophysical  data)  usually  did  not  produce  dispersion  values 
compatible  with  the  data  or  time  series  which  adequately  modeled  observed  seismograms. 

As  an  independent  test  of  crustal  thickness,  a  stochastic  inversion  of  the  dispersion  data 
was  conducted  with  a  generic  model  ( vp  =  7.0, v,  =  4.0)  to  find  the  depth  to  the  moho.  It 
was  found  in  each  case  that  the  structure  converged,  to  first  order,  to  a  layer  over  a  half¬ 
space.  The  depth  of  the  basal  crustal  layer  was  found  to  agree  with  the  values  inferred  from 
the  Mooney  &  Braile  (1989)  map.  In  the  case  of  JAS,  where  the  travel  path  transverses 
the  Sierra  Nevada,  it  was  unclear  what  the  average  crustal  thickness  should  be,  so  the 
preliminary  inversion  inferred  mantle  depth  of  37.5  km  was  taken  to  be  representative  of  the 
average  NTS-JAS  path,  and  the  Priestly-Brune  Basin  and  Range  based  starting  model  was 
modified  accordingly. 

Paths  from  NTS  to  PAS,  JAS,  and  RSTN  stations  were  modeled  in  this  study.  Elastic 
parameters  for  these  paths  were  determined  by  inverting  fundamental-mode  Rayleigh  wave 
dispersion  data  to  obtain  the  shear  wave  velocity  (/?)  structure  for  the  model.  The  sur- 
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face  wave  analysis  and  velocity  inversion  computer  codes  used  come  from  the  Computer 
Programs  in  Seismology  software  package  developed  at  the  University  of  St.  Louis  (Her¬ 
rmann,  1988).  The  shear-wave  velocity  (Vs)  structure  is  found  by  simultaneously  inverting 
for  group  and  phase  velocity  data  (for  Rayleigh  waves  measured  between  0.015  and  0.150 
Hz)  and  minimizing  the  integral  of  ||d/?/<iz||  over  the  structure,  except  across  discontinuities 
at  layer  boundaries.  This  linear  inversion  scheme  is  similar  to  the  one  described  in  Bache  et 
al.  (1978)  and  Stevens  (1986). 

Q  values  for  paths  to  RSTN  stations  were  taken  from  the  Stevens’  attenuation  model 
with  the  most  similar  path.  All  RSTN  stations  were  on  or  near  great  circle  paths  of  North 
American  WWSSN  or  CSN  stations.  The  PAS  and  JAS  structures  used  the  Patton  &  Taylor 
(1984)  attenuation  model.  These  two  stations  have  paths  to  NTS  which  traverse  a  portion 
of  the  Basin  and  Range  region. 

Although  Q  structure  was  not  determined  by  inversion,  it  was  modified  from  the  starting 
Q  model  by  forward  modeling  of  synthetic  time  series  to  be  compatible  with  observations. 
For  some  paths  it  was  necessary  to  create  narrow,  very  high  Q  (Q  «  1000)  zones  in  the 
upper  crust  (5  to  15  km  depth  )  in  order  to  model  the  data.  This  was  needed  when  the 
observed  Airy  Phase  amplitudes  (8  <  T  <  15)  were  considerably  larger  relative  to  the 
longer  period  data  than  the  synthetic  counterpart.  The  actual  cause  for  these  amplitude 
discrepancies  are  likely  to  be  lateral  variations  in  the  waveguide  which  tend  to  generate 
shorter  period  surface  waves,  rather  than  extreme  Q  values  in  the  crust  or  mantle.  However, 
to  model  these  paths  with  flat-layered  models,  amplitude  discrepancies  must  be  accounted 
for  in  the  attenuation  structure.  The  amplitude  discrepancy  between  relative  long  period 
and  shorter  period  Rayleigh  waves  is  in  all  likelihood  due  to  anomalous  short-period  energy, 
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as  it  is  much  harder  to  account  for  differences  in  Q  at  longer  periods,  since  the  long-period 
wave-field  completes  fewer  cycles  for  a  given  path  distance  than  do  the  shorter  period  waves. 

Surface  wave  dispersion  is  not  strongly  dependent  on  the  compressional  velocity  ( VP )  or 
density  (p),  hence  they  aren’t  directly  solved  for  in  the  inversion,  but  rather  treated  as 
functions  of  the  medium’s  shear  velocity,  to  which  surface  wave  velocities  are  most  sensitive. 
The  compressional  velocity  (V»  is  constrained  to  be  consistent  with  a  Poisson’s  ratio  of 
0.27,  while  the  density  is  constrained  by  a  Birch’s  law  formula: 

p  =  0.65  x/3  +  400. 

This  formula  is  an  empirical  relationship  developed  from  data  in  Dobrin  (1976)  and  is  in 
MKS  units. 

The  maximum  depth  of  the  inversion  model  was  set  to  the  approximate  longest  wavelength 
(A  =  c/f)  observations.  For  a  0.015  Hz  Rayleigh  wave  phase  velocity  (c  «  4.0),  this  depth 
would  be  270  km.  Little  resolution  is  to  be  expected  at  the  lower  bounds  of  this  depth 
estimate,  however.  Layer  thicknesses  in  the  crust  were  chosen  to  be  5  km,  with  the  top-most 
layer  being  divided  in  two,  allowing  for  some  resolution  of  possible  shallow  sedimentary 
basins.  The  upper  40  km  of  mantle  is  divided  into  10  km  thick  sections  and  below  this  to 
200  km,  it  is  divided  into  layers  20  km  thick.  These  layer  thicknesses  are  close  in  value  to 
the  resolution  length. 

Single-station  group  velocity  (£7)  dispersion  values  were  determined  from  multiple  filter 
analysis  (Dziewonski  et  al,  1969)  of  vertical  component  Rayleigh  waves.  After  obtaining 
group  velocities  (£7(u>)),  the  time  series  are  iteratively  phase  matched  filtered  (Herrin  & 
Goforth,  1977;  and  Stevens,  1986)  to  reduce  scattered  energy,  noise  and  higher  modes,  so 
that  the  signal  can  be  more  accurately  analyzed  for  phase  and  group  velocity  and  amplitude 
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information.  The  phase  information  obtained  from  the  phased  matched  filter  is  used  to 
obtain  a  phase  velocity.  Two  iterations  were  made  to  obtain  the  final  dispersion  data  used 
to  invert  for  a  shear  velocity  structure. 

Seismograms  used  for  dispersion  analysis  were  chosen  from  large  magnitude  explosions 
with  the  least  observed  tectonic  release.  Events  which  did  not  show  evidence  of  long-period 
tangential  energy  in  the  approximate  Love  wave  time  window  were  assumed  to  be  good 
candidates.  KEARSARGE  data  were  used  for  PAS  and  JAS  (CMB)  to  NTS  paths.  Dis¬ 
persion  analysis  was  done  for  other  events,  too.  For  all  events  examined,  coherent  spectral 
amplitudes  were  observed  up  to  30  sec.  At  periods  greater  than  this,  the  error  in  velocity 
measurements  are  considerably  larger.  The  propagation  paths  axe  less  than  400  km  in  length, 
so  longer-period  surface  waves  have  traveled  very  few  wavelengths  and  thus  axe  not  yet  fully 
developed  and/or  are  not  dispersed  enough  to  measure  their  group  energy  packets. 

For  the  RSTN  stations,  seismic  waveforms  came  from  a  variety  of  events.  For  these  stations 
at  greater  distances  (D  <  230  km),  coherent  signal  was  observed  up  to  50  sec.  The  RSSD, 
RSCP,  and  RSNT  paths  used  DARWIN  data.  The  DARWIN  group  velocity  curves  for 
RSON  and  RSNY  had  several  inflection  points  between  15  and  50  seconds.  This  feature  is 
believed  to  be  due  to  signal  noise  or  processing  artifacts.  These  two  stations’  deconvolved 
displacement  records  showed  long-period  noise  throughout  the  records.  For  other  larger 
events,  with  higher  signal  to  noise  ratios,  this  phenomenon  may  be  due  to  tectonic-release 
generated  signals.  These  deviations  from  the  correct  dispersion  curve  are  significant  and 
adversely  affect  inversion  results. 

The  group  velocity  dispersion  for  a  suite  of  seven  events  (both  from  Pahute  Mesa  and  Yucca 
Flat)  for  the  two  paths  were  compared.  The  group  velocity  curves  generally  correlated  quite 
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well  by  visual  inspection  between  8  and  40  sec.  Above  and  below  these  levels,  only  some 
events  had  coherent  smoothly  continuous  dispersion  curves.  In  the  period  range  between  15 
and  50  sec.  perturbations  from  a  smooth,  near-constant  slope  group  velocity  curve  varied 
in  shape  from  event  to  event  and  sometimes  were  absent,  implying  that  they  are  some  sort 
of  noise  effect  and  not  a  propagation  phenomenon.  The  final  dispersion  data  came  from 
BELMONT  for  RSON  and  HERMOSA  for  RSNY.  Deconvolved  displacement  seismograms 
for  these  events  show  low  long-period  noise  levels. 

The  best-fitting  forward-modeled  velocity  structure  was  then  used  as  the  initial  model 
for  the  shear-wave  inversion  procedure.  A  weighted,  least-squares  inversion  of  fundamental¬ 
mode  Rayleigh  wave  group  and  phase  velocity  is  performed.  Details  of  the  inversion  method 
are  given  in  Russell  (1986).  The  procedure  is  non-linear  in  nature,  so  iterative  solutions 
are  obtained.  Q  structure  was  not  inverted  for  these  paths.  Instead  the  attenuation  model 
associated  with  the  initial  inversion  model  was  used.  Here  the  assumption  is  that  the  gross 
attenuation  structure  is  the  same  between  nearby  paths.  To  robustly  solve  for  attenuation, 
one  must  employ  two-station  attenuation  measurements  or  joint  moment- Q  inversions.  The 
combination  of  data  and  paths  in  this  study  does  not  lend  it  to  the  two-station  method.  For 
the  shorter  paths  (JAS  and  PAS),  attenuation  effects  are  minimal  for  surface  waves,  as  they 
propagate  relatively  few  wavelengths  and  consequently  undergo  a  corresponding  low  number 
of  attenuation  cycles. 

The  shear  wave  inversion  is  a  non-linear  process,  so  many  iterations  were  performed.  The 
weighting  factors  used  for  the  dispersion  observations  were  proportional  to  the  inverse  of  the 
spectral  amplitude  at  a  given  period.  Approximately  every  fourth  iterated  model  was  saved 
in  case  later  iteration  models  converged  towards  a  ”  pathological”  one  with  unrealistically 


large  velocity-contrast  low  velocity  zones  in  the  crust.  Early  inversion  attempts  made,  using 
a  damped  least-squares  stochastic  inversion,  resulted  in  such  low  velocity  zones.  Applying  a 
weighted,  damped  differential  least-squares  inversion  produced  "smoother”  velocity  models 
which  were  geophysically  reasonable  ( i.e .,  did  not  have  oscillating  high  and  low  velocity  zones 
or  unreasonably  large  velocity  contrasts).  The  final  inverted  models  were  obtained  using  the 
latter  method. 

As  mentioned  earlier,  the  criterion  to  determine  the  goodness  of  fit  of  an  earth  structure 
were  visual  comparisons  of  the  observed  and  Green’s  Function  derived  dispersion  and  wave¬ 
form  fits  of  the  fundamental  Rayleigh  wave.  Figures  6  and  7  compare  observed  and  modeled 
group  and  phase  velocity  dispersion  for  the  seven  paths.  All  dispersion  curve  fits  are  quite 
good.  The  largest  mismatch  between  observed  and  inverted  model  dispersion  values  are  at 
the  longest  and  shortest  periods.  These  periods  had  the  lowest  spectral  amplitudes,  so  they 
are  weighted  least  in  the  inversion  scheme.  It  is  somewhat  surprising  that  such  good  fits 
were  obtained,  considering  the  significant  lateral  variations  in  the  waveguide,  particularly 
for  the  JAS  and  PAS  paths  which  each  traverse  several  mountain  ranges  and  valleys. 

Figure  4  compares  observed  to  synthetic  seismograms  for  the  fundamental-mode,  vertical 
component  Rayleigh  waves  for  the  seven  paths  modeled  in  this  study  as  well  as  for  the  other 
network  paths  used  in  this  study.  All  time  series  were  band-passed  filtered  between  60  and 
6  seconds.  Each  path  is  well  modeled  with  respect  to  waveform  and  arrival  time  (At  <  1.0). 
A  confirmation  of  how  well  the  models  reflect  the  average  earth  structure  of  the  path  are 
comparisons  of  observed  and  synthetic  fundamental-mode  Love  waves.  The  error  in  timing 
of  the  Rayleigh  wave-train  is  in  all  cases  small  and  is  estimated  to  be  less  than  2  seconds. 
Figure  5  compares  observed  and  synthetic  Love  waves.  Unlike  the  case  of  explosion  generated 
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Figure  6:  Fundamental-mode  Rayleigh  wave  dispersion  curves  for 
JAS,  PAS,  RSSD  and  RSON. 
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Figure  7:  Fundamental-mode  Rayleigh  wave  dispersion  curves  for 
RSCP,  RSNT  and  RSNY. 


Rayleigh  waves,  the  polarity  of  tangential  waves  is  ambiguous,  because  the  source  mechanism 
is  not  necessarily  known.  The  correct  polarity  or  phase  was  determined  by  comparing  the 
best  combined  waveform  fit/travel  time  residual  minimum  in  conjunction  with  assumptions 
on  the  radiation  pattern  from  previous  studies.  The  waveform  fits  are  sufficient  for  inversion 
purposes  and  aren’t  bad  predictions  of  the  observed  Love  waves,  although  the  wave-train 
arrival  time  errors  are  larger  (from  2  to  10  seconds),  with  the  largest  errors  correlating,  in 
general,  to  the  more  distant  paths.  Tables  2  lists  the  parameters  for  the  final  inverted  path 
models. 

Figure  8  displays  the  data/synthetic  seismograms  for  the  seven  paths  modeled  in  this  study. 
Both  fundamental-mode  Rayleigh  waves  (vertical  component)  and  Love  waves  are  plotted. 
As  discussed  previously,  the  Q  structures  are  based  on  pre-existing  Q  models,  which  are  then 
modified  by  forward  modeling  the  waveforms  and  only  (3  was  directly  inverted.  Estimated 
errors  (1  standard  deviation)  in  S  velocity  are  on  the  order  of  0.05  to  0.1  km/sec.  Errors  are 
smallest  in  the  upper  crust  and  largest  in  the  lower  crust  and  upper  mantle.  The  resolving 
kernel  widths  were  approximately  5  km  in  the  upper  crust  and  10  to  15  km  in  the  lower 
crust  and  upper  mantle.  Little  or  no  resolution  was  found  below  60  km.  The  P-wave  and 
S-wave  velocities,  Q  values  and  densities  for  these  structures  can  be  found  in  Woods  (1993). 
These  path  models  have  typical  continental  path  structures  with  crustal  thicknesses  in  line 
with  refraction  and  reflection  constraints.  No  large  low  velocity  zones  are  within  the  crust. 
Where  an  overlying  layer  does  have  a  greater  S  velocity,  the  contrast  is  small  (A Vs  <  0.2 
km/sec.),  and  the  errors  in  Vs  are  nearly  of  the  same  order  of  magnitude. 
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NTS  to  RSON  Earth  model 
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Table  2a:  NTS  path  models  for  JAS,  PAS,  RSSD  and  RSON 
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Synthetic  Seismogram  Generation  and  Moment  Inversion  Technique 

Since  the  explosions  are  shallow  sources  (ft  <  1  km)  and  tectonic  strain  release  is  also 
thought  to  be  not  much  deeper,  higher  modes  are  not  strongly  excited  as  the  fundamental¬ 
mode- Rayleigh  wave  at  the  periods  of  interest.  Patton  (1988)  found  that  spall  is  an  effectice 
generater  of  higher  modes.  However,  such  excited  energy  is  relatively  short-period,  so  that 
it’s  effect  on  the  long-period  function  is  minimal. 

In  the  following  we  give  the  expressions  for  the  Rayleigh  and  Love  surface  waves  due 
to  an  explosion  and  a  double  couple  (Ben-Menahem  &  Harkrider,  1964,  Harkrider,  1970, 
Mendiguren,  1977,  among  many  others),  These  epressions  have  been  modified  for  a  slowly 
varying  laterally  inhomogeneous  medium  as  in  Woods  &  Harkrider  (1995),  where  the  source 
region  is  one  verticaly  inhomogeneous  medium  and  the  propagation  path  and  receiver  are 
another. 

For  the  vertical  spectral  displacement  (positive  down),  the  Rayleigh  waves  yield 
w0  =  M#  |sin  A  sin  2 6  W45  —  ^cos  A  cos  6  cos  $  —  sin  A  cos  26  sin  <frj  W13 

+  —  sin  26  cos  2<})  +  cos  A  sin  6  sin  2^  Wi2^  +  M*  W*  (1) 

The  angles  A,  6,  and  <f)  are  the  standard  rake,  dip  and  azimuth  measured  from  the  fault 
strike.  The  double  couple  spectral  components  are 

W12  =  [j4#]x  ARm  H?]  f—r)  exp[— 72  r] 

Z  Cl  \C2  / 

W13  =  [C#]i  ARm  H[2)  f—r)  exp[— 72  r]  (2) 

2  Cl  \C2  / 

w45  =  -J—  [£#]i  ARm  H{02)  (~r)  exp[-72  r) 

4  Cl  VC2  / 
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and  the  explosion 


W*  =  —i—  ^  [Kr]i  ARm  exp[— 72  r] 

Cl  CL\  \C2  / 


where 


(3) 


Kr  =  y*W  ~  2^kMh)' 


A. #  =  -y3(h), 


B*  -  ~  1  1 3  *■  43  I  vM  + 


cr 


2  yz(fe)  1 

pa2A:  J  * 


or  in  terms  of  the  poisson  ratio,  <7, 


t?  f(1  +  <T)i/ m  i  2ya(&)] 


C'# 


y4(h) 


D 

k  =  - 
c 


and  the  are  the  cylindrical  Hankel  functions  representing  outward  propagation  over  a 
distance  r.  The  above  spectral  mixed  path  excitation  is 


ARm  —  [  Am  Ar2  ]^2 


(4) 


with  each  subscripted 

=  _ L_ 

~  2cUl 

where  the  energy  integral  is 

i-  f  piz)[yliz)  +  yl(z)]dz, 

Jo 

p(z)  is  the  local  density  distribution  in  the  medium  and  we  have  used  Saito’s  (1967)  Rayleigh 
wave  eigenfunction  notation,  y;(z).  The  eigenfunctions  are  normalized  in  such  away  that  the 
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vertical  displacement  eigenfunction,  y^z)  is  equal  to  1  at  the  free  surface,  z  —  0.  This 
results  in  the  horizontal  displacement  eigenfunction,  y3(z)  ,  being  equal  to  the  Rayleigh 
wave  surface  ellipticity  at  this  boundary.  y3{z)  and  j/4 (z)  are  the  normal  and  shear  stress 
eigenfuntions  associated  Rayleigh  modes.  U  and  c  are  respectively  the  local  group  and  phase 
velocities.  By  local  we  mean  the  eigenvalues  and  eigenfunctions  that  one  would  obtain  for 
a  laterally  homogeneous  half-space  consisting  of  the  vertical  elastic  and  density  distribution 
at  that  location.  The  subscripted  quantities  other  than  the  just  mentioned  eigenfunctions 
and  Hankel  functions  are  as  follows:  the  subscript  2  denotes  local  quantities  at  the  receiver 
location,  which  in  this  case  is  the  same  as  the  propagation  path,  and  1  is  the  point  source 
location  and  quantities  within  the  1  or  2  subscripted  square  brackets  are  evaluated  at  these 
locations.  7  is  the  frequency  dependent  attenuation  coefficient  due  to  the  anelastic  structure 
of  the  path,  ie.  7  =  u/(2QU),  where  Q  is  the  attenuation  quality  factor. 

For  the  azimuthal  surface  spectral  displacement  (positive  for  increasing  azimuth),  the  Love 
waves  yield 

v0  =  M#  |(2 cos  A  sin  8  cos  2<p  -  sin  A  sin  2 8  sin  2<f>)V\2 

+  (sin  A  cos  28  cos  ^  —  cos  A  cos  8  sin  <f>)Vi 3|  (5) 

where  the  double  couple  spectral  components  are 

^12  =  — 2 ~  lyi(h#)h  Alw.  exp[ — y2  r] 

V,3  =  ~  [G#].  Aun  (|r)  exp[— 7,  r]  (6) 

°*  = 


where 
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and  the  eigenfunctions  and  eigenvalues  are  those  for  the  fundamental  Love  wave  mode  and 
Aim  is  the  mixed  path  excitation  function  for  Love  waves  defined  as  its  Rayleigh  wave 
equivalent  above  with  the  exception  that 

I  =  [  p(z)[yl(z)]dz. 

Jo 

where  yi{z)  is  the  surface  azimuthal  displacement  eigenfunction  and  y2(z)  is  the  Love  wave 
shear  stress  eigenfunction,  both  normalized  to  surface  displacement.  Whenever  the  distinc¬ 
tion  between  Rayleigh  and  Love  wave  quantities  might  be  confusing,  we  will  denote  them 
with  R  or  L  subscripts  or  superscripts.  The  choice  between  subscript  or  superscript  is  made 
to  reduce  crowding. 

Assuming  that  at  the  periods  of  interest  {T  >  6s ec)  the  explosion  and  tectonic  release 
source  occur  simultaneously  and  their  respective  moments  can  be  considered  as  step  function 
moments,  z.e.M  =  M/(iu )  ,  the  spectral  far-field  vertical  components  of  the  Rayleigh  wave 
displacements  and  the  tangential  components  of  the  Love  waves  are  respectively 


/  2  \  */2  p~ 

w0  =  -ikRARm  exp[-ikR2r  +  t7r/4]  M*—KR(h *) 

+M*.  ji  sin  A  sin  2 S[B#(h#)  —  A#(h*)  cos  2 <f>] 


/32 


—  cos  A  sin  8  sin  2 <f>A#(h#) 


— i( cos  A  cos  8cos<f>-  sin  A  cos  28 sin  <i>)C#{h#)  j 


exp  [-7/12  r] 


and 


/  2  n1/2  M # 

Vo  =  kLALm  exp [-ikL2r  +  27r/4]— 


(7) 
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and 

/  2  \1/2  M& 

v0  =  kLALm  exp [~ikL2r  +  ix/4]— 

•(2  cos  A  sin  £  cos  2<p  —  sin  A  sin  28  sin  2 <p)  exp[7£2  r].  (10) 

At  NTS  the  assumed  mechanism  for  tectonic  release  is  believed  to  be  predominantly  ver¬ 
tical,  right-lateral  strike  slip  in  nature  (Toksoz  et  al.,  1965;d  Wallace  et  al.,  1983,  1985) 
(A  =  180°  and  8  =  90°).  This  assumption  is  made  in  part  to  help  simplify  the  inversion 
procedure  and  leads  to: 

/  o  \  1/2  „  Q2  f  M*  a2 

w0  =  -ikRARm  exp[-ifcfl2r  +  tx/4]yf (0)— ■  Mj  -  --^sm2<f>  (11) 

and 

/  2  \ 1/2  Mjt 

v0  =  -kLALm  :)  exp [~ikR2r  +  iir/4]  —  cos 2 <p,  (12) 

so  that  measured  moment  becomes  a  function  of  Mi,  M#  and  (p  (three  parameters),  where 
<f>  is  the  azimuth  measured  clockwise  from  the  strike  of  the  fault  plane.  In  the  above  we  have 
replaced  the  actual  explosion  moment  M *  with  it’s  estimate  Mi  based  on  these  assumptions. 

To  invert  for  the  combined  isotropic/double  couple  source  mechanism,  the  measured  funda¬ 
mental  Rayleigh  and  Love  wave  amplitudes  are  averaged  over  a  spectral  bandwidth  between 
10  and  60  seconds.  For  a  receiver  at  azimuth  6,  measured  clockwise  from  North,  and  with 
a  vertical  strike  slip  fault  with  strike  ip,  also  measured  clockwise  from  North,  the  observed 
Rayleigh  amplitude  will  be: 

Allay  OC  Mj  -  sin  2(6  -  ip), 

which  can  be  simplified  to: 

ARay  cm  Mi  -  ^  ^  [sin  26  cos  2 ip  -  cos  26  sin  2 ip] .  (13) 

Z  p 
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This  equation  can  be  formulated  as  a  linear  least-squares  inversion  problem  of  the  form: 

A  •  x  =  b, 
or 

.  sin  26  a2  cos  26  a2 

A^,=  +  -13  (14) 

where  X\  —  Mi ,  x2  —  M#  •  cos  2 ip,  and  S3  =  M#  •  sin  2 ip. 

Similarly  for  the  Love  wave  case,  the  amplitude  is 

Alov  oc  M#  cos  2<p  =  M #  cos  2(6  —  ip), 
or 


Aim  oc  A/#[cos  26  cos  2 ip  +  sin  26sin2ip)  (15) 

or  in  terms  of  a  least-square  problem  in  the  form 

Alo v  oc  cos  28  ■  x[+  sin  26  •  x'2  =  cos  26  ■  x2  +  sin  26  •  s3.  (16) 

The  errors  in  this  inversion  can  be  calculated  directly.  For  a  function  x  =  f(u,v),  the 
variance  (cr*)  is  defined  as 

2  o / \2  2 . >,  .dx..dx. 

’•  =  'W  +  *  3^  +  <17) 

For  the  case  of  u  and  v  being  uncorrelated  a^v  =  0.  ax  is  defined  as: 

with  N  being  the  number  of  observations  and  m  the  number  of  parameters 
For  the  Rayleigh  wave  inversion  (eq.  (14))  case,  the  variance  in  the  explosion  estimate,  xx 
is  crjj  =  a\.  For  the  tectonic  release  component  (TR  =  Jx\  -I-  x|), 


(19) 


29 


For  the  strike  azimuth  (V’  =  ^jr  arctan(x3/x2)), 


2  _  2  /■  90°  >2 

aV>  -  a2  •  (  — ) 


7T  '  (x^  +  x|) 


2  ,90°., 

+  ^-(—  f 


X  '  (x|  +  x|)‘ 


(20) 


Since  in  the  Love  wave  case  TJ2  =  yxf+xf  as  well,  the  variances  GjR  and  are  also 
given  by  eqs.(19)  and  (20),  respectively. 


Error  Analysis  of  Assumptions 

Since  reducing  the  error  in  explosion  moment  estimates  is  critical  to  their  worth,  it  is 
important  to  appreciate  the  margin  of  error  introduced  by  assumptions  in  the  analysis. 
One  such  question  is  how  reasonable  is  the  vertical  strike-slip  mechanism  assumption  for 
determining  the  tectonic  release  moment  and  what  are  the  associated  errors  in  the  explosion 
moment  estimate. 

A  theoretical  source  inversion  study  was  conducted  in  order  to  get  a  handle  on  such  errors. 
A  hypothetical  network  with  the  same  station  azimuth  make-up  as  that  of  the  observed  data 
set  was  used;  however,  all  paths  used  the  same  structure  and  were  of  the  same  path  length 
(A  =  40°).  An  explosion  of  given  moment  (1016iV — m)  is  superimposed  with  a  double-couple 
source  of  varying  focal  mechanism  and  F-factor. 

Time-domain  and  spectral  moments  and  isotropic  moments  and  their  standard  deviations 
are  calculated  from  the  synthetic  ’’data.”  The  results  are  shown  in  Table  3.  At  the  top 
of  each  box  the  rake  (A)  and  dip  (5)  are  given.  A  negative  A  denotes  a  normal  fault,  a 
positive  A,  a  thrust  fault.  The  left-most  column  gives  the  ratio  of  double-couple  to  isotropic 
moment.  The  first  column  gives  the  ratio  of  input  tectonic  release  moment  (M#)  to  explosion 
moment  (Mi).  For  comparison,  log-moments  are  converted  to  surface  wave  magnitudes  by 
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the  relation 

Ms  =  log  (Mo)  -  11.36,  (21) 

developed  in  Woods  &  Harkrider  (1995)  for  the  generic  path,  so  that  the  correct  Ms  (or 
converted  Mo)  of  the  explosion  alone  should  be  4.64  on  this  scale.  Only  Rayleigh  wave 
synthetic  data  were  used  to  determine  the  source  information.  Otherwise  moments  were 
determined  using  the  same  methods  used  on  the  data  here  and  in  Woods  Harkrider 
(1995).  The  M4’s  in  the  second  column  were  obtained  by  time  domain  measurements  of  the 
maximum  peak  to  trough  amplitude  using  the  von  Seggern  (1977)  formula  for  well  dispersed 
Rayleigh  waves  as  in  Woods  &:  Harkrider  (1995). 

Ms  =  log (A/T)  +  1.08  x  log(A)  +  4.38,  (22) 

Although  the  formula  was  chosen  because  the  coefficient  (1.08)  more  closely  approximated 
the  effect  of  attenuation  along  continental  paths,  it  should  be  noted  that  Herak  &  Herak 
(1993)  obtained  a  similar  result  for  a  set  of  250  earthquakes  mostly  in  the  distance  range  of 
20  to  150°.  Their  result 

Ms  =  log  (A/T)  +  1.094  x  log(A)  +  4.429,  (23) 

agrees  with  the  standard  Prague  formula  at  100°,  where  the  data  concentration  was  high¬ 
est  and  reduces  the  significant  distance  dependence  found  in  the  5514  Prague  formula  Ma 
readings  that  composed  their  data  set. 

The  F  factor  results  in  Table  5  of  the  next  section  imply  that  for  NTS,  the  M#/M/  ratio  is 
never  above  0.6,  so  the  maximum  moment  ratio  of  2.0  was  considered  an  extreme  value  for 
this  area.  The  rake  varied  between  -45  and  45  in  15  degree  increments  and  the  strike  varied 
between  90  (vertical)  and  45  degrees  in  15  degree  increments.  These  fault  parameters  are 
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reasonable  variations  from  the  assumed  vertical  (6  =  0),  pure  strike-slip  (A)  mechanism  for 
NTS  tectonic  release. 

Unlike  the  observational  results,  there  is  a  significant  reduction  in  variance  of  the  explosion 
moment  estimate  by  inverting  for  the  isotropic  moment  component  over  the  average  spectral 
moment  estimate.  For  all  but  the  most  extreme  deviations  in  source  parameters  and  F  factor, 
the  estimated  isotropic  moments  standard  deviations  are  small,  usually  under  0.01  magnitude 
units.  The  actual  data  errors  (Table  6)  are  somewhat  bigger.  This  can  be  attributed  to  the 
fact  that  in  the  synthetic  study  case  the  seismograms  have  no  noise  contamination,  hence 
there  is  no  problem  measuring  the  signal  accurately.  Also  in  the  synthetic  study  case,  far 
better  azimuthal  coverage  was  attainable,  since  only  a  small  proportion  of  the  network  was 
available  for  any  one  event.  In  the  theoretical  case  the  synthetic  seismograms  are  identical 
to  the  synthetic  ”data,”  so  that  there  will  be  no  error  in  the  spectral  ratio  measurements, 
whereas  in  the  observational  case,  errors  in  modeling  the  propagation  path  can  lead  to 
inaccuracies  in  moment  estimates. 

The  average  spectral  scalar  moment  (Af0)  standard  deviations  are  significantly  larger  than 
those  of  the  inverted  moment  estimates,  being  similar  in  size  to  errors  in  time-domain  moment 
estimates.  These  errors  are  significant  even  for  relatively  low  F-factor  events  < 

0.75),  being,  on  average,  0.07  for  M#/Mj  =  0.25  events  and  0.12  for  M#/Mj  =  0.5  events. 
For  large  F-factor  events  (M#/Mj  >  1),  the  variances  become  exceedingly  high,  implying 
large  scatter  in  the  moment  estimate. 

The  inverted  isotropic  moment  is,  in  most  all  cases,  closer  in  estimating  the  explosion 
moment  than  either  average  moment  estimate.  For  all  but  the  most  extreme  cases  (say, 
X  _  _45°  and  g  =  60°  and  M#/Mj  >  0.75),  the  difference  between  the  estimated  and 
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actual  isotropic  moment  is  less  than  0.1  magnitude  units  and  is  often  less  than  0.01  log 
units,  whereas  for  the  M0(ptpk)  and  M0(w),  this  difference  is  larger  sometimes  as  great  as 
0.24  magnitude  units.  For  a  pure  vertical  strike-slip  tectonic  release  mechanism,  the  average 
estimated  moment  or  Ms  determined  from  a  network  with  complete  azimuthal  coverage 
should  be  that  of  the  explosion,  with  the  variations  due  to  the  tectonic  release  canceling 
itself  out,  on  average  (Helle  t  Rygg,  1984).  However,  the  network  used  does  not  have 
equal  coverage  for  each  lobe  of  the  sin  2 6  radiation  pattern,  so  that  the  apparent  moment  or 
magnitude  differs  from  the  actual  value  for  even  the  smaller  F  factor  cases.  For  the  case  of 
M#/Mi  =  0.75,  this  difference  is  0.17  (48  percent).  These  canonical  test  results  show  that 
estimating  the  explosion  moment  by  inverting  for  the  isotropic  component,  in  conjunction 
with  a  double-couple  component,  should  yield  the  most  accurate  measurements. 

The  error  in  moment  estimate  due  the  difference  between  the  modeled  source  depth  (600 
m)  and  actual  source  depth  (between  300  and  750  m)  is  minimal  (0.25  percent),  unless  the 
difference  in  depth  should  place  the  shot  in  a  medium  with  a  significantly  different  a2  IP 2 
ratio.  The  error  in  using  one  generic  centralized  NTS  location  to  calculate  the  distance  for 
all  station-path  synthetic  is  also  minimal,  being  at  the  most  0.5  percent  for  the  extreme  case 
of  short  regional  paths  ( da  400  km,  8d  =  20  km). 

Another  source  of  error  in  explosion  moment  determinations  is  caused  by  mis-modeling 
the  near-source  region,  in  particular  the  shot-point  medium  elastic  properties,  since,  by  eq. 
(3),  the  displacement  generated  by  an  explosion  is  inversely  proportional  to  a2/P2  which,  in 
turn,  is  related  to  Poisson’s  ratio  (<j).  Various  studies  have  determined  the  P-wave  velocity 
structure  for  the  various  NTS  sub-regions;  however,  the  S-wave  velocity  structure  is  far  less 
well  constrained.  Depending  on  the  shot  depth  and  particular  model,  a2/p2  varies  between 
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3.25  and  2.75  (3  for  a  Poisson  solid),  leading  to  potential  systematic  errors  on  the  order  of 
8  percent  (0.035  log  units). 

This  effect  is  significant  and  could  be  corrected  for  in  terms  of  sub-site  coupling  coef¬ 
ficients  or  calibration  constants.  In  a  global  nuclear  monitoring  environment,  this  factor 
may  only  be  ascertainable  by  indirect  means  (such  as  estimating  a2//32  from  the  ratio  of 
M DC(Ray)/Mi)c(Lov ))•  This  effect  is  investigated,  as  in  Woods  &  Harkrider  (1995),  by  examin¬ 
ing  the  moment-mj,  and  moment-yield  scaling  relationships  for  various  explosion  populations 
and  sub-populations  at  NTS.  Unfortunately  the  spectral  moment  data  set  is  significantly 
smaller  and  less  comprehensive  than  the  data  set  used  in  Woods  &  Harkrider  (1995),  since 
many  of  those  moments/Ms’s  were  determined  from  historical  analog  data.  However,  the 
available  data  set  does  contain  events  from  Rainier  Mesa  as  well  as  smaller  events  from  the 
other  sub-sites,  down  to  mb  =  4.9,  and  in  a  few  cases  down  to  the  mb  =  4.4  level,  whereas 
previous  (Stevens,  1986;  Given  &  Mellman,  1986)  studies  have  only  examined  mj,  =  5.5  and 
larger  events.  Hence  the  moment  scaling  relationships  can  be  better  constrained  for  Pahute 
Mesa  and  Yucca  Flats,  and  the  Rainier  Mesa  long-period  scaling  relationships  can  be  es¬ 
tablished  using  this  data  set.  Including  the  large-yield  explosions  from  the  previous  studies 
makes  the  data  set  even  more  comprehensive  in  scale  range.  By  comparing  the  relationships 
obtained  from  these  three  NTS  sub-sites,  one  can  get  a  handle  on  the  effect  of  near-source 
structure  and  shot  medium  effects  on  these  seeding  relationships. 

Source  structure  significantly  affects  absolute  surface  wave  amplitudes,  hence  surface  wave 
magnitudes  or  moments.  For  all  of  the  synthetic  seismograms  generated,  we  used  the  Stevens 
(1986)  NTS  source  elastic  structure,  which  was  also  used  by  Given  and  Mellman  (1986)  and  is 
basically  a  Pahute  Mesa  velocity  structure.  Since  there  are  known  shot  point  source  velocity 
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structure  variations  at  NTS,  we  estimated  the  error  in  log  moment  which  would  occur  by 
assuming  the  common  NTS  source  structure  instead  of  the  actual  structure.  Using  the  same 
approximations,  synthetic  Rayleigh  waves  were  generated  for  various  shot  point  structures 
reported  for  NTS.  These  were  compared  with  synthetics  using  the  Stevens  (1986)  NTS  source 
region  to  obtain  relative  moment  corrections.  The  NTS  sub-site  source  structures  and  the 
log  moment  corrections  are  given  in  Table  4.  The  correction  values  are  all  positive.  Thus, 
if  the  event  actually  occured  in  this  structure,  the  log  moment  obtained  using  the  Stevens’s 
NTS  source  structure  would  be  too  small  by  the  amount  given  in  the  tables.  On  the  average, 
the  corrections  for  Yucca  are  about  0.10  units  greater  than  Pahute.  This  implies  that  the  log 
explosion  moment  obtained  for  these  Yucca  structures  using  the  NTS  source  region  is  0.10 
units  less  than  the  Pahute  structures.  This  does  not  explain  the  entire  difference  observed 
by  Given  and  Mellman  but  it  is  in  the  right  direction.  It  is  interesting  to  note  that  South 
Yucca  model  yields  moment  values  closer  to  those  of  the  Pahute  models  than  the  North 
Yucca  model. 

In  order  to  reduce  the  effect  of  differing  shot  point  velocity  ratios,  Stevens  (1986)  suggested 
a  new  explosion  moment,  M'0,  defined  by 

K  =  (24) 

For  a  shot  point  medium  with  Poisson’s  ratio  of  0.25  (a2//?2= 3),  the  value  of  the  moment  is 
unchanged.  If  we  had  used  the  same  value  of  M'0  for  each  source  model  instead  of  M0  ,  we 
would  have  obtained  log  moment  corrections  that  are  the  difference  in  value  between  the  two 
entries  for  each  source  structure.  The  variation  of  this  value  over  the  test  site  is  small  and  is 
a  measure  of  the  effect  of  the  vertical  velocity  structure  and  not  the  shot  point  medium.  At 
long  periods  the  difference  in  spectra  from  the  standard  NTS  structure  is  given  by  this  ratio 


Table  4:  Various  NTS  sub-area  models  and  their  log  moment  corrections 
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and  if  we  had  used  the  frequency  range  of  Given  and  Mellman  (1986)  instead  of  Stevens, 
who  uses  more  high  frequency  values  to  average  the  spectra,  we  would  have  obtained  an  even 
smaller  effect  using  M'0.  The  use  of  M'Q  for  these  source  structure  models  effectivly  removes 
the  effects  of  differing  local  shot  point  properties.  None  of  this  however  addresses  coupling 
or  how  yield  is  related  to  Mq  through  the  shot  point  material  properties. 

Thus,  by  numerical  simulations  using  a  variety  of  different  NTS  structures,  we  found  that 
for  the  frequencies  of  interest  and  sources  in  the  upper  6  kilometers,  the  primary  effect  was 
due  to  the  difference  in  shot  point  velocity  ratios.  The  size  of  the  effect  can  be  predicted 
extremely  well  from  their  explicit  presence  in  the  mixed-path  expression,  eq.  (3).  As  an 
example,  our  Green’s  functions  are  computed  for  an  explosive  source  at  a  depth  of  600 
meters.  In  the  NTS  (Stevens)  source  structure,  the  second  layer  starts  at  a  depth  of  500 
meters.  There  is  a  significant  difference  between  the  Poisson’s  ratio  of  the  surface  and 
second  layer  in  the  source  earth  structure.  The  log  difference  between  the  square  of  their 
compressional  to  shear  velocities  would  predict  from  eq.  (3)  Ms  difference  of  0.17.  The 
actual  difference  between  the  Ms  of  a  surface  explosion  and  our  Green’s  function  is  0.16  with 
the  near  surface  explosion  smaller  as  predicted. 

Results 

Several  variations  of  spectral  source  inversions  were  conducted.  First  a  weighted  average 
spectral  moment  was  calculated  assuming  only  a  scalar  moment  from  an  explosive  source. 
No  station  corrections  were  included  in  these  calculations.  Variance  reciprocals  were  used  as 
the  weighting  coefficients.  These  source  parameters  are  given  in  Table  5;  for  each  event  the 
upper  entry  is  the  source  parameter  and  the  lower  one  is  its  standard  deviation.  Moment 
values  are  in  units  of  1015  Newton-meters  ( N  —  m)  and  the  strike  angle  is  in  degrees.  In  the 


36 


second  column  the  upper  entry  is  the  number  of  Rayleigh  observations  and  the  lower  entry 
is  the  number  of  Love  wave  observations  used. 

Next,  the  isotropic  (Mj)  and  double-couple  (M#)  moment  components  and  strike  angle  (<f>) 
were  then  determined  using  the  least-squares  inversion  method  described  in  the  previous  sec¬ 
tion.  Separate  inversions  were  performed  for  the  Rayleigh  waves  (M*(ji),  Mdc(r),  and  <f>(R)), 
Love  waves  Mdc(l)  and  <f>(L))>  and  the  combined  data  (M,(r+l),  Mdc(r+l)>  and  4>(r+l))- 
The  values  and  standard  deviations  are  also  listed  in  Table  5.  Again  variance  reciprocals 
were  used  as  the  weighting  coefficients  in  the  inversion.  In  the  last  column  is  the  F  factor, 
defined  by  Toksoz  et  al.  (1965)  as: 


M#  a2  _  3  M# 
“  Mj  2/?2  ~  2  Mj 


(25) 


for  near-Poisson  solids  (a2//?2  ~  3). 

The  different  inversion  schemes,  with  a  few  exceptions,  all  yield  similar  explosion  moment 
estimates  for  a  given  event.  Also,  most  strike  angles  for  the  double-couple  source  are  in  the 
Northwest  quadrant  (280°  E  <  (f>  <  5°  E).  The  15°  -  20°  NW  strike  angle  inferred/observed 
by  other  studies  of  NTS  tectonic  release  lies  within  the  2cr  confidence  levels  for  most  events, 
implying  that  the  results  obtained  here  are  consistent  with  others’  observations.  In  general, 
for  a  given  event,  4\r)  and  4>{L)  correlate  reasonably  well  and  Mdc(R)  and  Mdc(L)  are  within 
a  factor  of  two  of  each  other.  Notable  exceptions  to  these  observations  are  PILEDRIVER 


(66153),  HURONLANDING-DIAMONDACE  (82267)  and  MIZZEN  (75154).  PILEDRIVER 
is  the  only  event  in  this  data  set  detonated  at  Climax  Stock,  a  granite  pluton  and  is  be¬ 
lieved  to  have  generated  tectonic  release  with  a  thrust-like  mechanism  (Stevens,  1986),  so 
the  inversion  scheme  employed  here  would  not  accurately  model  it.  HURONLANDING- 
DIAMONDACE  was  a  double  event,  which  is  a  possible  explanation  for  its  anomalous  radi- 


Table  5  (a):  Inverted  Explosion  and  Tectonic  Release  Source  Parameters 
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0 
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14 

1 
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Table  5:  Source  parameters  for  explosions  and  associated  tectonic  release. 


Table  5  (c):  Inverted  Explosion  and  Tectonic  Release  Source  Parameters 


Table  5  (d):  Inverted  Explosion  and  Tectonic  Release  Source  Parameters 


Table  5  (e):  Inverted  Explosion  and  Tectonic  Release  Source  Parameters 
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ation  pattern.  It  is  also  very  close  in  location  to  the  PILED  RIVER  shot,  so  that  the  tectonic 
release  faulting  mechanisms  may  be  similar.  The  Rayleigh  and  Love  wave  double-couple 
moments  vary  by  a  factor  of  seven  and  the  strikes  are  off  by  90  degrees  for  this  Rainier  Mesa 
event.  However,  the  Love  wave  source  solution  is  poorly  constrained  having  only  a  minimum 
number  of  observations.  It  may  be  that  the  tectonic  release  mechanism  for  this  event  may 
vary  significantly  from  the  assumed  vertical  strike-slip  one.  MIZZEN  belongs  to  subgroup 
of  Yucca  shots  north  of  latitude  37.09  N,  many  of  which  yield  double-couple  mechanisms 
with  strikes  near  NS0°E.  Given  &  Mellman  (1986)  found  an  offset  in  the  M0  vs.  mh  curves 
between  northern  and  southern  Yucca  Flat,  which  was  interpreted  as  a  difference  in  Pois¬ 
son  ratios  between  shot  sites  or  as  having  different  stress  regimes  which  would  give  rise  to 
differing  tectonic  release  mechanisms. 

In  most  cases,  the  double-couple  moment  determined  from  Love  wave  observations  is  within 
a  factor  of  two  of  the  moment  determined  from  Rayleigh  waves  alone.  The  relative  errors 
of  the  double-couple  moment  estimates  are  appreciably  larger  than  those  of  the  isotropic 
moment  estimates,  so  the  2a  confidence  levels  of  the  two  independent  double-couple  moment 

estimates  overlap  and  are  thus  compatable.  The  joint  inversion  results  in  significantly  better 
constrained  M#  estimates. 

Many  of  the  Rainier  Mesa  shots  were  relatively  small  (mb  «  5.0)  and  consequently  not 
enough  observations  were  available  to  invert  for  a  source,  so  it  is  hard  to  ascertain  the 
predominant  tectonic  release  mechanism  at  this  sub-site.  However,  barring  the  double 
event  HURONLANDING-DIAMONDACE,  all  events  with  sufficient  observations  for  inver¬ 
sion  yielded  strikes  appropriate  for  a  vertical  strike-slip  fault  with  the  strike  angle  in  the 
North  West  quadrant.  Patton  (1988)  speculates  that  block  thrust  or  rotation  may  be  the 


M_tr  (10*10  H-ra) 


Figure  Tectonic  release  moment  plotted  vs.  isotropic  moment. 
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prime  mechanism  at  Rainier  Mesa;  however,  the  results  here  imply  that  strike-slip  motion 
can  explain  the  fundamental-mode  surface  wave  radiation  pattern  observations  adequately. 

Figure  9  plots  the  estimated  tectonic  release  (from  combined  Rayleigh  and  Love  wave  inver¬ 
sion)  vs.  the  isotropic  moment.  Events  with  no  observed  Love  waves  or  too  few  observations 
to  invert  for  a  double-couple  source  are  treated  as  having  no  tectonic  release  (M#  =  0).  For 
all  events  Mj  >  M#.  All  Pahute  events  show  evidence  of  tectonic  release.  Only  about  half 
of  the  Yucca  events  do  so.  Otherwise,  there  are  no  clear  trends  in  the  data.  The  amount  of 
tectonic  strain  release  does  not  appear  to  depend  on  whether  the  shot  takes  place  above  or 
below  the  water  table,  either. 

The  explosion  moments  and  errors  were  converted  to  log10  units  and  listed,  along  with 
time  domain  moments  (and  their  errors)  for  overlapping  events  determined  in  Woods  & 
Harkrider  (1995),  in  Table  6.  The  time-domain  moments  are  plotted  versus  the  spectral 
moments  in  Figure  10.  Spectral  domain  moment  errors  are  significantly  smaller  than  those 
of  time-domain  moments  by  over  0.1  magnitude  units  on  average,  resulting  in  a  reduction  in 
standard  deviation  from  41.3  to  12.2  percent.  However,  the  time-domain  log-moments  are 
in  close  agreement  with  the  scalar  spectral  log-moments;  yielding  estimates  within  15  to  30 
percent  of  one  another.  There  is  no  discernible  improvement  in  error  reduction  between  the 
different  types  of  spectral  moment  estimates.  Differences  between  constrained  explosion  (Mo) 
and  inverted  isotropic  (Mi)  log-moments  are  generally  less  than  0.1  log  units  (25  percent); 
often  the  difference  is  much  smaller. 

Figure  11a  plots  all  final  explosion  moments  from  this  study  vs.  teleseismic  mt  taken 
from,  or  calculated  using  the  method  of  Lilwall  &  Neary  (1985).  Vertical  error  bars  on 
the  data  points  are  the  estimated  error  (standard  deviation,  s.d.)  of  the  data  (note  that 


Figure  :  Time  domain  log  moments  regressed  against  spectral  domain  moments. 
In  the  top  figure  spectral  moments  were  determined  assuming  an  isotropic  source 
only,  while  in  the  bottom  figure  the  spectral  moments  were  determined  by  inverting 
for  an  isotropic  source  +  a  double-couple  source.  The  regressions  were  constrained 
to  a  slope  of  unity. 
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Figure  :  Mo  vs.  mj,  regression  curves  for  the  aggregate  NTS  population  and  for 
Rainier  Mesa. 
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for  most  solid  points  the  scaled  errors  are  smaller  than  the  symbols  and  hence  are  not 
viewable).  A  weighted  linear  least-squares  regression  was  performed  to  obtain  the  Mo:mb 
scaling  relationship.  Below  the  curve  is  the  regressed  scaling  relationship  and  below  that, 
in  parentheses,  is  the  standard  deviation  in  the  fit  of  the  data  to  the  curve.  Given  & 
Mellman  (1986)  formulate  the  M0:mb  (or  M/:mk)  scaling  law  as  a  function  of  M0,  so  for 
consistency  in  comparisons  the  same  is  done  here.  Here,  unlike  in  their  study,  the  slope 
is  not  set  to  0.9,  an  artificial  constraint,  needed  because  of  the  narrow  range  of  data  and 
based  on  assumptions  in  mb:yield  and  moment :yield  scaling  laws.  The  solid  line  is  the  fit  to 
the  data  and  the  dashed  lines  are  the  2cr  error  bars  for  this  line.  The  scaling  relationship 
is  well  constrained;  however,  the  standard  deviation  is  large,  being  0.2  log-moment  units 
(58  percent).  Separating  the  events  into  populations  detonated  below  and  above  the  water 
table  (Figures  lib  and  11c,  respectively)  does  not  reduce  the  scatter  significantly,  but  it 
does  change  the  slope  of  the  curve  some  (10-15  percent),  giving  an  indication  of  how  well 
the  regression  curves  are  constrained. 

Figure  lid  is  a  regression- curve  plot  of  Rainier  Mesa  explosions.  Although  the  curve  is 
not  well-constrained,  the  slope  is  close  in  value  to  that  of  the  other  curves.  Rainier  shots  are 
detonated  in  large  horizontal  tunnels,  unlike  other  shots  which  are  detonated  in  boreholes. 
So  there  is  some  question  as  to  the  complexity  of  Rainier  sources.  As  discussed  in  Woods 
&  Harkrider  (1995),  the  Ms- yield  relationship  is  similar  to  that  of  well-coupled  Pahute 
detonations  below  the  water  table  (BWT)  from  which  it  is  assumed  that  the  Rainier  shots 
can  be  treated  as  being  in  saturated  rock.  This  inference  is  substantiated  by  Taylor  (1983), 
who  reports  that  the  Mesa  sports  a  perched  aquifer.  Treating  the  Rainier  detonations  as 
having  been  below  the  water,  the  entire  NTS  data  set  is  again  separated  into  shots  occurring 


below  and  above  the  water  table  populations  in  Figures  lie  and  Ilf,  respectively.  The 
inclusion  or  exclusion  of  Rainier  events  does  not  appreciably  change  either  curve,  so  it  not 
clear  from  these  results  whether  Rainier  shots  are  detonated  in  water-saturated  rock  or  not. 

Figure  12a  is  the  moment -mb  regression  curve  for  Yucca  Flats  explosions.  Figures  12b  and 
12c  separate  the  Yucca  population  into  shots  below  and  above  the  water  table.  All  three 
of  these  curves  are  well-constrained  and  have  slopes  near  0.8,  similar  to  the  complete  NTS 
population  curves.  The  standard  deviations  of  these  curves  tire  slightly  smaller  them  that 
of  the  aggregate  NTS  data  curves,  indicating  that  the  sub-sites  do  have  relative  excitation 
levels  for  surface  waves  and  body  wave. 

Figures  12d,  12e  and  12f  are  analogous  curves  for  Pahute  Mesa.  These  curves  are  poorly 
constrained,  because  of  the  limited  data  set,  and  have  quite  different  slopes.  To  obtain 
useful  scaling  relationships  for  Pahute  Mesa,  more  data  is  necessary.  In  an  effort  to  better 
constrain  these  curves,  Rainier  shots  were  added  to  this  population  under  the  premise  that 
the  two  source  regions  have  similar  geology  (volcanic  mesas  of  tuff  and  rhyolite).  Although 
the  resultant  curves  are  well-constrained  (Figures  13a,  13b  and  13c),  their  slope  values  are 
quite  different  from  the  other  populations.  It  appears  that  the  Rainier  shots  have  a  larger  mb 
for  a  given  moment  than  Pahute  events.  This  point  will  be  further  examined  after  examining 
the  moment:log-yield  scaling  laws. 

In  order  to  better  constrain  the  77i(,:moment  scaling  relationships,  data  from  the  Stevens 
(1986)  and  Given  &  Mellman  (1986)  studies  were  included  after  correcting  for  their  modeled 
shot  points.  Figures  14  and  15  are  regression  curve  plots  for  the  combined  data  sets  analogous 
to  Figures  11  and  12.  The  main  improvement  is  to  better  constrain  the  Pahute  Mesa  curves 
(Figures  15d  and  15e)  which  have  slopes  in  line  with  the  regression  curves  for  the  other 
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Figure  3:  Mo  vs.  mj  regression  curves  for  Pahute  and  Rainier  events 
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Figure  .4:  Mo  vs.  mj  regression  curves  for  all  NTS  events  and  for  Rainier  events. 
Data  from  Stevens  (19S6)  and  Given  and  Mellman  (1986)  are  included. 
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Figure  :  Comparison  of  observed  and  synthetic  seismograms  for  paths  modeled 
in  this  study. 
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Figure  5:  Mo  vs.  mh  regression  curves  for  Yucca  and  Pahute  events.  Data  from 
Stevens  (1986)  and  Given  and  Mellman  (1986)  are  included. 


populations  (~  0.8).  The  Pahute  awt  curve  is  no  better  constrained  however. 

From  these  results  it  appears  that  the  mt,:  moment  scaling  relationship  curve  slope  is  near 
to  0.8  for  all  well-constrained  cases.  The  data  sets  were  then  regressed  again  with  a  fixed 
slope  of  0.8.  The  difference  in  off-sets,  relative  to  the  aggregate  NTS  curve,  are  given  in 
Table  7,  with  a  positive  value  meaning  that  for  a  given  moment,  the  mb  value  will  be  larger 
by  that  many  magnitude  units.  The  standard  deviation  of  the  data  fit  to  the  curve  is  also 
given,  so  that  goodness  of  fit  of  the  unconstrained  and  constrained  curves  can  be  compared. 
In  all  cases,  except  for  Pahute  events  alone,  this  difference  is  negligible.  For  Pahute  the  errors 
are  slightly  smaller  for  the  unconstrained  curve.  For  the  entire  NTS  data  set  the  fixed-slope 
scaling  relationship  is 

mb  =  0.8  x  Mo  —  7.0.  (26) 

Here  the  M0  can  be  interchangeably  used  with  Mj.  There  is  no  appreciable  off-set  between 
shots  detonated  above  or  below  the  water  table  at  an  individual  site  implying  that  P-wave  and 
Rayleigh  wave  coupling  effects  are  very  similar.  There  are  differences  between  sites,  however, 
the  largest  difference  being  between  Pahute  and  Rainier.  From  these  results  Pahute  Mesa 
will  have  the  smallest  mb  for  a  given  moment,  followed  by  Yucca  and  then  Rainier,  implying 
that  it  is  the  richest  in  long-period  energy  excitation.  The  scaling  curve  off-set  difference 
between  shots  above  and  below  the  water  table  are  negligible. 

In  order  to  determine  seismic  coupling  effects  and  moment-scaling  independent  of  another 
seismic  scale,  the  seismic  moments  were  also  regressed  against  yield.  Figures  16,  17  and  18 
are  log-moment  vs.  log-yield  regression  plots  analogous  to  the  M0  vs.  mb  plots  of  Figures  11, 
12  and  13,  respectively. 

For  the  complete  NTS  data  set  (Figure  16a)  the  curve  is  well-constrained,  with  a  slope  near 


Table  7:  Offsets  in  Mo:mb  regression  curves  for  NTS  subsites 
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unity  (1.04).  The  group  of  data  points  grouped  away  from  the  regression  curve  are  Yucca 
Flat  explosions  detonated  above  the  water  table  (AWT),  about  half  of  these  were  detonated 
in  alluvium  as  well.  Werth  et  al.  (1963),  in  a  previous  study  of  P-wave  coupling  found  that 
alluvium  had  the  lowest  seismic  coupling  coefficient  amongst  the  shot  medium  rock  types 
that  were  studied.  It  is  also  known  that  saturated  rock  tends  to  be  a  more  effective  medium 
for  seismic  coupling.  Compressional  velocity  tends  to  increase  with  decreasing  dry  porosity, 
as  the  bulk  modulus  increases,  so  by  equation  (3)  an  increase  in  a  lowers  the  observed 
displacement  for  a  given  moment,  Mj.  So  it  is  not  surprising  that  bombs  detonated  in 
such  muffling  material  will  give  lower  apparent  moments.  Splitting  the  NTS  population  into 
shots  detonated  in  water-saturated  or  unsaturated  rocks  leads  to  similar  sloping  curves,  with 
values  near  unity  (0.98,  1.03  and  1.00  for  curves  b,  c  and  e,  respectively)  but  with  different 
intercepts.  The  scaling  relationship  for  Rainier  events  (Figure  16d)  is  poorly  constrained; 
however,  a  slope  of  unity  is  within  the  2(7  confidence  level  of  the  curve. 

For  the  case  of  AWT  explosions,  excluding  Rainier  shots  (Figure  16f),  the  slope  of  the 
curve  is  considerably  larger  than  for  any  of  the  other  cases.  This  data  set  falls  into  two 
clusters  (above  and  below  log(Mj)  =  15.5),  the  upper  group  being  predominantly  Pahute 
events  and  the  lower  being  all  Yucca  events.  If  there  is  a  difference  in  coupling  between  the 
two  sites,  there  should  be  a  corresponding  offset  in  scaling  curves.  So  the  best-fitting  curve 
for  the  combined  data  set  will  give  a  very  skewed,  incorrect  curve  (slope=1.35). 

Separated  Yucca  and  Pahute  population  log-moment:log-yield  curves  are  shown  in  Figures 
17a  through  17f.  The  variance  is  significantly  smaller  for  the  Pahute  populations,  particularly 
for  BWT  events  (Figure  17e).  The  slope  for  this  curve  (0.86)  is  heavily-constrained  by  the 
one  low-yield  event  (REX),  where  as  the  Pahute  and  Pahute  AWT  curves  both  have  slopes 
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Figure  -8:  Mo  vs.  log-yield  for  combined  Yucca  Flat  and  Pahute  Mesa  data  sets. 
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Figure  .9:  Mo  vs.  log- yield  for  the  Rainier  and  complete  NTS  data  sets  which 
include  moment  data  from  Stevens  (1986)  and  Given  and  Mellman  (1986). 
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near  0.69  in  value.  In  fact,  the  best-constrained  curve  is  for  Yucca  BWT  shots  (Figure  17)b, 
slope=0.97).  The  Yucca  AWT  regression  curve  (slope=1.27)  is  not  well-constrained,  due  to 
the  highly  scattered  nature  of  the  data. 

Combining  the  Pahute  and  Rainier  data  leads  to  a  better  constrained  curve  for  BWT 
events  by  making  the  assumption  that  all  of  the  Rainier  shots  are  fired  in  saturated  rock. 
For  all  events  above  and  below  the  water  table,  the  slope  is  0.95  (Figure  18a).  For  the  case 
of  BWT  shots,  the  slope  is  1.06  (Figure  18b).  Combining  Pahute  AWT  shots  with  Rainier 
shots  gives  a  slope  of  0.90  (Figure  18c).  The  fact  that  the  Pahute  AWT  +  Rainier  slope  is 
furthest  from  unity  is  interpreted  to  mean  that  the  Rainier  events  couple  intermediate  to 
Pahute  shots  detonated  below  and  above  the  water  table. 

In  order  to  better  constrain  the  larger-yield  portion  of  the  regression  curves  (particularly 
for  Pahute  events),  the  other  studies’  moments  are  also  added  to  the  data  set,  bringing  the 
total  number  of  events  to  155  from  109.  Figure  19a  is  for  all  NTS  BWT  events;  the  slope  of 
the  scaling  relationship  is  1.01.  Including  Rainier  shots  leads  to  a  slope  of  1.03  (Figure  19b). 

None  of  the  AWT  curves  are  appreciably  changed  by  the  addition  of  only  a  few  points  (1 
to  3)  at  most.  The  most  improvement  in  constraining  the  data  is  for  those  data  sub-sets 
including  Pahute  events.  Figure  19c  is  for  all  Pahute  alone  and  the  slope,  1.09,  is  better 
constrained  than  that  of  the  curve  in  Figure  17d.  Including  Rainer  data  gives  a  slope  of  1.01 
(Figure  19d).  The  Pahute  BWT  data  and  corresponding  scaling  curve  are  shown  in  Figure 
19e  (slope=1.14),  unfortunately  the  lower  end  of  the  scale  is  still  poorly  constrained.  Figure 
19f  is  for  Yucca  BWT  data.  Its  slope  (0.96)  is  close  to  the  case  for  this  study’s  data  only 
(Figure  17b). 

From  these  results  the  estimated  average  slope  of  the  moment  :yield  scaling  curves  is  near 


unity,  with  a  2a  margin  of  error  of  0.05.  A  constrained  (fixed-slope)  least-squares  regression 
for  each  data  set  was  also  performed.  The  tabulated  results  for  this  study’s  data  set  and 
the  combined  data  set  are  in  Table  8  and  Table  9,  respectively.  The  first  column  of  each 
table  describes  the  data  set  sub-population.  The  second  column  is  the  off-set  relative  to 
the  complete  NTS  data  set  determined  in  this  study  (£D).  The  last  column  is  the  standard 
deviation  for  each  fix-slope  regression  curve,  which  are  generally  not  appreciably  larger  than 
for  the  free-slope  regression  curves.  The  standard  error,  a%,  defined  as 

=  7^’ 

with  a  being  the  standard  deviation  and  n  being  the  number  of  observations,  is  between  0.01 
and  0.04  for  £D  in  all  cases.  So  £D’s  larger  than  this  amount  are  statistically  significant. 

Using  either  data  set  (this  study’s  or  the  combined  data  set)  a  large  offset  is  seen  in  the 
log-moment :log-yi eld  intercept  between  shots  detonated  above  and  below  the  water  table. 
For  Pahute  this  offset  is  0.22  log  units  (a  difference  of  66  percent).  For  Yucca  the  offset  is 
0.41  log  units  (a  275  percent  difference).  The  difference  in  apparent  coupling  between  Pahute 
and  Yucca  shots  detonated  in  saturated  rock  is  only  0.08  log  units  or  20  percent.  The  Rainier 
scaling  curve  intercept  appears  to  be  intermediate  to  that  of  the  Pahute  BWT  and  Pahute 
AWT  curves.  Yucca  AWT  shots  have  the  lowest  apparent  coupling,  being  0.3  log  units  (200 
percent)  lower  than  that  of  the  aggregate  NTS  data  curve.  The  difference  in  offset  between 
composite  (BWT+AWT)  data  sets  is  small  for  all  subsites,  being  less  than  0.07  (17  percent) 
in  all  cases.  Hence,  site  source  effects  are  relatively  small;  however,  whether  or  not  the  shot 
occurs  in  saturated  material  is  a  significant  effect. 


Table  8:  Offsets  in  constrained  Mo:yield  regression  curves  for  NTS  sub-sites. 
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Table  9:  Offsets  in  constrained  Mo:yield  regression  curves  for  NTS  sub-sites  using 
the  combined  data  set. 
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Conclusion 

Spectral  moments  were  determined  for  109  NTS  explosions  using  regional  surface  wave 
data.  This  data  set  is  significantly  larger  and  more  comprehensive  than  previous  published 
studies  of  long-period  explosion  sources.  Whereas  these  studies  (Stevens,  1986  and  Given 
&  Mellman,  1986,  for  example)  are  confined  to  explosions  greater  than  ra&  =  5.4;  in  this 
study,  events  down  to  mb  =  4.9  were  consistently  observed  and  analyzed  (in  some  cases 
even  smaller).  With  this  extensive  data  set  it  is  possible  to  make  more  robust  inferences 
concerning  long-period  source  scaling  for  underground  nuclear  explosions. 

The  reduced  magnitude  threshold  is  due,  in  part,  to  the  inclusion  of  near-regional 
(d<700  km)  surface  wave  records  for  source  inversions.  Only  with  such  data  were  the  smallest 
events  (raj,  <  4.9)  measurable.  Systematic  errors  due  to  incorrectly  modeling  path  attenua¬ 
tion  effects  are  minimized  using  such  stations  as  well.  It  was  also  found  that  regional  digital 
long-period  stations,  specifically  RSTN  and  DWWSN  stations,  were  superior  for  observing 
the  smaller  events  than  were  the  earlier-model  analog  WWSSN  stations. 

Moment  estimates  assuming  a  pure  explosion  source  (i.e,,  the  average  of  the  station  mo¬ 
ments)  do  not  vary  significantly,  in  almost  all  cases,  from  estimates  determined  from  joint 
isotropic  +  double-couple  inversions.  Moment  standard  errors  (la)  are  also  comparable, 
being  on  average  approximately  7  percent  (0.03  log  units),  with  the  smaller  events  tend¬ 
ing  to  have  larger  errors.  This  is  a  significant  improvement  over  the  time  domain  moments 
determined  in  Woods  &  Harkrider(1995)  ( a  =  0.15).  Past  investigations  of  surface  wave  mag¬ 
nitudes  and  seismic  moments  did  not  achieve  this  level  of  accuracy,  even  when  secondary 
path  corrections  (or  more  exactly  station  corrections)  were  incorporated.  The  improvement 
in  accuracy  in  the  moment  estimates  found  in  this  study  are  attributed  to  using  near-regional 


stations  and  having  better  azimuthal  coverage. 

The  results  of  the  canonical  source  inversion  study  predict  errors  for  the  joint  source 
inversion  compatible  with  observations;  however,  the  observed  errors  for  the  constrained 
explosion  source  estimates  are  smaller  than  the  predicted  values.  A  possible  explanation 
for  this  fact  is  that  the  path  Green’s  functions  were  inverted  from  bomb  data  contaminated 
with  tectonic  release  energy,  so  that  azimuthal  radiation  effects  are  incorporated  into  them. 
Thus  azimuthal  amplitude  variations  will  be  smeared-out,  making  the  observed  non-isotropic 
component  appear  smaller. 

The  constrained  tectonic-release  moment  results  are  not  as  well  resolved,  with  standard 
deviations  varying  from  10  to  50  percent  (0.04  to  0.18  log  units).  Moment  estimates  deter¬ 
mined  separately  by  Rayleigh  wave  and  Love  wave  varied  by  as  much  as  a  factor  of  2.5.  The 
inferred  strike  angles  obtained  from  the  constrained  moment  inversion  generally  were  found 
to  confirm  earlier  studies  that  found  the  predominant  tectonic  release  mechanism  could  be 
modeled  as  a  near-vertical  strike-slip  double-couple  source  with  a  strike  near  N20°W. 

As  a  goal  of  this  study  was  to  re-evaluate  the  utility  of  surface  wave  measurements  to 
quantify  explosions  for  yield  and  discrimination  purposes,  particular  attention  was  paid  to 
error  analysis  of  the  moment  estimates.  To  this  end  it  is  important  to  determine  how  much 
near-source  effects  can  modify  apparent  moments.  Moment  scaling  relationships  for  the 
various  NTS  sub-sited  were  determined  using  this  study’s  analyzed  events  in  conjunction 
with  larger  yield  explosion  data  from  other  studies  in  order  to  better  constrain  their  curves. 
It  was  found  that  raj,  oc  0.8(±0.05)  x  log(M0),  while  M0  oc  y(i-°±°-°5) 

The  off-set  in  scaling  curve  intercept  between  sub-sites  is  at  most  0.07  log  units  (17  percent) 
for  the  super  data  set  and  only  0.02  log  units  (5  percent)  when  only  data  from  this  study 
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was  considered.  A  much  greater  difference  in  scaling  curve  intercepts  exists  between  events 
detonated  above  and  below  the  water  table,  being  between  0.22  log  units  for  Pahute  Mesa 
and  0.41  log  units  for  Yucca  Flats.  This  coupling  affect  appears  to  be  similar  for  P-wave 
amplitudes.  If  this  coupling  effect  cannot  be  accounted  for  in  practice,  then  predicted  yield 
estimates  could  off  by  a  factor  of  two.  Gupta  et  al.  (1989)  found  that  the  spectral  content  of 
P-waves  could  be  used  to  determine  whether  a  shot  occurred  above  or  below  the  water  table. 
In  an  actual  monitoring  environment  where  only  remote-sensing  may  be  possible,  employing 
this  technique  would  enable  one  to  correct  explosion  moment  estimates  for  the  shot  medium 
coupling,  thus  significantly  improving  surface  wave  moment  estimates,  as  well  as  body  wave 
magnitudes. 

The  improvement  in  estimating  the  long-period  source  spectrum  of  explosions  found  in  this 
study  makes  such  measurements  more  useful  for  seismic  nuclear  monitoring.  The  resolution 
attained  in  this  study  makes  surface  wave  moment  estimates  a  prime  candidate  as  compan¬ 
ion  measurement  to  Lg  magnitude  (Nuttli,  1986;  Patton,  1988).  The  magnitude 

threshold  for  Lg  is  lower  than  that  of  surface  wave  measurements,  being  near  mb  =  3.1  for 
comparable  distances;  however,  for  paths  in  which  the  Lg  wavetrain  is  ’’pinched  out,”  a  term 
for  which  is  ”Lg  blockage,”  other  methods  are  necessary. 
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